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Executive  Summary 

This  MUR  I  project  is  aimed  at  designing  a  set  of  new  energy-harvesting  and  storage 
materials  (EHSM)  and  their  integration  of  load-bearing  airborne  structures.  The  EHSM  that  were 
covered  are,  (i)  organic  solar  cells  based  on  two  designs  ,(a)  dye-sensitized  solar  cells  (DSSC) 
and  (b)  solid  state  semi-conductor  polymer  based  solar  cells(SSSPS),  (ii)  thermoelectric 
materials  and  modules,  (iii)  piezoelectric  materials  and  (iv)  electromagnetic  generator. 

Research  on  organic  solar  cells  is  still  in  its  infant  stage  and  the  power  conversion 
efficiency  (PCE)  of  those  is  still  lower  than  the  commercial  Si  based  solar  cells.  But  the  cost  and 
weight  of  the  organic  solar  cells  are  lower^  thus,  providing  a  good  potential  for  airborne  energy¬ 
harvesting  materials,  Jenekhe  designed  selected  SSSPS  with  PCE  being  close  to  5%  while  Taya 
and  Cao  groups  designed  DSSC  with  PCE  of  7-10%  for  liquid  base  electrolyte,  and  PCE  of  5-6% 
for  gel  electrolyte.  The  DSSC  designed  by  Taya  group  has  two  versions;  one  based  on  thin  glass 
substrate,  the  other  based  on  polymer  substrate,  the  PCE  of  the  former  with  gel  electrolyte  is  7%, 
while  that  based  on  PET  substrate  is  5%,  Cao  group  focused  on  the  designing  nanostructured 
nanoparticles  with  aim  of  increasing  PCE. 

Taya  group  designed  new  thermoeiectric  materials  based  on  low-cost  and  low-weight  TE 
materials,  i,e  iVlg2Si  and  MnSi  for  use  in  the  intermediate  temperature  range,  thus  higher  specific 
figure-of-merit  (ZT/p,  where  ZT  is  the  flgure-of-merit,  p  is  density)  while  they  used  BiTe  for 
lower  temperature  range  and  assembled  segmented  TE  module  whose  power  density  is  0,75 
W/cin^,  Ju  designed  thermal  switch  by  which  the  heat  flow  can  be  controlled  for  more  enhanced 
or  more  retarded  flow  control.  Feraboli  group  at  UW  studied  extensively  the  processing  and 
characterization  of  the  battery  and  thin  film  solar  ceils  integrated  in  composite  panels.  The 
results  of  the  Feraboli 's  group  are  complimentary  to  those  of  Hahn’s  group  work  at  UCLA, 

Hahn  group  studied  experimentally  the  durability  of  commercial  EHSM  embedded  in 
structural  composite  panels  which  are  subjected  cyclic  bending,  then  measured  the  residual 
performance  of  solar  cells  and  Li-ion  battery  film. 

Inman  designed  a  set  of  new  energy-harvesting  and  storage  systems  (EHSS)  integrated  in 
structural  composite  panels  while  the  main  emphasis  was  on  the  mechanical  energy  harvesting 
based  on  piezoelectric  and  also  magnetic  system  of  Faraday  principle. 

University  of  Colorado(CU)  team  made  multi-scale  modeling  and  design  optimization  of 
EHSM;  (i)  Electrochemical  modeling,(ii)  Nano-scale  heat  propagation  model,  (iii)  optimal 
design  of  piezoelectric  energy  harvesting  systems  and  (iv)  UAV  test  bed  design,  CU  team  also 
studied  experimentally  the  optimum  cathodic  material  of  Li-ion  battery. 

University  of  Illinois  at  Urbana-Champaign  made  micro-characterization  of  the  damages 
observed  in  the  integrated  EHSM  components. 

The  subjects  covered  by  this  MURI  team  are  extensive,  spanning  over  nanostructured 
airborne  energy-harvesting  and  storage  materials  and  their  devices,  their  integration  into  load- 
bearing  structures.  Due  to  this  broad-band  research  area  under  the  five-year  period,  we  strongly 
feel  that  several  more  research  subjects  remain  to  be  performed  in  the  near  future.  They  include 
(i)  integration  of  a  thermal  energy- harvester  into  a  UAV  combustion  chamber  and  (ii)  integration 
of  airborne  sensors,  actuators,  energy-harvestor  and  battery  systems  into  a  single  cohesive 
system  that  can  be  durable  under  themo-mechanical  use  loading. 
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Chapter  1:  Introduction 


This  annual  report  covers  the  progress  made  on  our  MURI  for  the  period  of  May  1,  2006  - 
December  31, 201 1,  and  covers  on  the  following  subject  areas. 

Task  A,  Design  Tools  for  Multifunctional  Materials  and  Structures 
(Dunn,  Maute  and  Yang) 

Task  B.  Energy  Harvesters  based  on  thermoelectrics  (Taya  and  Yang) 

Task  C.  Energy  harvester  based  on  polymer  solar  cells  (Jenekhe,  Taya,  Cao,  Xu) 

Task  D,  Integrated  EHSS  components  in  load-bearing  structures  (Ju,  Hahn,  Feraboli  and 

Chasiotis) 

Task  E*  Design  of  Energy  Storage  System  (Cao,  Taya,  Kuga,  Lee  and  Xu) 

Task  F.  Electric  System  Design  (Kuga) 

Task  G.  Multimode  Energy  Harvesting  for  ISR/MAV  Missions  (Inman) 

Table  1  shows  all  the  tasks  that  were  performed  by  the  above  PI  and  Co-Pls 


Group 

Names 

PI,  Co-Pls 

Year  1 

Year  2 

Year  3 

Year  4 

Year  5 

UW 

Taya 

Kuga 

Xu 

Cao 

Jenekhe 

Feraboli 

- - - 

cu 

Dunn 

Maute 

Yang 

Lee 

— 

VT 

Inman 

UI-UC 

Chaosiotis 

Table  L  The  periods  of  PI  and  Co-Pls  working  for  this  MURI  project 


The  semi  “detailed  statements  of  the  above  tasks  are  given  in  the  following  pages* 


Chapter  2:  Achievements  of  various  tasks 

Task  A:  Multi-scale  Modeling  and  Design  Optimization  of  EHSS  (Dunn,  Maute  and  Yang) 
Objective: 

Conventionally,  batteries  are  designed  to  provide  a  single  function:  to  store  energy.  The  battery 
cell  is  shielded  from  mechanical  loads  by  enclosure  in  a  stiff  containment  and  the  integration  of 
batteries  into  structural  components,  subsystems,  and  devices  is  mainly  driven  by  thermal 
management  requirements.  Integrating  batteries  in  the  load  bearing  strueture  of  small  air  vehicles 
is  a  promising  concept  to  reduce  the  weight  and  storage  volume  requirement  for  on-board  energy 
storage  systems.  The  goal  of  this  research  thrust  was  to  develop  an  in-depth  understanding  of 
mechanical  phenomena  in  Li-ion  batteries  across  multiple  length  scales.  This  insight  enables  the 
rigorous  design  of  multi-functional  battery  architectures,  using  formal  optimization  approaches. 
CU  team  studied  also  three  more  sub-tasks,  (i)  nano-scale  heat  propagation  model,  (ii)  optimal 
design  of  piezoelectric  energy  harvesting  system,  and  (iii)  UAV  test  bed  design. 


Al,  Introduction 

This  report  summarizes  the  research  conducted  by  the  co-PIs  Dunn,  Maute  and  Yang  under  the 
MURI-EHSS  grant  at  the  University  of  Colorado.  The  research  focused  on  the  development  of 
design  tools  for  advanced  energy  harvesting  and  storage  materials  and  the  development  of  a 
UAV  test  bed  which  allows  assessing  the  system-level  performance  of  the  EHSS  materials  and 
devices.  The  work  on  analysis  and  design  methods  for  advanced  EHS  materials  was  concerned 
with: 


Electrochemical  modeling 
Nano-scale  heat  propagation  model 

Optimal  Design  of  Piezoelectric  Energy  Harvesting  Systems 
UAV  Test  Bed  design 


A2.  Electrochemical  modeling 

Integrating  batteries  in  the  load  bearing  structures  requires  to  model  batteries  as  “material”.  To 
evaluate  the  electro-chemical  and  the  mechanical  performance  of  such  materials  we  have 
developed  a  computational  simulation  framework  that  couples  electro-chemical  and  mechanical 
models  and  applied  it  to  the  analysis  of  a  representative  unit-cell. 

Electro-chemical  model  and  mechanical  coupling 

We  consider  a  simple  layout  similar  to  a  typical  Li-ion  battery,  that  is,  a  separator,  or  electrolyte, 
and  a  cathode.  For  the  sake  of  simplicity  and  as  a  first  step,  we  ignore  the  effects  of  time  in  the 
model  and  solve  the  steady-state  equations.  Since  we  are  also  interested  in  the  mechanieal 
performance  of  the  system,  we  model  electrochemical  stresses  resulting  from  the  diffusion  of  the 
metal  Li  into  the  porous  cathode  structure.  This  is  a  one-way  coupling,  from  electrochemistry  to 
mechanical  performance. 
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In  an  electrolyte  part,  we  are  solving  the  following  non-linear  diffusion  equations: 

DFc  ,  , 

DV^c  +— — -0 

e  Ij 

DFc  ,  ^ 

‘  V^c  ,  +\:  V'O  .  =  0 

JYf  !■'  ‘  I.i 

where  clh  is  the  concentration  of  Li-ions  in  the  electrolyte,  (t'u+  is  the  electrical  potential  in  the 
electrolyte,  De  and  Ke  are  the  diffusivity  and  the  conductivity  of  the  electrolyte,  respectively.  F,  R, 
and  T  are  Faraday  constant,  universal  gas  constant,  and  temperature,  respectively. 

Within  the  cathode,  the  electric  conduction  is  decoupled  from  the  (neutral)  Li  diffusion: 

C  Lt 

=0 


where  D,  are  the  diffusivity  and  the  conductivity  of  the  cathode,  cli  is  the  concentration  of  the 
metallic  lithium  in  the  cathode,  and  Oq,  is  the  electrical  potential  in  the  cathode. 


“Global”  boundary  conditions,  that  is,  the  conditions  on  the  “outer  boundary”  of  the 
computational  domain  are  of  no  critical  importance,  and  are  chosen  rather  arbitrarily  so  far.  For 
example,  on  certain  boundaries  we  prescribe  Li+  concentrations,  potentials,  LiH-  inflows,  e- 
inflows,  or  periodic  boundary  conditions.  “Local,”  or  electrolyte/cathode  interface  conditions, 
are  much  more  important  since  they  determine  the  electrode  kinetics.  This  part  of  the  model  will 
be  refined  in  the  future.  The  interface  dynamics  is  modeled  using  the  Butler-Volmer  relations, 
which  impose  the  continuity  of  mass/charge  fluxes  across  the  interface;  that  is  no  sources  or 
sinks  on  the  boundary.  The  fluxes  depend  on  the  local  concentrations  of  species  near  the 
boundary,  and  some  empirical  material  properties.  Namely, 


/ 


( a  F  ^  (  a  F  \ 


where  io  is  the  exchange  current  density,  r\s  is  the  surface  overpotential,  and  tta  and  etc  are 
respectively  empirical  anodic  and  cathodic  constants.  In  our  model,  we  use  aa=cxc=0.5, 
r|s=constant.  The  appropriateness  of  the  choice  of  these  values  remains  to  be  verified.  The 
exchange  current  density  is  computed  by 


where  Cum  is  the  limiting  Li  concentration,  and  kr  is  the  reaction  rate  (5e-7[m/s]). 

Finally,  everywhere  in  the  domain  the  “usual”  linear  elasticity  equations  are  solved.  In  the 
cathode  part  they  are  coupled  to  the  concentration  of  Li  via  eiectrochemically  induced  strains: 
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where  is  the  electrochemicaily  induced  strain,  ;9is  a  material  dependent  Vegard  coefficient, 
and  Cp  is  the  strain-free  Li  concentration. 

Unit  cell  analysis  is  the  concentration  of  Li-ions  in  the  electrolyte,  <t*Li+ 

We  apply  the  above  model  to  the  analysis  of  a  3-D  unit  cell  and  consider  a  3D  cube  aligned  with 
coordinate  axes.  Periodic  boundary  conditions  are  imposed  in  y  and  z  directions.  The  cube  is 
divided  into  two  subdomains,  representing  the  electrolyte  and  cathode.  At  x=0  we  prescribe  the 
inflow  of  Li-ions  in  the  electrolyte  part  [thus,  Neumann  b.c.  on  CLi+  and  and  concentration 
of  metal  Li  and  electric  potential  in  the  cathode  [Dirichlet  b.c.  for  Cli  and  Oe-].  At  x=H  we 
prescribe  concentration  of  Li+  and  electric  potential  in  the  electrolyte  [Dirichlet  b.c.  for  CLi+and 
Oli+J,  and  concentration  of  metal  Li  and  electric  current  in  the  cathode  [Dirichlet  b.c.  for  criand 
Neumann  b.c.  for  The  considered  geometry  and  computational  results  are  displayed  below: 


Figure  Al.  Computational  domain.  Left:  cathode;  right:  electrolyte. 
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Figure  A2.  Concentration  of  Lithium,  Left:  cathode;  right:  electrolyte. 
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Figure  A3.  Electric  potentials.  Left:  cathode;  right:  electrolyte. 


Figure  A4.  Von  Mises  stress  distribution  in  the  electrode. 

We  developed  a  multi-scale  finite  element  model  of  lithium  batteries  to  study  electrochemical- 
mechanical  interaction  phenomena  at  macro  and  micro  scales.  We  developed  a  multi-scale 
approach  to  analyze  the  surface  kinetics  and  electrochemical-mechanical  phenomena  within  a 
single  spherical  particle  of  the  active  material.  Homogenization  techniques  relate  parameters  in 
the  micro-scale  particle  model  to  those  in  the  macro-scale  mode!  describing  the  lithium  ion 
transport,  electric  potentials  and  mechanical  response  based  on  porous  electrode  theory.  At  the 
macro-scale,  we  have  extended  Doyle  and  Newman’s  electrochemical  porous  electrode  model  to 
account  for  elastic  deformations.  At  the  micro-scale  we  have  accounted  for  differences  in  Li+  ion 
flux  into  the  particles  due  to  interfacial  surface  conditions  using  the  Butler-Volmer  equation  and 
for  surface  pressures  exerted  on  the  particles  by  the  composite  electrode  matrix.  We  have 
introduced  a  meso-scale  aggregate  model  to  relate  micro-  and  macro-scale  mechanical  effects. 
This  model  was  implemented  into  a  finite  element  framework  and  used  to  study  the  influence  of 
loading  and  design  parameters  on  the  electrochemical  and  mechanical  response  of  the  battery 
cell.  The  macro-  and  micro-scale  models  are  discretized  in  time  by  an  implicit  Euler  backward 
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scheme  and  in  space  by  a  Galerkin  finite  element  method.  The  nonlinear  macro-  and  micro-scale 
subproblems  are  solved  by  Newton’s  method  using  analytically  derived  Jacobians. 


Initia[  design 


Optimized  design 


Li*  conoentralion  Lt  conoentnatton 

in  electrolyte  in  partides 


Fig.  A5  Maximizing  the  useable  energy 
densUy  of  a  Li-ion  cell  for  given  discharge 
rate:  The  top  figure  shows  the  layout  of  the 
initial  and  optimally  graded  electrodes, 
varying  the  porosity,  Le.  electrolyte  volume 
fraction,  and  the  size  of  the  active  particles 
across  anode  and  cathode.  The  plots  at  the 
bottom  show  the  Li-ion  and  Li  concentrations 
in  the  electrolyte  and  active  particles.  This 
example  is  for  a  rather  high  discharge  rate 
which  leads  to  a  poor  Li  intercalation  into  the 
active  particles  for  the  uniform  porosity  and 
size  design.  The  optimized  layout  allows  for  a 
significantly  better  utilization  of  the  active 
material. 


To  improve  the  energy  and  power  density  of  the  battery  cell  we  integrated  the  above  multi-scale 
model  into  a  computational  framework  for  optimizing  the  layout  of  electrodes.  A  key  component 
of  the  optimization  methodology  was  the  formulation  of  the  adjoint  sensitivity  equations  of  the 
multi-scale  battery  model.  The  efficient  solution  of  the  adjoint  equations  relies  on  the 
decomposition  of  the  multi-scale  problem  into  multiple,  computationally  small  problems 
associated  with  the  individual  realizations  of  the  micro-scale  model.  This  decomposition  method 
was  shown  to  significantly  reduce  the  computational  time  needed  for  sensitivity  analysis  versus 
numerical  finite  differencing.  The  potential  of  the  proposed  optimization  framework  was 
illustrated  with  numerical  problems  involving  both  macro-  and  micro-scale  performance  criteria 
and  design  variables.  The  usable  capacity  of  a  lithium  ion  battery  cell  was  maximized  while 
limiting  the  stress  level  in  the  electrode  particles  through  manipulation  of  the  local  porosities  and 
particle  radii.  We  further  applied  to  this  optimization  framework  to  tailor  the  electrode  layout  to 
specific  energy  versus  power  requirements. 

Accomplishitients/New  findings 

We  expanded  multi-scale  battery  models  to  account  for  external  mechanical  loading;  we 
developed  a  mathematically  rigorous  computational  optimization  approach  to  systematically 
design  the  layout  of  electrodes  in  multi-functional  batteries. 

We  found  that  the  load  bearing  capacity  of  Li-ion  batteries  is  negligible  if  liquid  electrolytes  are 
used  but  can  be  significantly  improved  using  an  all  solid-state  architecture;  we  found  that  the 
stresses  generated  through  intercalation  of  Li  into  active  materials  dominates  over  the  ones 
generated  by  reasonable  levels  of  external  loading;  we  found  that  functionally  grading  electrodes 
is  an  efficient  approach  to  optimize  both  the  power  and  energy  density  of  electrodes  and  to  tailor 
the  power-density  behavior  to  particular  application  scenarios. 

A3,  Nano-scale  heat  propagation  model 
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At  sub-micro-  and  nano-scales  the  Fourier’s  law  of  heat  conduction  is  no  longer  valid.  This  is  not 
especially  surprising  considering  the  fact  that  the  Fourier  relation  between  the  temperature  and 
heat  flux  is  empirical  and  is  not  based  on  any  sound  physical  principles.  Therefore,  we  resort  to 
modeling  heat  conduction  using  statistical  mechanics  principles,  in  which  we  consider  statistical 
properties  of  heat  carriers,  such  as  the  electrons  and  phonons.  As  a  first  step  in  this  direction,  we 
have  chosen  the  simplest  possible  kinetic  model,  namely,  the  gray  body  model.  Under  the  gray 
body  assumption,  the  only  heat  carriers  considered  are  the  acoustic  phonons. 

Nano-scale  optimal  design 

In  order  to  evaluate  the  potential  of  using  the  phonon  Boltzmann  equation  for  the  purpose  of  the 
nano-scale  optimal  design,  we  have  implemented  a  pilot  numerical  simulation  code  using  rapid 
prototyping  packages  MatLab  and  Comsol  Multiphysics.  Given  the  boundary  heat  fluxes,  we 
distribute  two  given  materials  in  a  prescribed  control  volume  in  order  to  minimize  a  desired 
performance  functional.  In  our  preliminary  study  we  have  chosen  to  maximize  the  temperature 
difference  between  two  points;  one  may  for  example  envision  attaching  a  thermo-electric 
material  at  these  points,  thus  maximizing  the  electric  output  of  the  system. 

Numerical  tests 

Importance  of  the  phase  space  resolution:  Many  works  on  nano-scale  heat  propagation  employ 
the  lattice  Boltzmann  method  (LBM)  as  an  underlying  numerical  scheme  for  solving  the 
underlying  governing  Boltzmann  equations.  The  lattice  Boltzmann  method  is  a  finite-difference 
type  scheme,  in  which  the  velocity  space  discretization  is  fixed  and  is  determined  by  the 
topology  of  the  considered  lattice.  In  our  numerical  experiments  we  found  out  that  proper 
discretization  of  the  velocity  space  is  essential  fur  obtaining  well  resolved  solutions,  and  that 
numerical  solutions  generated  by  LBM  generally  suffer  from  non-symmetries  due  to  the 
existence  of  “dedicated”  non-physical  space  directions  aligned  with  lattice  directions.  To 
illustrate  this,  we  have  performed  the  following  experiment.  A  point  heat  source  was  placed  in  a 
square  piece  of  Si.  We  solved  the  problem  on  a  sequence  of  finite  element  discretizations, 
keeping  the  spatial  discretization  fixed  while  increasing  the  resolution  of  the  velocity  space 
discretization  (see  Figures  A6,  A7).  One  can  clearly  see  gradually  improving  quality  of  the 
temperature  resolution. 
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Figure  A6.  Importance  of  the  proper  velocity  space  resolution;  temperature  fields. 


Figure  A7.  Importance  of  the  proper  velocity  space  resolution:  phase  space  view. 

Length-scale  effects:  To  justify  the  computational  costs  associated  with  solving  the  multi¬ 
dimensional  integro-differential  equations,  we  perform  the  following  numerical  experiment.  We 
model  a  square  piece  of  Si  with  a  square  SiO  inclusion,  heated  on  the  left-hand  side.  We  solve 
the  problem  for  three  different  length-scales  of  the  domain:  1  Ojim  x  I  Opm,  1  pm  x  1  pm,  and 
300nm  x  300nm.  Assuming  the  phonon  mean-free  path  in  Si  is  of  the  order  300nm,  we  span 
multiple  heat  propagation  regimes,  namely  diffusive  (Fourier),  transient,  and  ballistic.  We 
simulate  the  problem  using  the  standard  continuum  mechanics  steady  state  heat  conduction 
equations,  and  steady  state  gray  body  phonon  Boltzmann  equations.  As  expected,  continuum 
mechanics  equations  do  not  “see”  the  length-scale  changes,  whereas  the  Boltzmann  equation 
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starts  to  capture  small-scale  phenomena  (such  as  the  temperature  slip,  shading,  reflections),  see 
Figure  A8. 


Preliminary  optimization  results:  In  order  to  evaluate  the  potential  of  using  the  phonon 
Boltzmann  equation  for  the  purpose  of  the  nano-scale  optimal  design,  we  have  implemented  a 
pilot  numerical  simulation  code  using  rapid  prototyping  packages  MatLab  and  Comsol 
Multiphysics.  Given  the  heat  fluxes,  we  distribute  two  given  materials  (Si  and  Ge)  in  order  to 
minimize  a  desired  performance  functional.  In  our  preliminary  study  we  have  chosen  to 
maximize  the  temperature  difference  between  two  points;  one  may  for  example  envision 
attaching  a  thermo-electric  material  at  these  points,  thus  maximizing  the  electric  output  of  the 
system. 


Figure  A8.  Length-scale  effects:  square  piece  of  Si  with  SiO:  inclusion,  heated  on  the 
right  side.  Temperature  isolines  are  shown.  Upper  row:  continuum  mechanics  modeling; 
lower  row:  gray  body  phonon  Boltzmann  equation.  Domain  size:  left  column:  lOpm  x 
10pm;  center  column:  I  pm  x  I  pm;  right  column:  300nm  x  300nm. 
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Figure  A9.  Preliminary  optimization  results:  maximizing  the  temperature  difference 
between  two  points  in  the  computational  domain  with  heated  left  side  by  distributing  two 
materials  (Ge  and  Si).  Top  row:  10(^m  x  lOpm  domain;  bottom  row:  I  pm  x  1pm 
domain;  left  column:  optimal  material  distribution  (black:  Ge,  white:  Si),  center  column: 
optimal  temperature  fields;  right  column:  phase  space  view  of  phonon  energy  density. 
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A4.  Optimal  Design  of  Piezoelectric  Energy  Harvesting  Systems 


Objective: 

There  is  substantial  promise  in  the  use  of  piezoelectric  energy  harvesting  systems  to  power 
electronics  and  devices  through  the  conversion  of  ambient  vibrations  into  electrical  energy. 
Accurately  predicting  the  energy  harvesting  capability  in  these  systems  is  of  paramount 
importance  as  the  harvesting  structure  and  the  associated  circuit  must  be  “tuned”  in  concert  to 
the  signature  of  the  ambient  vibrations.  To  this  end,  our  objectives  were  to  develop  modeling, 
simulation,  and  computational  design  approaches  for  piezoelectric  energy  harvesting  systems, 
with  goals  of:  i)  exploiting  the  properties  and  layout  of  piezoelectric  energy  harvesting  elements 
embedded  within  or  mounted  on  a  composite  structure,  and  li)  incorporating  a  faithful 
representation  of  the  nonlinear  electrical  circuit. 

Brief  sum  mary  of  the  CU  team  progress  toward  the  above  objectives 

We  successfully  developed:  i)  a  computational  approach  to  analyze  and  design  piezoelectric 
energy  harvesting  systems  composed  of  layered  plates  and  shells  connected  to  an  electrical 
circuit;  ii)  an  approach  to  simulate  the  response  of  realistic  nonlinear  energy  harvesting  circuits 
that  can  be  coupled  with  a  computational  design  approach,  and  iii)  a  new  approach  to  harvest 
mechanical  energy  from  propagating  waves. 

Accomplish  men  ts/New  findings 

We  used  the  finite  element  method  to  model  the  coupled  electromechanics  of  the  piezo-electric 
harvesting  structure  and  a  lumped  parameter  model  for  the  dynamics  of  the  electrical  circuit.  We 
assumed  the  harvester  is  subjected  to  a  prescribed  harmonic  base  excitation  and  that  the 
structural  and  electrical  responses  are  linear.  We  developed  a  topology  optimization  approach  to 
design  the  layout  of  a  multilayer  structure  consisting  of  structural,  piezoelectric,  and  electrode 
layers,  as  well  as  the  electrical  circuit.  The  flexibility  of  our  formalism  admits  the  definition  of 
specific  system-level  objectives,  e,g.,  maximize  the  power  harvested,  in  an  algebraic  fashion.  In 
Rupp  et  aL  (2009)  we  describe  our  analysis  and  design  approaches  in  detail  and  present  examples 
that  demonstrate  the  versatility  of  our  approach  and  show  how  it  can  be  used  to  explore  general 
behavior  of  and  to  develop  overarching  design  principles  for  piezoelectric  energy  harvesting 
devices*  For  an  objective  of  maximizing  the  power  harvested,  we  investigated:  (i)  optimal 
designs  for  various  piezoelectric  to  substrate  thickness  ratios,  (ii)  the  effect  of  mass  loading  on 
optimal  design,  and  (iii)  the  sensitivity  of  designs  to  shape  variations*  Figure  A 10  shows  an 
example  of  a  basic  problem  setup  and  Fig*  A1 !  shows  results. 


Figure  AIO.  Design  problem  setup  of  a  square  plate  subject  to  harmonic  loading* 
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Fig.  A)  I .  Comparison  of  the  strain  pattern  for  a  plate  without  a  piezoelectric  layer  with  optimal 
distribution  of  piezoelectric  material  (red)  on  an  aluminum  substrate  (blue)  for  six  different  ratios 

of  piezoelectric  to  substrate  thickness 

We  also  developed  a  numerical  simulation  methodology  for  piezoelectric  harvesting  systems  that 
allows  for  high-fidelity  finite  element  models  of  the  harvester  electromechanics  as  well  as 
realistic  nonlinear  circuits  models;  a  typical  example  is  shown  in  Fig.  A 12.  Our  approach 
synthesizes  conventional  harmonic  finite  element  analysis  and  a  harmonic  balance  approach  that 
accounts  efficiently  for  the  nonlinear  circuit  behavior,  bypassing  the  need  for  computationally 
expensive  time  marching  schemes.  The  proposed  methodology  solves,  via  an  iterative  residual 
reduction  scheme,  for  the  coupled  electrical  output  and  structural  dynamics  of  the  piezoelectric 
harvesting  structure,  which  is  analyzed  in  the  frequency  domain.  The  proposed  method  is  well 
suited  for  integration  into  model-predictive  design  methods.  We  illustrated  the  capability  of  the 
proposed  approach  to  predict  previously  unreported  coupling  phenomena  between  the  nonlinear 
circuit  response  and  piezoelectric  structural  behavior.  The  increased  accuracy  in  the  circuit 
model  and  flexibility  of  this  method  over  current  models  outweighs  the  increased  numerical 
costs. 


Fig.  A 12.  Schematic  of  the  coupled  electromechanical  energy  harvesting  system  with  a  finite 
element  model  and  full-bridge  rectifier  with  voltage  (  and  current  ( 

interface  conditions 

Finally,  we  demonstrated  the  ability  to  harvest  energy  from  the  propagation  of  phononic  (elastic, 
acoustic)  waves  in  two-dimensional  piezoelectric  solids  by  spatially  patterning  the  polarization 
distribution.  We  simulated  the  wave  fields  by  the  finite  element  method  and  demonstrated  the 
ability  to  dynamically  alter  the  wave  propagation  by  switching  on/off  the  piezoelectric  behavior 
by  operating  the  electrodes  in  a  closed  or  open  circuit  configuration.  The  piezoelectric 
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polarization  patterns  were  non- intuitive  and  were  determined  by  topology  optimization.  We 
showed  the  interesting  response  of  optimally  patterned  phononic  devices  with  four  examples;  a 
filter,  a  waveguide,  an  energy  harvester,  and  a  wave  actuator.  Fig.  A13  shows  a  device  we 
designed  to  harvest  energy  from  elastic  waves  propagating  in  a  solid  via  the  optimal  layout  of 
piezoelectric  polarization  throughout  its  volume. 


Fig.  A13.  Piezoelectric  polarization  pattern  for  the  wave  (a)  harvester  and  (b)  actuator  with 
corresponding  distributions  of  the  (c-d)  Poynting  vector  magnitude 


A5.  Design  of  UAV  Test  Bed 

To  evaluate  the  performance  of  the  EHSS  technology  developed  in  this  MURI  project,  we  have 
developed  a  virtual  and  a  hardware  UAV  test-bed.  This  includes  a  detailed  high-fidelity 
aeroelastic  model  of  a  mid-size  UAV  which  has  been  developed  and  built  at  the  Research  Center 
for  Unmanned  Vehicles  at  the  University  of  Colorado  (see  A-Figure  A14).  The  numerical  model 
will  allow  the  MURI  team  to  evaluate  quickly  the  improvements  of  the  system-level 
performance  through  novel  EHS  materials  and  devices  and  guide  the  development  of  the 
fundamental  research.  It  will  further  allow  us  to  optimize  the  integration  of  one  or  multiple 
EHSS  devices  into  the  airframe.  The  UAV  hardware  will  mainly  serve  as  verification  platform  to 
test  the  EHSS  in-situ  and  as  technology  demonstrator. 

In  the  first  year  of  this  MURI  project  we  developed  a  large  aeroelastic  model  and  performed 
mesh  and  time-step  convergence  studies.  In  collaboration  with  Dan  Inman  we  are  currently  using 
the  virtual  UAV  test  bed  to  quantify  the  kinetic  energy  available  for  mechanical  harvesting  due 
to  representative  wind  variations. 
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half  model 


Mach  contours 
on  deformed  structure 


Figure  A 14.  Hardware  and  virtual  UAV  test  bed.  The  figures  in  the  upper  row  show  the  CAD 
model  of  the  HAVUC  UAV  developed  by  senior  aerospace  students  at  the  Research  Center 
for  Unmanned  Vehicles,  as  well  the  aircraft  on  the  ground  and  in  flight.  The  figures  in  the 
middle  row  show  the  finite  element  model  of  the  HAVUC  aircraft  developed  in  the  past  year 
under  this  MURI  grant.  The  FE  mode!  captures  accurately  the  CAD  geometry  and  contains 
shell  and  volume  elements.  The  half  model  has  59,124  nodes  leading  to  302,253  degrees  of 
freedom.  The  lower  row  shows  the  CFD  model  and  aeroelastic  simulation  results.  In  the  past 
year,  we  have  developed  a  3-D  inviscid  model  with  373,690  nodes  resulting  in  about  1.8- 10^ 
unknowns.  The  CFD  Euler  and  the  FE  structure  models  are  coupled  using  our  in-house 
aeroelastic  simulation  platform.  The  figure  in  the  lower  right  comer  shows  the  Mach 
contours  on  the  aeroelastically  deformed  vehicle  for  a  nominal  flight  speed  of  Ma=0.04  at  sea 
level. 


Task  B.  Energy  Harvesters  based  on  Thermoelectrics  (Taya,  Yang  and  Ju) 

Here  we  report  on  the  progress  in  the  following  three  subtasks 

Bl:  Design,  synthesis,  characterization  and  modeling  based  on  low-cost  high  specific  tlgure-of- 
merit  TE  materials  and  modules, 

B2:  Thermoelectric  characterization 
B3:  Thential  switch  design 
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BI:  Design,  synthesis,  characterization  and  modeiing  based  on  low-cost  high  specific  figurc- 
of-merit  TE  materials  and  modules  (Taya). 


Introduction 


Thermoelectric  phenomena  involve  the  direct  conversion  between  thermal  and  electrical  energy. 
Due  to  its  solid-state  technology,  it  has  unique  advantages  as  an  electrical  power  generator  such 
as  high  reliability,  no  vibration  and  quiet  operation.  With  these  advantages,  thermoelectric 
devices  have  been  used  in  a  wide  range  of  applications  from  waste  heat  recovery  to  refrigeration 
(Goldsmid,  1964).  The  efficiency  of  the  thermoelectric  materials  is  determined  by  the 
dimensionless  figure-of-merit  (Z7),  defined  as 

ZT=  T  (Bl) 


where,  S,  <j,  k,  and  T  are  the  Seebeck  coefficient,  the  electrical  conductivity,  thermal 
conductivity  which  is  a  sum  of  lattice  (k^/,)  and  carrier  component  (k^/),  and  the  absolute 
temperature  at  which  the  properties  are  measured,  respectively.  For  the  enhancement  of 
thermoelectric  performance,  higher  values  of  Seebeck  coefficient  and  electrical  conductivity, 
which  give  rise  to  maximized  power  factor  (5^ a),  and  lower  thermal  conductivity  are  required. 
However,  the  interdependence  of  these  physical  properties  makes  the  improvement  of  ZT 
difficult  because  an  increase  of  electrical  conductivity  usually  decrease  Seebeck  coefficient,  and 
increase  thermal  conductivity  by  Wiedemann-Franz  law. 


The  majority  of  research  on  thermoelectrics  are  aimed  to  increase  the  figure-of-merit  (ZT)  value 
without  considering  other  factors;  weight  of  TE  device,  cost-effectiveness,  environmental  issue, 
and  ease  of  mass  production.  The  weight  of  the  TE  modules  can  he  a  critical  parameter 
considering  specific  energy  efficiency  particularly  airborne  applications  where  specific  figure  of 
merit,  defined  as  ZT  divided  by  mass  density,  should  be  emphasized.  Recently,  thermal  energy 
harvesting  based  on  TE  generator  modules  has  been  applied  to  automobiles  (Saqr  et.  al.,  2008) 
and  unmanned  aerial  vehicles  (UAV)  or  micro-aero- vehicle  (MAV)  (Pogue  et.  al.,  2005),  see  Fig. 
Bl.  In  automobile  applications,  use  of  TE  modules  is  to  improve  fuel  efficiency  by  utilizing  the 
wasted  heat  dissipation  which  is  about  40%  in  gasoline  engine  system  (Yu  et.  al.,  2009). 
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Wing  span  Fuel  tank 


Fig.  B 1 .  Conceptual  view  of  combustion  engine  system  and  TE  generator  in  UAV  (Pogue  et.  al., 

2005) 


In  this  MURI  project  we  focused  on 

(i)  design,  synthesis,  characterization  and  validation  of  TE  materials  and  modules  based  on  low 
cost  and  high  specific  figure-or-merit,  and  (ii)  integration  of  such  TE  modules  into  UAV/MAV 
combustion  chamber. 

(I)  Synthesis  of  bulk  n-  and  p-type  TE  alloys  based  on  IVIgiSi  compounds 
Mg2Si  based  TE  alloys  are  light-weight,  toxic-free,  and  cost-effective  materials.  Popular  TE 
compounds  working  at  intermediate  temperature  up  to  500°C  are  based  on  Lead  Telluride  (PbTe) 
(Crane  et.  al.,  2009),  AgPb,„SbTe2+m  (LAST)  (D’Angelo  et.  al.,  2011),  Te/Ag/Ge/Sb  (TAGS) 
(Crane  et.  al.,  2009)  and  Co2Sb3  (El-Genk  et.  al.,  2006).  Even  though  these  TE  materials  provide 
higher  ZT  values,  they  have  large  mass  density  of  more  than  8.0  [g/cm‘^]  which  increases  the 
weight  of  TE  modules.  In  addition.  Lead  (Pb)  is  toxic  element,  and  Tellurium  compound  is 
classified  as  poison  materials.  Regarding  the  cost  effectiveness.  The  Clarke  index  of  Si  and  Mg 
has  higher  rank  of  2"''  and  8*,  respectively  (Lotka,  1923).  Te,  one  of  the  key  components  in  TE 
system,  is  rare  metal  element  and  the  common  materials  such  as  Ag  in  LAST  is  very  expensive, 
resulting  in  increasing  base  cost  of  final  module  device. 

Recently,  environmentally- friendly  energy  harvesting  system  attracts  a  great  deal  of  interest. 
Mg-Si  system  thermoelectric  (TE)  materials  are  one  of  the  promising  candidates  in  an 
intermediate  temperature  range  (Noda  et  al,  1992;  Aizawa  and  Song,  2006;  Riffel  and  Schilz, 
1998)  thanks  to  its  light  weight,  low-cost,  scalable  to  larger  size  and  non-toxic  materials,  best 
suited  for  future  air  vehicles.  The  specific  figure  of  merit  (ZT/density)  of  Mg2Si  is  one  of  the 
highest  values  as  shown  in  Table  Bl,  though  its  performance  is  not  higher  than  existing  bulk  TE 
materials  (Fig.  B2).  Our  goal  is  synthesis  of  high  specific  figure  of  merit  Mg2Si  and  MnSi 
thermoelectric  (TE)  materials  suitable  for  airborne  energy  harvesting  through  “nanocomposite”. 

Most  recently,  significant  advances  for  increasing  ZT  have  been  made  based  on  electron  and 
phonon  engineering  using  nanostructure  (Venkatasubramanian,  2001 ;  Halman  et  al,  2002;  Hsu  et 
al,  2004;  Sales  et  al,  1996).  The  extensive  research  in  thermoelectric  properties  of  superlattices 
(Venkatasubramanian,  2001;  Halman  et  al,  2002)  shows  that  ZT  enhancement  comes  mainly 
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from  thermal  conductivity  reduction  due  to  incoherent  phonon  scattering  at  the  interfaces  of 
periodic  superlattice  sites  white  concurrently  improving  or  maintaining  electron  performance. 


TE  Material 

ZT/density 

Mg:Si 

0.380 

SiGe 

0  331 

MnSi 

0.090 

0  178 

0.214 

Ago  j  gSbTejq 

0  210 

Table  Bl.  Specific  ZT  values  of  selected  bulk 
TE  materials. 


Tempefature,  T  /  K 

Figure  B2.  Figure  of  Merit  (ZT)  as 
a  function  of  temperature  for  some 
bulk  TE  materials 


Bl-1  Synthesis  and  characterization  of  TE  materials 
Processing  route 

In  order  to  obtain  high  ZT  value,  minimizing  thermal  conductivity  is  one  of  the  key  methods  for 
bulk  TE  alloy  system.  This  is  required  to  introduce  grain  size  reduction  up  to  nano  scale  to 
enhance  intensive  phonon  scattering  resulting  in  the  reduction  of  lattice  thermal  conductivity. 
However,  conventional  processing  technology  has  difficulty  in  fabricating  bulk  thermoelectric 
materials  with  very  fine  grains  because  of  grain  growth  during  solidification  process.  This  work 
enlightened  that  the  high  energy  mechanical  alloy  enables  large  quantities  of  nanoparticles,  and 
advanced  hot-pressing  technique  can  suppress  the  grain  growth  compared  to  traditional  sintering 
method.  The  idea  is  nanograins  produced  by  high  energy  mechanical  alloying  (planetary  ball 
milling)  and  rapid  sintering  technique  (spark  plasma  sintering)  can  introduce  a  high  density  of 
grain  boundary  interfaces,  which  scatter  phonons  more  effectively  than  charge  carrier,  and  result 
in  reduction  of  the  lattice  thermal  conductivity.  This  combination  process  of  mechanical  alloying 
and  Spark  Plasma  Sintering  (SPS),  see  Fig.  B3,  has  an  advantage  of  controlling  the  exact  molar 
ratio  of  Mg  to  Si  while  the  volatilization  of  Mg  via  conventional  approach  such  as  direct  melting 
method  is  inevitable  since  the  evaporation  temperature  of  Mg  is  lower  than  the  melting 
temperature  of  Si. 

Mechanical  alloying  (MA) 

In  this  study,  planetary  ball  mill,  see  Fig.  B4,  was  used,  which  generates  strong  centrifuge  forces 
to  create  impacts  of  high-energy  mechanical  alloy  actions  inside  the  container  by  introducing  the 
opposite  rotating  directions  of  grinding  jar  and  supporting  disc.  The  loading  and  unloading  of 
raw  materials  in  the  milling  jar  with  balls  were  performed  in  the  glove  box  system  at  Ar 
atmosphere  to  avoid  oxidation  of  the  starting  materials,  and  then  the  milling  jar  was  set  in  the 
planetary  ball  milling  equipment  at  dry  condition  followed  by  wet  grinding  milling  with  hexane 
as  the  process  control  agent  (PCA).  The  milling  liquid  was  used  because  oxygen-free  n-hexan 
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(C6Hi2)  has  the  additional  advantage  to  protect  the  material  from  oxidation.  After  finishing 
alloying  process,  the  milled  powders  were  transferred  to  the  glove  box  and  dried  in  vacuum 
environment  in  antechamber  of  the  glove  box,  and  then  finally  nanostructured  thermoelectric 
powder  was  obtained.  Fine  magnesium  particles  are  flammable  elements  when  it  is  exposed  to 
air,  and  it  is  even  more  dangerous  when  high  energy  mechanical  alloying  yields  nano  sized 
powders.  The  glove  box  system  is  introduced  to  protect  flammable  powers  from  oxidation  and 
other  contamination. 


Spark  Plasma  Sintering  (SPS) 

Spark  Plasma  Sintering  (SPS)  is  a  pressure  sintering  method  based  on  high  temperature  plasma 
(spark  plasma)  generated  by  electrical  discharge  ON-OFF  DC  pulse  energizing.  The  key  feature 
of  SPS  is  a  high  thermal  efficiency  due  to  the  direct  heating  of  the  graphite  mold  and  stacked 
powder  materials  by  the  large  spark  pulse  current.  Fig  B3  shows  the  configuration  of  our  SPS 
system  located  at  UW-CIMS  laboratory.  The  process  of  SPS  in  the  graphite  mold  can  be 
described  with  5  steps  such  as  spark  discharge,  a  local  high  temperature  state,  evaporation, 
melting,  and  necking,  see  Fig.  B4,  From  a  spark  discharge  at  the  contact  point  of  particles  a  local 
high  temperature  state  is  generated  and  then,  evaporation  and  melting  on  the  surface  of  powders 
occur,  and  necking  is  formed.  This  noble  sintering  technique  enables  nano  powders  to  be 
consolidated  with  a  high  densification  and  keep  still  nano  sized  grains  by  providing  rapid 
sintering  procedure,  in  the  process,  mechanically  processed  nanopowders  were  set  into  graphite 
die,  and  compressed  by  single-axis  pressurization  system  in  the  chamber  with  vacuum  state,  and 
then  sintering  followed  by  cooling  performed.  During  the  sintering  process  carbon  die  set  should 
be  inspected  frequently  because  nanopowder  can  be  melt  below  its  melting  point  at  a  certain 
condition  or  fracture  of  graphite  die  can  be  occurred  by  heavy  pressure  at  high  temperature. 


(a)  _  (b) 

Fig.  B3.  SPS  system  configuration:  (a)  overall  SPS  system,  and  (b)  a  graphite  mold  (shown  in 
dark  color)  in  the  chamber 
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>  Ar  atmosphere 

>  Oxygen  level:  <  0,2ppm 

>  Humidity:  <0.1  ppm 

>  Powder  storage 

>  MA  preparation 

>  Vacuum  drying 

>  Graphite  die  assemblyfor 
SPS  process 


1.  Process  setting  in  glove  box  system 


>  Rapid  sintering  process 

>  Minimizing  the  grain 
growth  resuiting  in 
intensive  phonon 
scattering 

>  Reduction  of  iattice 
thermai  conductivity 


2.  Alloy  process  (Planetary  Ball  Mill) 


3.  Bulk  disc  synthesis  (Spark  Piasma  Sintering) 


>  low  temperature  alioy  ing  process 

>  Dry/Wet  miiiing 

>  Nano  sized  grains  in  ailoyed  powder 


Fig  B4,  Powder  metallurgy  route  for  TE  material  synthesis 


Fig.  B4  illustrates  the  processing  steps:  (I)  the  preparation  of  ball  milling  process  in  the  glove 
box  (Innovative  Technology)  by  putting  raw  powers  and  grinding  ball  into  the  jar  and  sealing, 
(2)  mechanical  alloying  process  (Retsch,  PM  1 00)  in  Ar  environment  in  the  milling  jar,  (3) 
transferring  the  alloyed  powder  in  the  glove  and  preparing  graphite  die  assembly  for  SPS  process, 
and  (4)  running  SPS  (Sumitomo  Coal  Mining  Co.,  Ltd.,  Dr.  Sinter  1020S). 


Non-doped  Mg  based  TE  alloy 
n-type  of  Mg2Si 

Pure  Mg2Si  of  n-type  TE  alloy  was  synthesized  via  MA  followed  by  SPS  process.  The  initial  Mg 
and  Si  powders  were  mixed  in  alumina  milting  jar  with  alumina  ball  with  the  ball  to  powder 
weight  ratio  of  40  to  1,  and  the  jar  was  clamping  in  the  glove  box  in  Ar  atmosphere.  The  milling 
condition  was  100  rpm  for  24  hrs.  Fig.  B5  (a)  is  XRD  pattern  of  alloyed  compound  which  shows 
the  mixture  of  Mg  and  Si  together.  The  Mg-Si  compounds  were  transferred  to  glove  for  precise 
measurement  of  its  weight  and  control  of  Mg  and  Si  contents,  followed  by  encapsulation  of  the 
Mg-Si  powders.  The  compounds  were  pressurized  by  SPS  under  SOMPa  in  vacuum  chamber  and 
sintered  at  700C  for  3 min  holding  time  with  the  heating  rate  of  30C/min.  The  XRD  pattern  of 
solid  bulk  Mg2Si  shown  in  Fig.  B5  (b)  illustrates  pure  Mg2Si  alloy. 
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29  degree 

Fig.  B5,  XRD  pattern  of  n-type  Mg2Si:  (a)  Ball  milled  Mg-Si  composite  powder,  and  (b)  as- 

SPSed  Mg2Si  bulk  disk 


p-type  of  Mg2Ge 

As  to  Mg  based  p-type  alloy,  Mg2Ge  behaves  p-type  TE  materials.  Similar  to  synthesizing  n-type 
Mg2Si,  the  combination  of  MA  and  SPS  process  was  performed.  The  initial  Mg  and  Ge  powders 
were  mixed  in  alumina  milling  jar  with  alumina  bail  with  the  ball  to  powder  weight  ratio  of  40  to 
1,  The  mechanical  alloying  was  run  at  100  rpm  for  24  hrs.  Fig.  B6  (a)  is  XRD  peaks  of  the 
mixture  of  Mg  and  Ge.  The  graphite  die  setting  was  done  in  the  glove  box  to  prevent  oxidation 
and  contamination.  The  SPS  was  carried  out  under  50MPa  in  vacuum  and  sintered  at  500C  for 
3min  holding  time  with  the  heating  rate  of  30C/min.  Fig.  B6  (b)  shows  the  XRD  pattern  of  solid 
bulk  Mg2Ge. 


Fig.  B6.  XRD  pattern  of  p-type  Mg2Ge:  (a)  Ball  milled  Mg-Ge  composite  powder,  and  (b)  as- 

SPSed  Mg2Ge  bulk  disk 
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Both  n-type  of  MgjSi  and  p-type  of  MgiGe  were  synthesized  successfully  via  the  combination  of 
MA  and  SPS  technique.  However,  MA  process  generated  only  the  mixture  of  initial  raw  powders, 
not  alloy  compounds  after  milling  process,  see  Fig.  B6  (a)  which  means  MA  process  did  not 
generate  nanostructures  in  the  powders  as  well  as  grain  boundaries  to  increase  intensive  phonon 
scattering  to  minimize  thermal  conductivity.  Low  milling  energy  of  100  rpm  would  be  the  reason 
of  incomplete  alloying  during  MA. 

In  order  to  increase  the  milling  energy  (speed  and  time),  different  types  of  jar  and  ball  were 
carried  out.  Various  jar  and  balls  can  be  selected  such  as  alumina,  stainless  steel,  tungsten 
carbide,  agate  and  zirconium  oxide  depending  on  its  purpose  and  alloying  elements.  In  this  study 
of  Mg2Si  alloying,  three  types  of  alloying  medium  have  been  tried  for  the  suitable  selection. 
Fig.B7  is  XRD  peak  analysis  of  alloyed  powders  by  using  grinding  medium  of  alumina,  steel  and 
tungsten  carbide.  Alumina  and  steel  based  milling  jar  and  balls  are  widely  used  for  crushing  to 
powders  and  mechanical  alloying  process,  however,  iron  peak  was  observed  in  alloyed  powder  at 
300  rpm  for  8  hr,  Fig.  B7  (a),  and  AI2O3  peaks  were  also  detected  in  synthesized  powder  at  200 
rpm  for  24  hr,  Fig.B7  (b).  This  contamination  came  from  the  reaction  between  particles  and 
grinding  medium  due  to  the  high  energy  induced  by  high  rotating  speed  and  longer  time.  Fig.  B8 
is  (a)  SEM  image  and  (b)  EDAX  analysis  of  iron  contamination  of  Mg-Si-Ge  sintered  bulk 
sample  at  which  bright  spot  indicates  metal  contents  of  iron  coming  from  alloy  jar  and  balls. 


Fig.B7.  XRD  patterns  of  three  different  milling  mediums:  (a)  Stainless  steel,  (b)  alumina,  and  (c) 

tungsten  carbide 


The  EDAX  result  in  Fig.  B8  (b)  shows  the  amount  of  iron  contamination.  For  this  reason, 
tungsten  carbide  jar  and  balls  were  selected  thereafter  for  high  energy  mechanical  alloy  process. 
In  experiments  with  tungsten  carbide  jar  and  balls,  no  contamination  from  grinding  medium 
occurred,  which  is  evidenced  from  Fig.  B7  (c).  Therefore,  all  mechanical  alloying  processes  in 
this  study  were  performed  with  tungsten  carbide  medium  under  Ar  atmosphere. 
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Fig.  B8.  Fe  contamination  in  MA  of  stainless  steel  jar  and  balls:  (a)  SEM  image  of  sintered  Mg- 

Si-Ge,  (b)  EDA  X  data  of  the  surface 


Doped  Mg2Si  TE  alloy 

The  typical  TE  alloys  based  on  solid  solution  of  Mg2B’'^  compounds  are  Mg2Si  -  Mg2Ge,  MgiGe 
-  Mg2Sn,  and  Mg2Si  -  Mg2Sn.  Compared  to  Si  element,  Ge  and  Sn  have  higher  mass  density,  and 
Ge  is  very  expensive  element,  so  Mg2Ge  and  Mg2Sn  compounds  are  not  suitable  for  airborne  TE 
materials.  To  meet  the  requirements  of  light-weight,  non-toxic,  and  cost-effective  TE  materials, 
doped  Mg2Si  was  selected  and  studied. 

By  using  tungsten  carbide  milling  medium  for  higher  energy  MA  process,  complete  alloyed 
compounds  were  synthesized  after  MA.  Fig.  B9  is  the  XRD  peaks  of  alloyed  powder  with  350 
rpm  for  5h  (Fig.  B9a)  and  25hr  (Fig.  B9b)  milling,  in  which  complete  alloyed  Mg2Si  compounds 
are  directly  obtained  after  M A  with  25h  followed  by  SPS. 


Fig.  B9.  XRD  peaks  of  alloyed  powder  and  sintered  bulk  sample 
Bi  doped  n-type  Mg2Si 

Bismuth  (Bi)  was  selected  as  doping  element  because  the  density  of  Bi-doped  Mg2Si  alloy  is 
lower  than  that  of  other  solid  solutions.  To  examine  the  solid  solubility  of  Bi  in  Mg2Si  alloy  for 
the  maximum  performance,  MA  processes  with  various  Bi  contents  in  commercial  Mg2Si 
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powder  offered  by  Union  Material  Inc.  were  carried  out  under  350  rpm  -  25  hr  of  dry  MA 
followed  by  250  rpm  -  3  hr  of  wet  alloy,  and  were  sintered  by  SPS  method  with  80MPa— 600C- 
3  min  condition. 


Fig.  BIO  Microscope  images  of 
pulverized  powders:  (a)  SEM  image 
before  MA  (b)  SEM  image  after  MA 
(400rpm-l5hr)  (c)  TEM  image  after 
MA  (d)  TEM  image  of  nano  sized 
particle 


Fig.BlO  shows  the  microscopy  images  of  mechanically  alloyed  particles  of  Mg-Si.  Particle  size 
was  in  the  range  of  sub  micron  to  hundred  microns  (Fig.  BlOa)  before  alloying  process.  Under 
increasing  alloying  energy,  the  size  of  as-milled  particles  became  finer  less  than  20  micron  at 
400rpm  for  15hr  in  MA-2  (Fig.  BlOb,  c).  Fligh  resolution  transmission  electron  microscopy 
(HRTEM)  image  in  Fig.  Bll  (b)  shows  that  the  grain  boundaries  are  high  angled  randomly 
oriented  boundaries.  These  adjacent  grains  with  random  orientation  enable  more  phonon 
scattering  effectively,  which  gives  rise  to  the  reduction  of  lattice  thermal  conductivity. 

Therefore,  high  energy  mechanical  alloy  process  successfully  yields  nano  sized  alloyed  particles 
with  range  of  20nm  to  200nm  and  sub  boundary  interfaces  in  the  particles  with  sub  nanometers. 


Fig.  B1 1  HRTEM  images  of  «-type  Bi  doped  Mg2Si:  (a)  nano-particle  of  diameter  of  30  nm  (b) 
which  as  grain  boundaries  with  different  orientation 

The  electrical  conductivity  and  Seebeck  coefficient  were  measured  by  commercial  equipment 
(Ulvac,  ZEM-3),  and  thermal  conductivity  was  measured  by  a  laser  flash  system  (Ulvac,  TC- 
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9000)  from  room  temperature  to  600C.  The  results  of  the  measured  thermoelectric  properties  are 
shown  in  Fig.  B12.  As  larger  amount  of  Bi  contents  (larger  x  contents  in  figures),  electrical 
conductivity  was  increased  (Fig.  BI2a)  while  Seebeck  coefficient  was  decreased  (Fig.  B12b). 
This  is  because  the  carrier  mobility  was  increase  by  the  effect  of  Bi  dopant.  Thermal 
conductivity  are  little  influenced  by  Bi  (Fig.  B12c),  From  the  result  of  the  figure-of-merit  (Fig. 
B]2d),  3  at  %  Bi  doped  Mg2Si  gives  the  best  performance,  in  which  the  highest  ZT  was  obtained 
of  0.61  atbOOC. 
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Fig.  B12  Thermal  properties  of  Mg2Si:Bi,v:  (a)  Electrical  conductivity,  (b)  Seebeck  coefficient, 
(c)  Thermal  conductivity  and  (d)  The  figure-of-merit 

Since  Mg  is  ductile  particles,  there  can  be  a  loss  caused  by  mechanical  alloy  process.  In  order  to 
investigate  the  effect  of  excessive  Mg  contents  to  the  thermoelectric  properties,  mechanical  alloy 
with  three  different  amount  of  excessive  Mg  powder,  0  w/%,  2  wt%  and  4  w/%,  was  carried  out. 
In  this  case,  Mg  (99.9%),  Si  (99.9%),  and  Bi  (99.999%)  powder  were  used  as  the  starting 
materials  instead  of  pre-synthesized  commercial  Mg2Si  powder.  Fig.  B13  shows  the 
thermoelectric  properties  of  electrical  conductivity  (Fig.  B13a),  Seebeck  coefficient  (Fig.  B13b), 
thermal  conductivity  (Fig.  B13c)  and  the  figure-of-merit  (Fig.  B13d),  in  which  2%  surplus  Mg 
gave  the  highest  value  of  figure-of-merit,  0.44,  but  the  difference  is  less  than  3%  each  other.  This 
means  that  MA  process  can  control  exact  molar  ratio  by  avoiding  the  loss  of  Mg  contents  during 
alloying  process. 
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Fig.  B13  Thermal  properties  of  Mg2Sio  97Bioo3  with  excessive  Mg  powder:  (a)  Electrical 
conductivity,  (b)  Seebeck  coefficient,  (c)  Thermal  conductivity  and  (d)  The  figure-of-merit 

Bi-In  doped  n-type  Mg2Si 

In  order  to  improve  electrical  properties  of  Bi  doped  Mg2Si,  new  doping  element.  Indium,  was 
investigated.  In  addition,  to  verify  the  suppression  of  grain  growth  of  synthesized  alloys  via  MA 
and  SPS  processes,  two  specimens  were  prepared.  First,  Bi  doped  Mg2Si,  where  Mg2Si  powder 
was  offered  by  Union  Material  Inc.,  was  prepared  as  a  reference  alloys  via  direct  melting  process 
in  electric  furnace  followed  by  crushing  the  ingot  into  powders.  Then,  the  powder  is  consolidated 
via  SPS  process.  Second,  Bi-In  doped  lVIg2Si  was  synthesized  via  MA  followed  by  SPS. 
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Fig.  B14  TEM  images:  (a)  and  (b)  are  two  beam  bright  field  images,  and  (c)  diffraction  pattern 


Fig,  B15  TEM  images:  (a)  two  beam  bright  field  TEM  image,  (b)  diffraction  pattern  and  (c) 

STEM  image 

In  Fig.  B14,  two  beam  bright  field  TEM  images  show  grains  and  their  boundaries  with  grain  size 
of  more  than  2  micron.  This  large  grain  is  shown  in  Fig.  B14  (a)  and  diffraction  pattern  of  its 
grain  is  in  inlet  figure.  Fig.  BI4  (c).  One  of  the  smallest  grains  which  is  about  1  micron  is  shown 
in  Fig.  B14  (b).  Fig.  B15  shows  TEM  images  of  Mg-Si-Bi-In  alloy  and  its  diffraction  pattern.  Fig. 
B 1 5(a)  and  (b)  are  TEM  and  STEM  images,  respectively,  which  are  showing  the  size  of  grains. 
The  size  of  most  grains  is  smaller  than  500nm.  The  ring  pattern  of  the  diffraction  pattern  in  inlet 
figure  tells  different  orientation  of  grains  in  the  selected  area,  which  means  small  grains  of  this 
sample.  Based  on  TEM  observation,  it  is  proved  that  the  combination  of  MA  followed  by  SPS 
processing  can  generate  nano  size  powders  and  suppress  grain  growth  due  to  rapid  sintering 
process  via  SPS.  This  method  for  TE  material  synthesis  is  good  for  decreasing  thermal 
conductivity  resulting  from  the  reduction  of  lattice  thermal  conduction. 
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Fig,BI6  Thermoelectric  properties  of  Mg2Sio,97Bioo3  and  Mg2Sio,%Bio,o3lnooi:  (a)  electrical 
conductivity  (b)  Seebeck  coefficient  (c)  thermal  conductivity,  and  (d)  the  figure-of-merit 


Fig.  B16  shows  electrical  and  thermal  conductivities  and  Seebeck  coefficient  of  the  two  alloys. 
The  negative  sign  indicates  «-type.  While  Seebeck  coefficient  is  almost  identical  each  other  for 
the  measured  temperature  range,  electrical  conductivity  (Fig.  B16a)  of  the  bulk 
Mg2Sioy6Bio,o3lno,oi  is  increased  by  factor  of  2  from  room  temperature  to  around  350C  compared 
with  Mg2Sio  97Bio  o3f  and  the  difference  is  smaller  from  350C  to  600C.  The  enhancement  of  the 
electrical  conductivity  is  due  to  the  increase  of  mobility  by  doping  with  In.  The  thermal 
conductivity  of  Mg2Sio  ^Biooslno.oi  was  decreased  compared  with  Mg2Sio.97Bi0.03  due  to  longer 
alloying  time  with  higher  alloying  speed,  which  gives  rise  to  finer  grains.  Therefore,  the  figure- 
of-merit  of  Bi-ln  doped  Mg2Si  was  enhanced  up  to  0.68,  which  is  1 .5  times  larger  than  that  of 
Mg2Sio,97Bio.o3 alloy,  see  Fig.  BI6  (d). 


X  doped  MgjSi  (X=Na,  K,  Ag) 

Various  doping  elements  in  Mg2Si  system  were  studied  such  as  Na,  K  and  Ag.  MA  was  carried 
out  under  350  rpm  for  25  hr,  and  SPS  process  was  performed  at  800C  for  3min,  Fig,  B17  shows 
TE  properties  of  various  doped  Mg2Si.  Na  is  p-type  doping  element,  and  K  is  n-type  dopant 
while  Ag  doped  Mg2Si  has  a  transition  from  n-type  to  p-type  at  350C.  ZT  of  these  TE  alloys  are 
much  lower  than  that  of  Bi  doped  or  Bi-ln  doped  Mg2Si.  Thus  we  have  not  pursued  to  use  these 
alloys. 
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Fig,B!  7  Thermoelectric  properties  of  X  doped  IVfg2Si  (X=^Na,  K,  Ag) :  (a)  electrical  conductivity 
(b)  Seebeck  coefficient  (c)  thermal  conductivity,  and  (d)  the  figure-of-merit 


Bl-2  TE  generator  module 

Preparation  of  TE  materials  and  characterization 

Collaborating  with  Prof  lida  at  Tokyo  University  of  Science  for  Bi  doped  Mg^Si  (n-type),  and 
Prof  Itoh  at  Nagoya  University  for  higher  manganese  silicide  (HMS,  MnSi2-;f,Jc^0.250-0,273)  for 
p-type  materials,  segmented  TE  generator  was  fabricated  where  the  lower  temperature  TE  leg  of 
commercial  Bi-Te  based  TE  compounds  (Tellurex  Inc*)  was  used*  Fig.  B18  (a)  through  (d)  show 
temperature  dependent  TE  properties  of  n-type  Mg-Si  and  p-type  HMS  for  high  temperature 
segments,  and  n-  and  p-type  Bi-Te  for  low  temperature  segments.  TE  properties  of  all  TE 
materials  used  in  our  TE  module  are  given  in  Fig*  B18  where  the  high  temperature  TE  materials 
used  are  Mg2Si  (n-type)  and  HMS  (p-type),  and  low  temperature  TE  materials  are  n-type  BbTe^. 
xSe>t  and  p-type  Bi2..xSbxTe3» 
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Fig.  B18  Temperature  dependent  TE  properties  of  n-type  MgiSi  and  p-type  HMS  as  high 
temperature  segment,  and  n-type  B^Tea-xSCx  and  p-type  Bi2.xSbxTe3  as  low  temperature  segment: 
(a)  electrical  conductivity,  (b)  Seebeck  coefficient,  (c)  thermal  conductivity,  and  (d) 
Dimensionless  the  figure  of  merit,  ZT 


Module  assembly 

For  TE  module  assembly,  high  temperature  TE  legs  of  Mg2Si  and  HMS  are  connected  to  Cu 
electrode  by  brazing  technique,  and  segments  of  TE  materials  are  soldered  each  other.  Before 
brazing  and  soldering  process,  it  is  essential  that  Ti  {50nm)  and  Ag  (Ipm)  layers  are  coated  by  E- 
beam  evaporator  to  increase  the  interface  bonding.  Ti  is  active  element,  ideal  adhesion  metallic 
material  between  semiconductor  blocks  and  metal  coating  layer,  and  Ag  has  low  wetting  contact 
angle  which  helps  the  surface  stick  to  bonding  alloys.  Table  B2  shows  all  solder  materials  for 
bonding  these  interfaces.  All  bonding  process  is  carried  out  under  the  gas  flowing  (2  I/min) 
environment  of  Ar  balance  +  5%  hydrogen  to  prevent  joint  interfaces  from  oxidizing  and 
contaminating  during  high  temperature  process.  This  process  with  Ti  and  Ag  adhesive  layers 
ensures  lower  contact  resistance  and  strong  bonding.  After  the  assembly  of  tc  shape  TE  module, 
aluminum  nitride  (AIN)  substrate  is  attached  to  the  top  and  bottom  sides  of  pre-assembled  TE 
generator  by  using  ceramic  bond.  Since  AIN  is  thermally  conductive  and  electrically  insulating 
material,  it  transfers  heat  flux  from  heat  source  to  TE  segments,  while  preventing  from  short 
circuit  between  adjacent  n-p  pair  as  well  as  a  heat  source. 
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Alloys 

Soldering  materials 

Bonding  layer  (A-B) 

Processing  condition 

Alloy  1 

Ag-Cu-Sn-Zn 

Cu  -  {rtupi) 

680°C-5min 

Alloy  2 

Ag-Cd 

nrn2,P!-P2 
(n,,pi)  -  Cu 

425'’C-3min 

Alloy  3 

Sn-Ag 

fti-n},  P1-P2 
(«aP/)  -  Cu 

250°C-3min 

Alloy  4 

Pb-Ag-Sn-ln 

P1-P2 
(«/./>/)  -  Cu 

350°C  -SiTiin 

Table  B2  The  list  of  all  solder  materials  for  bonding  interfaces:  alloy  4  is  the  best  solder  among 
alloy  2  -4  for  low  temperature  junctions 

Contact  resistance 

Electrical  contact  resistance  plays  an  important  role  in  a  segmented  TE  module  generator.  High 
contact  resistance  between  TE  segments  as  well  as  at  Junctions  of  metal  electrodes  would 
decrease  the  generator’s  performance  by  reducing  output  power  generation.  Thus,  the 
minimization  of  contact  resistance  at  all  interfaces  is  highly  desired.  The  contact  resistance 
measurement  is  carried  out  based  on  the  four-point  probe  set  up  where  two  current  probes  are 
connected  to  both  ends  of  a  sample  to  generate  current  flow  through  the  sample.  One  voltage 
probe  is  located  at  the  end  of  the  Cu  electrode  while  the  other  voltage  probe  is  scanning  along  the 
sample  length  to  detect  voltage  drop  depending  on  distance  between  two  voltage  probes,  so  the 
electric  contact  resistance  is  obtained  from  ohm’s  law  based  on  measured  current  and  detected 
voltage  drops  across  the  bonding  interfaces  between  the  fixed  and  the  scanning  probe.  Fig.  B19 
shows  the  measured  resistance  associated  with  the  cross  section  area  of  TE  segments  depending 
on  the  distance  of  two  voltage  probes.  The  measured  contact  resistance  at  the  interface  of  Cu 
electrode  bonded  with  n-  and  p-type  TE  segments  are  4.8  [pfl  cm^j  and  3.2  [pfT  cm^]  at  hi^h 
temperature  side,  and  those  at  low  temperature  side  are  43.2  [pfl  cm^]  and  73.6  [pH  cm*"], 
respectively.  The  contact  resistance  across  the  interface  of  high  temperature  and  low  temperature 
TE  segments  are  75.2  [p£T  cm^]  for  n-type  and  94.4  [p£l  cm^]  for  p-type.  The  measured  average 
electrical  contact  resistance  across  the  junctions  is  49  [pST  cm^],  and  the  measured  electrical 
resistance  of  2-pair  segmented  module  is  25  [mQ  cm^]. 


Fig.  B19  Electrical  resistance  as  a  function  of  distance  through  segmented  legs  of  p-  and  n-type 

TE  column 
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Bl-3  Power  measurement 

The  assembled  TE  module  generator  is  tested  for  its  output  power  performance  under 
temperature  difference  (AT)  of  up  to  500°C.  A  heat  source  located  at  the  top  of  the  TE  module 
generator  is  made  of  two  cartridge  heaters  insulated  by  magnesium  oxide,  which  are  inserted  in 
metal  block.  It  gives  rise  to  hot  side  temperature  {Th)  on  top  of  TE  generator  module  after 
steady-state  temperature  is  reached,  and  heat  sink  is  Cu  solid  block  with  water  flow  channel 
inside  for  tab  water  circulation  to  maintain  the  bottom  side  of  the  TE  device  at  constant  cold 
temperature  (Tc).  Under  the  heater,  Fe  plate  is  inserted  to  simulate  combustion  chamber  wall, 
into  which  TE  generator  is  integrated  to  utilize  temperature  gap  between  chamber  temperature 
and  outside  air.  Then,  TE  module  device  is  stacked  between  Fe  plate  and  Cu  block  on  the  heat 
sink.  The  stacks  are  surrounded  with  2-inch  thick  insulation  wall  made  of  compressed  fiberboard 
to  ensure  1-D  heat  flow  vertically,  and  thermocouples  are  inserted  at  each  stacked  interface  as 
well  as  Cu  electrode  to  monitor  temperatures  through  heat  flow  direction. 

The  measurement  tests  are  carried  out  under  various  temperature  differences  {AT),  AT=  100°C  ~ 
500"C.  For  zl7’=500°C,  the  hot  side  temperature  of  top  AIN  substrate  of  the  module  is  set  to 
520‘’C  while  the  cold  side  temperature  of  bottom  AIN  substrate  is  kept  below  23®C  during  the 
whole  measurement  process.  The  measurement  at  each  step  of  temperature  difference  is  recorded 
after  waiting  for  the  steady-state  heat  flow  reached  typically,  in  22  hours  for  Jr=500°C.  The 
maximum  output  power  is  generated  at  a  half  of  open  circuit  voltage  under  the  impedance 
matching.  Fig.  B20  (a)  shows  the  output  power  generation  of  the  module  for  various  JFs.  At 
AT  of  498°C,  the  maximum  output  power  measured  is  0.48  [W]  at  the  load  current  of  3.4  [A] 
under  load  resistance  of  41  [mfi].  Since  the  net  area  of  TE  materials  in  2-pair  module  device  is 
0.64  [cm^],  the  power  density  of  the  TE  modules  (output  power  per  unit  area)  becomes  0.75 
[W/cm^].  Fig-  20  (b)  shows  the  power  generation  and  its  power  density  as  a  function  of 
temperature  difference,  which  grow  parabolical ly  with  increasing  temperature  gap. 


Fig.  B20  The  experimental  results  of  2-pair  segmented  TE  module  test:  (a)  Output  power 
generation  as  a  function  of  load  current  with  various -dT,  and  (b)  The  maximum  output  power 
generation  and  power  density  as  a  function  of  AT 
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Bl  -4  Modeling  of  segmented  TE  generator 

T 


Fig,  B21  General  segmented  TE  module  with  the  number  of  N„  for  n-type  leg  and  Np  for  p- 

type  column 


The  generalized  model  of  segmented  TE  module  is  shown  in  Fig.  B21,  and  it  is  made  of 
combined  TE  legs  consisting  of  the  number  of  n-  and  p-type  segments,  and  Np,  respectively. 
In  addition,  typical  TE  generator  is  made  of  connecting  N  pairs  electrically  in  series  and 
theimally  in  parallel. 

Open  circuit  voltage  of  N  pairs  is, 

Voc=N(S‘;  ~S^){T„~T,)  (B2) 

where,  and  5",^^  are  the  effective  Seebeck  coefficient  of  Np  and  Nn  number  of  TE  elements, 
respectively. 
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where,  m  is  the  ratio  of  load  resistance  to  module  resistance,  R[/R.  Ap  and  A„  are  cross  sectional 
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area  of  p-  and  n-type  TE  element,  respectively.  ,  p^^  and  are  electrical  resistivity  of 

f'  p-  and  n-type  element  and  contact  resistance  between  and  (i+lf'  of  p-  and  n-type  elements, 
respectively,  is  electric  resistance  of  lead  wires  measured  experimentally.  Generally,  the 
effect  of  Rw  diminishes  when  N  becomes  larger  since  internal  resistance  is  relatively  much  larger 
than  that  of  wires,  but  it  cannot  be  ignored  when  N  is  small.  The  load  voltage  can  be  obtained 
from  Eq.  (B2)  and  (B3), 


JJ 


+  R..  1  +  V, 


(B4) 


Thus,  output  power  is  obtained  as. 
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In  this  study,  it  is  noted  that  2  pairs  TE  module  is  assembled,  i.e.,  N=2  and  two  segments  are 
used  at  each  type,  i.e.,  Np^N„=2,  and  Ap^A„=A.  The  maximum  output  power  can  be  obtained 
under  /?£=/?,  i.e.,  m=\.  Therefore,  output  power  generation  is  expressed  as. 


1=1  1=1  1=1  1=1  ; 


The  efficiency  of  the  module  is  estimated  as. 


rj  = 


4 


(B7) 


Bl-4  ComparisoQ  of  experimental  data  with  predictions  by  model 

The  measurement  results  of  2-pair  module  are  compared  to  the  predicted  results  for  the 
validation.  Fig,  B22  (a)  shows  the  measured  and  predicted  I-V curves  under  d7’=498°C,  in  which 
the  solid  line  is  the  predicted  results  by  Eq.  (B4)  while  the  plotted  markers  are  the  measured  data 
points.  The  measured  open  circuit  voltage  [  Voc]  is  285  [mV]  and  the  predicted  value  is  295  [inV], 
which  is  in  good  agreement  within  about  3%  error.  The  measured  short  circuit  current  [/sc]  of 
about  6.7  [A]  is  close  to  the  predicted  value  of  6.65  [A]  within  1%  error.  This  good  agreement  is 
mainly  due  to  the  low  contact  resistance  at  bonding  interfaces.  The  measured  internal  resistance 
of  2  pair  module  at  room  temperature  is  23.1  [m£2],  and  calculated  value  is  22.3  [mQ],  This 
explains  that  reasonably  low  contact  resistance  at  the  bonding  interfaces  are  achieved,  and  the 
measured  current  is  in  a  good  agreement  with  theoretical  value  by  minimizing  the  current  drop. 
The  maximum  power  output  in  Fig.  B22  (a)  is  0.48  [W]  measured  at  the  load  current  of  3.4  [A], 
and  the  predicted  maximum  output  power  is  0.49  [W],  This  measured  value  and  the  theoretical 
calculation  for  output  power  generation  are  in  good  agreement  within  error  of  2%.  As  indicated 
in  Fig.  B22  (b),  the  generated  power  and  power  density  are  a  parabolic  function  of  JT  (^Tm-Tc) 
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in  an  agreement  with  the  predicted  power  equation,  Eq*  (B6).  The  maximum  efficiency  of  the 
module  as  a  function  of  zJT  is  illustrated  in  Fig,  B22  (c)  where  the  solid  line  is  the  prediction  by 
model  and  the  square  markers  are  the  measured  points.  The  maximum  module  efficiency  is  4,6% 
measured  at  zir=498°C,  and  the  predicted  efficiency  by  Eq,  (B7)  is  4,9%,  A  comparison  between 
the  measured  data  and  predictions  by  the  model  demonstrates  that  the  segmented  TE  module 
fabrication  and  its  measurement  process  are  carried  out  correctly  while  validating  the  model. 


Fig.  B22  Comparison  of  experimental  data  with  prediction  by  the  model:  (a)  FV  curve  and 
output  power  as  a  function  of  load  current,  (b)  the  maximum  output  power  and  power  density  as 
a  function  of  Al\  (c)  the  maximum  efficiency  of  the  module  as  a  function  of  AT.  Solid  lines  are 
analytical  modeling  and  marks  are  experimentally  measured  data,  and  the  maximum  output 
power  generation  as  a  function  of  N.  (Kim  et.  al.  2012) 

Bf-5  Study  of  TE  module  design 

In  this  project,  n  shaped  planar  TE  design  was  built  and  perfonnance  test  was  carried  out,  n 
shape  TE  module  is  a  conventional  design  concept  and  widely  used  for  TE  applications. 
However,  it  is  unstable  if  TE  legs  are  not  exact  parallel  each  other,  and  weak  to  shear  stress  at 
the  interfaces  between  TE  segment  and  metal  electrode  due  to  large  thermal  expansion  mismatch. 
Moreover,  planar  design  requires  a  large  number  of  metal  electrodes  and  big  pieces  of  ceramic 
plate  such  as  AIN  as  a  thermal  conductor  and  electrical  insulator.  These  parts  will  increase  the 
total  weight  of  TE  module  device.  Therefore,  use  of  '‘linear  TE  module’'  is  better  suited  for 
airborne  TE  generators. 

Table  B3  shows  four  possible  linear  TE  module  designs.  The  last  row  of  this  Table  B3  is  the 
module  based  figure  of  merit  by  accounting  for  cost,  density  and  energy  conversion  efficiency 
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which  can  be  termed  as  “specific  energy  conversion  per  cost  (SECC)  (defined  as 
efficiency/ (cost* density)”.  It  follows  from  Table  B3  that  linear  TE  module  design  (d)  gives  rise 
to  the  highest  power  density  while  design  (c)  provides  the  highest  SECC. 


Design  type 

Properties 

(a)  Linear  single 

(b)  Linear 
segmented 

(c)  Linear 
segmented 

(d)  Linear 
segmented 

Materials 

n:  Bi-Te 
p:  Bi-Te 

n:  Mg-Si  /  Bi-Tc 
p:  Si-Ge  /  Bi-Te 

n:  Mg-Si /Bi-Te 
p:  Mn-Si  /  Bi-Te 

n:  Mg-Si  /  Bi-Te 
p:  Zn-Sb  /  BLTe 

I’emperalure 

condition 

'rh=300C 

Tc-lOOC 

Tli=500C 

Tin=300C 

Tc=100C 

Th=500C 

Tin=300C 

Tc=l00C 

Th=500C 

Tin=300C 

Tc=100C 

Power  [W] 

0,21 

0,51 

0,38 

0.57 

Power  density 
[W/cm^l 

0.66 

1.60 

i.l8 

t.78 

Specific  Power 
[W/g] 

0,28 

1.39 

0.85 

1.14 

Efficiency  [%] 

3.2 

4.2 

3.3 

5.3 

Material  cost 
[dollar/g] 

0.091 

0.195 

0.046 

0.056 

Specific  power 
[W/g/dollarll 

3.08 

7,12 

18.48 

10.18 

Table  B3  Comparison  of  various  TE  module  designs  and  their  performance  in  terms  of  power 
density,  energy  conversion  efficiency,  TE  module  figure  of  merit,  specific  efficiency,  material 
cost  of  TE  module,  and  the  specific  energy  conversion  per  cost  (SECC) 

B2-2  Thermoelectric  Characterizafioii(Yaiig) 

We  have  built  setups  for  temperature-dependent  thermoelectric  characterization  of 
nanocomposites.  The  sample  setup,  shown  in  Fig.  B23(a),  is  housed  in  a  high  temperature 
furnace,  which  can  be  used  for  a  temperature  range  between  -70C  -llOOC.  The  electrical 
conductivity  and  Seebeck  coefficient  is  measured  using  the  four-wire/probe  measurement  and 
the  thermal  conductivity  is  measured  using  the  Angstrom  method. 

To  study  the  details  of  interface  effect  on  phonon  and  electron  transport,  we  started 
building  a  two-color  optical  pump-and-probe  system  from  Summer  2006,  shown  in  Fig.  B23(b). 
The  technique  could  also  probe  the  ultra- fast  non-equi librium  hot  electrons,  electron  phonon 
interaction  and  phonon  transport  to  gain  fundamental  understanding  of  the  electronic  and  thermal 
properties  at  the  material  interface.  Better  understanding  of  the  electron  and  phonon  dynamics  at 
the  interface  will  in  turn  help  in  synthesizing  better  thermoelectric  nanocomposites.  Figure 
B23(c)  shows  some  preliminary  signals  we  obtained  from  the  system. 

A  new  auto-assembly  process  has  been  proposed  at  CU  beyond  the  original  pick  and 
place  process  for  fabricating  miniature  thermoelectric  devices  which  we  termed  as  “universal 
fabrication  process”.  Such  a  unique  process  enables  fabrication  of  miniature  thermoelectric 
devices  using  ultrasonic  auto-assembly  for  mass  production.  We  have  designed  the  masks  and 
started  the  fabrication  process. 
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Figure  B23,  a)  sample  setup  for  temperature-dependent  thermoelectric  characterization  of 
thermoelectric  nanocomposites,  b)  two  color  optica!  pump-probe  system  for  studying  electron 
and  phonon  dynamics  at  interfaces,  c)  preliminary  signals  using  the  optical  pump-probe  system. 


B3  Design  of  thermal  switches  (Ju) 

Thermal  conductance  switches  based  on  three  physical  mechanisms  were  developed.  The  first 
switch  based  on  the  EWOD  (eiectrowetting  on  dielectric)  phenomenon  relies  on  electric-field 
induced  changes  in  the  shape  of  liquid  droplets  to  make  or  break  thermal  contacts  between  two 
planar  substrates  (Fig.  B24a).  Thermocapillary  flows  were  found  to  significantly  enhance  heat 
transfer  even  for  dielectric  liquids  with  low  thermal  conductivities.  The  second  thermal 
conductance  switch  is  based  on  reversible  morphological  transformation  of  microscale  liquid 
bridges  between  a  thin  continuous  liquid  film  (on  state)  or  discrete  droplets  to  form  or  break 
thermal  contacts  (Fig.  B24b).  It  can  achieve  on-state  thermal  eonductance  comparable  to  that  of 
direct  solid-solid  contacts  at  loading  pressures  several  orders  of  magnitudes  smaller  while 
avoiding  various  reliability  challenges  (cold  welding,  fracture/debris,  wear).  The  third  switch 
utilizes  the  preferential  alignments  and  chain  formation  of  magnetic  particles  dispersed  in  an 
inert  liquid  under  magnetic  fields.  This  design  eliminates  the  need  for  mechanically  moving 
mechanical  supports/substrates  but  does  not  provide  as  high  on-state  thermal  conductance  as  the 
second  switch  design. 

We  explored  an  alternative  integration  strategy  where  we  exploit  the  unique  mechanical 
characteristics  of  microscale  liquid  elements  to  enable  more  versatile  multifunctional  integration 
of  brittle  components  into  mechanical  structures.  Our  basic  idea  is  to  mechanically  decouple 
load  bearing  structures  from  energy  harvesting/storage  or  other  functional  devices  by  using 
liquid-based  mechanical  elements  (Fig.  B25a).  Liquid  elements  offer  many  intriguing 
mechanical  characteristics.  They  may  break  (rupture)  but  they  can  be  '‘perfectly”  healed.  They 
can  withstand  significant  cyclic  deformations  without  suffering  from  fatigue  and  also  be 
designed  to  change  shapes  and  slip,  merge,  or  split  on  solid  surfaces  in  response  to  external 
stimulations.  We  developed  a  rigorous  modeling  approach  based  on  the  surface  energy 
minimization  and  the  virtual  work  mode!  to  predict  the  shapes  and  mechanical  characteristics  of 
microscale  liquid  elements  (Fig.  B25b)  and  experimentally  validate  it.  We  experimentally 
demonstrated  the  feasibility  of  such  flexible  mechanical  interfaces.  A  silicon  chip  (thickness  of 
500  pm)  was  mounted  on  a  substrate  that  was  subjected  to  tension  and  bending.  The  liquid 
bridges  deformed  to  accommodate  the  various  mechanical  loading  while  effectively  protecting 
the  brittle  silicon  chip  (Fig,  B25c). 
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Fig.  B24:  Thermal  switches  based  on  (a)  EWOD  (b)  surface-tension  driven  reversible 
morphological  transformation  (c)  magnetorheological  instability. 
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Fig.  B25:  A  conceptual  illustration  of  the  mechanical  deformation  of  a  functional  component 
(red)  integrated  onto  a  load  bearing  structure  (black)  using  liquid  bridges.  Liquid  bridges  (blue) 
deform  or  rupture  in  response  to  various  mechanical  loadings,  (b)  Comparison  bet\veen  the 
experimentally  obtained  and  predicted  shapes  of  a  liquid  bridge  mounted  on  a  substrate  under 
bending,  (c)  Optical  images  of  the  integrated  structures  illustrating  intact  silicon  chips  mounted 
on  PDMS  substrates  under  three  loading  conditions. 


Task  C:  Energy  Harvesters  based  on  organic  solar  cells 
Here  we  report  on  three  subtasks 

Cl:  Solid  state  semi-conductor  polymer  based  solar  cells  (Jenekhe) 
C2:  Dye-sensitized  solar  cells  (DSSC)  based  on  ZnO  (Cao) 

C3:  Design  of  airborne  DSSC  (Taya) 


Cl:  Solid  state  semi-conductor  polymer  based  solar  cells  (Jenekhe) 

Cl  a.  Objective 

The  objective  of  this  task  of  the  EHSS  MURl  effort  was  to  develop  high-efficiency  polymer 
or  hybrid  organic-inorganic  solar  cells  for  integration  into  aero  vehicles.  Towards  this  goal  we 
performed  the  following  studies:  (i)  designed,  synthesized,  and  characterized  new  donor-acceptor 
(D-A)  conjugated  copolymer  semiconductors  that  combine  broad  absorption  spectra  with  high 
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charge  carrier  mobilities;  (ii)  fabricated  and  evaluated  bulk  heterojunction  (BHJ)  photovoltaic 
cells,  and  (iii)  attempted  to  integrate  the  organic  photovoltaic  (OPV)  cells  with  other  EHSS 
components.  Our  accomplishments  during  the  period  we  participated  in  the  EHSS  MURI  are 
summarized  below. 

Ob.  Accomplishments 

1.  Highly  Efficient  Solar  Cells  Based  on  Poly(3-butylthiophene)  Nanowires. 

We  studied  the  self-assembly,  charge  transport,  and  photovoltaic  properties  of  poly(3- 
butylthiophene)  (P3BT)  nanowires,  using  TEM  and  AFM  imaging,  field-effect  charge  transport, 
and  photovoltaic  measurements  (Figure  Cl).  We  have  shown  that  the  fullerene/P3BT-nanowire 
nanocomposites  exhibit  an 
electrically  bicontinuous  morphology 
in  P3BT  nanowire/PC7iBM  bulk 
heterojunction  (BHJ)  solar  cells.  Our 
results  demonstrate  the  achievement 
of  highly  efficient  solar  cells  from 
P3BT  nanowires  and  that  the 
photovoltaic  properties  of  poly(3- 
butylthiophene)  can  be  comparable  to 
the  much  studied  P3HT.  For  example, 
power  conversion  efficiencies  of 
3.3%  PCE  could  be  achieved  in  un¬ 
optimized  solar  cells.  This  discovery 
expands  the  scope  of  promising 
materials  and  architectures  for 
efficient  bulk  heterojunction  solar 
cells. 

We  have  also  studied  the  in- 
silu  self-assembly  of  P3BT  nanowires 
(NWs)  as  the  donor  component  in 
BHJ  solar  cells.  The  P3BT  NWs, 
self-assembled  in  the  presence  of  the 
fullerene  acceptor,  formed  an 
interconnected  network  in  the  blend 
films.  TEM  and  AFM  imaging 
revealed  P3BT  NWs  with  1 1-15  nm  width  and  several  micrometer  length,  fligh  hole  mobilities 
on  the  order  of  1  x  iO'^  cmVV  s  were  observed  in  the  P3BT  nanowires/PC61BM  nanocomposites 
by  both  SCLC  and  field-effect  transistors.  The  photovoltaic  properties,  especially  the  short- 
circuit  current  density,  fill  factor,  and  power  conversion  efficiency  were  found  to  depend 
strongly  on  blend  composition.  A  power  conversion  efficiency  of  2.52%  was  achieved  from  the 
in  situ  1:0.5  P3BT:  PCf,iBM  blend  in  ambient  air.  The  new  in  situ  self-assembly  approach 
provides  a  means  for  the  rational  control  of  the  film  morphology  in  bulk  heterojunction  polymer 
solar  ceils.  The  performance  of  the  P3BT  nanowire/PC6iBM  solar  cells  is  mainly  limited  by  the 
monomolecular  recombination  and  the  high  series  resistance,  which  should  be  addressed  in 
future  studies. 
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Figure  Cl.  (a)  Molecular  structure  of  poly(3-butylth]ophene) 
(P3BT)  and  schematic  illustration  of  nano  wire  network  of 
P3B  lYPCBM  composites,  (b)  TEM  image  of  P3BT  nanovvires, 
(c)  J-V  curA^es  of  P3BT-nw;PCBM  and  P3BT;PCBM  BHJ  solar 
cells,  (d)  J-V  curves  of  P3BT-nw:PCBM  and  P3HT:PCBM  solar 
cells. 
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2.  Regioregular  Poly(3-peiitylthiophene);  Synthesis,  Self-Assembly  of  Nanowires,  and 
Eftlcient  Photovoltaic  Ceils. 


Figure  C2.  (a)  Molecular  structure  of  poly(3-pentyl thiophene)  (P3P'r).  (b)  TEM  image  of  P3PT  nanowires,  (c) 
J-V  curves  of  P3PT  NW:PC7|BM  BHJ  solar  cells,  (d)  IPCE  spectra  of  P3PT  NWiPCtiBM  BHJ  solar  cell. 


We  have  synthesized  high  molecular  weight  (Mw  =  61800-77  000  g/mol,  PD1=1.4-1.5) 
regioregular  poly(3-pentyithiophene)  (P3PT)  and  investigated  its  electronic  and  optoelectronic 
properties  as  well  as  its  self-assembly  in  solution  (Figure  C2).  P3PT  combines  good  solubility  in 
organic  solvents  with  facile  solution-phase  self-assembly  into  crystalline  nanowires.  P3PT  films 
were  found  to  be  highly  crystalline,  exhibiting  a  lamellar  structure  and  molecular  packing  that 
are  characterized  by  15.1  A  interlayer  dioo  spacing  and  3.74  A  doio  7t-stacking  distance.  The 
P3PT  nanowires  (NWs)  assembled  from  dichlorobenzene  solution  had  a  width  of  16-17  nm  and 
aspect  ratios  of  70-465,  depending  on  the  solution  concentration.  P3PT  had  an  optical  band  gap 
of  1 .9  eV  and  thin  film  absorption  spectrum  that  are  almost  identical  to  those  of  the  well-known 
P3HT,  To  demonstrate  the  potential  of  P3PT  as  a  p-type  semiconductor  in  organic  electronics, 
we  have  fabricated  and  evaluated  bottom -contact  thin  film  field-effect  transistors  and  found 
average  and  maximum  mobility  of  holes  to  be  0.04  and  0.10  cm^/(Vs),  respectively,  for  the 
highest  molecular  weight  sample.  Bulk  heterojunction  solar  cells  based  on  either  P3PT/PC71 BM 
blend  thin  films  or  P3PT  nanowire/PC7iBM  thin  films  were  found  to  have  a  maximum  3.33- 
3.70%  power  conversion  efficiency  under  100  mW/cm^  AMI, 5  sunlight  illumination  in  ambient 
air.  The  similarity  of  the  electronic,  optical,  and  photovoltaic  properties  of  P3PT  to  those  of 
P3HT  when  combined  with  the  higher  melting  temperature  and  shorter  interchain  distances  of 
P3PT  make  regioregular  poly(3-pentylthiophene)  a  promising  p-type  semiconductor  for  organic 
electronics. 
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3.  Efficient  solar  cells  based  on  a  new  phthalimide-bascd  donor-acceptor  copolymer 
semiconductor. 

Bulk  heterojunction  solar  cells  based  on  a  new  small  bang  gap  donor- 
acceptor  copolymer ,  poly(A/-(dodecyl)'3,6-bis(4-dodecyloxythiophen-2-yl)phthalimide) 

(PhBTI2),  have  been  investigated  (Figure  C3).  PhBTI2/fullerene  blend  films  without  high 
temperature  (thermal)  annealing  were  shown  to  exhibit  ordered  nanoscale  morphology 
by  AFM  and  X-ray  diffraction  experiments.  A  space-charge- limited  current  hole  mobility  as  high 
as  4.0  X  10”“*  cm^A^s  was  measured  in  the  PhBTl  2/ful lerene  blend  films.  A  power  conversion 
efficiency  of  2.0%  was  achieved  from  the  1:1  PhBT!2:PC7iBM  blend  bulk  heterojunction 
solar  cells  with  a  current  density  of  6.43  mA/cm^,  an  open  circuit  voltage  of  0.56  V,  and  a  fill 
factor  of  0.55.  In  contrast,  thermally  (120  ®C)  annealed  PhBTl 2/fuI lerene  blend  films  were  found 
to  have  negligible  photovoltaic 
properties  due  to  severe  phase 
separation  which  is  likely  to  be  caused 
by  the  long  side  chains  and  resultant 
PCBM  diffusion.  The  performance 
of  phthalimide-based  donor- 

acceptor  copolymers  likely  has  large 
room  for  further  improvement  by 
modification  of  polymer  side  chain 
and  device  processing  methods  to 
control  film  morphology  as  well  as  the 
absorption  efficiency.  These  results 
demonstrate  that  phthalimide-based 
donor-acceptor  copolymers , 
exemplified  by  PhBTl  2,  are 
promising  low  band  gap  candidates 
for  developing  efficient  bulk 
heterojunction  solar  cells . 

4.  Organ  ometailic  Don  or- Acceptor  Polymer  Semiconductors  for  Solar  Cells. 

Eleven  Pt-bridged  organometallic  donoracceptor  conjugated  polymers,  incorporating 
various  electron  acceptors  have  been  synthesized  by  a  Sonogashira-type  coupling  polymerization 
and  characterized  (Figure  C4).  Most  of  the  polymers  had  weight-average  molecular  weights  (Mw) 
exceeding  10000,  while 
some  had  Mw  as  large  as 
21000-58000,  with 

polydispersity  index  around 
2.  By  varying  the  electron- 
accepting  strength  in  the 
donoi^acceptor  architecture, 
we  demonstrate  molecular 
engineering  of  the  optical 
absorption  bands,  electronic 
band  structure 

(HOMO/LUMO  levels). 
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Figure  C4.  (a)  Molecular  structures  of  organometallic  OA  polymer 
semiconductors,  (b)  J-V  curves  of  BT-BTPt:PC7|BM  BHJ  solar  cells. 


Figure  C3*  (a)  Molecular  structures  of  poiy(N-(dodecyi)“3,6-bis(4“ 
dodecyloxythiophen-2-yl)phthalimide)  (PhBTl  2),  and 

PC71BM.  (b)  The  current  density- voltage  characteristics  of 
PhBTl 2/PC71BM  (1:1)  and  PhBTl 2/PQri^M  (1:1)  blend  solar  cells 
with  and  without  annealing. 
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charge  transport,  and  photovoltaic  properties.  Although  the  absorption  spectra  of  thienopyrazine- 
containing  copolymers  (MTP-BTPt,  TP-BTPt,  and  PTP-BTPt)  were  broader  than  BT-BTPt  due 
to  the  stronger  effect  of  ICT,  the  improved  light  harvesting  did  not  translate  into  improved 
photovoltaic  properties.  The  observed  ionization  potentials  of-^E3  eV  facilitated  the 
injection  and  transport  of  holes  in  the  polymers,  which  showed  field-effect  hole  mobilities  of 
3.87  X  10“^  to  3.32  x  10“'^  cmV(V  s).  Electron  transport  was  not  observed  and  this  is  likely  a 
result  of  the  small  electron  affinities  (2.95-3.28  eV).  The  highest  photovoltaic  power  conversion 
efficiencies  were  observed  in  blends  of  the  polymers  with  PC71BM,  ranging  from  0.68-0.71% 
for  HPP-BTPt  and  BT-T50  to  2.41%  in  BT-BTPt.  Our  observed  2.41%  PCE  for  bulk 
heterojunction  solar  cells  based  on  BT-BTPt  as  the  donor  is  significantly  lower  than  the 
previously  reported  4.4-5.0%  PCE  but  is  much  closer  to  the  theoretically  estimated  2.2%  PCE 
for  this  polymer. 

5.  Don  or- Acceptor  Copolymer  Semiconductors  for  Solar  Cells. 

Several  D-A  copolymers  based  on  thieno[3,4- 
^Ipyrazine  as  the  acceptor  were  synthesized  by  Stille 
and  Suzuki  copolymerizations(Zhu  et  al,  2006).  In 
addition,  thiophene-quinoxaline  and  thiophene- 
pyridopyrazine  D-A  copolymers  were  also 
synthesized  by  Stille  copolymerization(Zhu  et  al, 

2008;  Kulkarni  et  al,  2008).  The  molecular  structures 
of  these  D-A  copolymer  semiconductors  are  shown 
in  Figure  C5.  As  expected,  Irong  intramolecular 
charge  transfer  between  the  electron  donor  and 
acceptor  moieties  in  the  copolymer  chains  resulted  in  broad  absorption  in  the  visible-near-IR 
region  (Figure  C6).  Some  of  these  copolymers,  e.g.  BTTP,  BTTP-F,  and  BTTP-P,  had  broad 
optical  absorption  bands  that  extend  into  the  near  infrared  region  with  absorption  maxima  at  667- 
810  nm,  and  optical  band  gaps  of  1.1-1.6  eV(Zhu  et  al,  2006).  The  other  copolymers  with 
quinoxaline  and  pyridopyrazine  acceptor  moieties  (PTHQx  and  PTDPPz)  also  have  broad  optical 
absorption  with  absorption  maxima  at  630-657  nm,  and  optical  band  gaps  of  1.75-1.83 

eV(Kulkarnl  et  al,  2008).  The  cyclic  voltammograms 
of  the  copolymers  showed  reversible  redox  peaks, 
implying  that  the  materials  may  be  capable  of 
ambi  polar  charge  transport. 

We  have  evaluated  the  charge  transport 
properties  of  all  the  D-A  copolymers  by  using  the 
field-effect  transistor  as  a  platform  and  gold  source 
and  drain  electrodes(Zhu  et  al,  2006;  Zhu  et  al,  2008; 
Kulkarni  et  at,  2008).  The  results  showed  that  BTl’P- 
F,  PTHQx  and  PTDPPz  have  good  field-effect 
mobilities  of  holes  with  values  of  1.6  xl0■^  8.5xl0'\ 
and  4.4x10'^  cm^/(Vs),  respectively.  However,  n- 
channel  field-effect  characteristics  and  thus  electron 
transport  were  not  observed  in  any  of  the  copolymers. 
This  means  that  the  materials  are  not  ambipolar 
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Figure  C6.  Absorption  spectra  of  the  D-A 
copolymers  as  thin  films. 
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Figure  C5.  Molecular  structures  of  some  the 
D-A  conjugaled  ^polymers  synthesized* 
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semiconductors  when  using  the  gold  source/drain  electrodes  for  charge  injection  under  ambient 
air  conditions.  Nevertheless,  the  high  carrier  mobility  the  copolymers  suggest  that  they  are 
promising  candidates  for  plastic  solar  cells  if  combined  with  a  suitable  acceptor  material. 

The  photovoltaic  properties  of  the  new  D-A  copolymers  and  those  of  poly(3- 
hexylthiophene)  (P3HT)  were  investigated  by  making  bulk  heterojunction  OPV  cells  through 
blending  with  the  [60]-fullerene  derivative  PCBM.  The  OPV  cells  were  tested  at  AM  1 .5  and  1 

sun  under  ambient  air  conditions.  The 

Table  Cl.  Solar  Cell  Parameters  of  Solvent  Annealed 
BTTP-F/PCBM  OPV  Cells  Tested  at  1  Sun. 


SVATime 

Voc  (V^ 

(mA/cm^) 

FF 

PCE(%) 

Omin 

0.743 

2,52 

0,32 

0.597 

30min 

0-748 

3.02 

0.41 

0.931 

eOmin 

0,743 

3.35 

0,44 

1,11 

QOmin 

Q-745 

_ 

Q.42 

0.984 

Figure  C7  [-V  curves  of  the  best  Bll'P- 
F;PCBM  (1:3)  cells  with  60  min  solvent  vapor 
annealing  (blue)  and  without  annealing  (red). 


P3HT:PCBM  cells  had  a  power 
conversion  efficiency  (PCE)  of  2.38  % 
((Jsc^  7.1  mA/cm^  Voc  =  0.6  V,  FF  = 
0.55).  Among  the  D-A  copolymers,  we 
found  the  BTTP-F:PCBM  (1:3)  cells  to 
be  the  most  promising  because  of  the 
superior  thin  film  quality  and  higher 
molecular  weight  of  the  copolymer.  The  thermally 
annealed  BTTP-F:PCBM  (1:3)  devices  had  a  0,84% 
PCE  (Jsc  =3.20  mA/cm^  Vo,  =  0.69  V,  FF  =  0.38). 
We  discovered  that  solvent  vapor  (chlorobenzene) 
annealing  of  the  photovoltaic  cells  significantly 
improved  the  performance  (Figure  C7).  The  solvent 
vapor-annealed  BTTP-F:PCBM  (1:3)  devices  had  a 
1 .21%  PCE  (Jso  =3.35  mA/cm^  Vo,  =  0.743  V,  FF  = 
0.44).  The  device  performance  at  different  solvent 
vapor  annealing  stages  are  shown  in  Table  Cl  and 
their  corresponding  surface  morphology  was 
observed  under  the  optical  microscope,  as  shown  in 
Figure  C8.  For  example,  the  I-V  curves  in  Figure 
C7  show  that  60  min  of  solvent  vapor  annealing 
process  has  doubled  the  power  conversion 


efficiency  of  the  solar  cells  compared  to  the  untreated 
devices,  from  0.6%  to  1.2%.  Upon  examination  of 
the  thin  film  morphology,  the  appearance  of  ~10  pm 
(diameter)  and  70  nm  (height)  regions  formed  during 
the  initial  stages  of  solvent  vapor  annealing.  At 
longer  times,  crystallites  appeared  and  after  2  hrs,  the 
entire  film  consisted  of  these  large  PCBM  crystallites. 

The  device  parameters  for  the  2  hr-anneaied  device 
are  not  shown  in  the  Table  1  because  the  crystallites 
were  so  tail  that  they  reached  the  cathode,  resulting  in 
a  much  lower  PCE  (0.4%).  We  conclude  that  the 
optimal  solvent  vapor  annealing  time  for  BTTP- 

F:PCBM  thin  films  in  OPV  cells  is  around  60  min.  In  contrast,  the  usual  thermal  annealing 
approach  reported  in  the  literature  did  not  improve  the  device  performance  in  the  BTTP- 
F:PCBM  cells. 


Figure  C8.  Optical  imagess  of  at  different 
solvent  annealing  stages:  A)  0.5hr,  B)  ihr,  C) 
l.Shr,  D)  21ir, 
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C2:  Dye-sensitized  solar  cells  (DSSC)  based  on  ZnO  (Cao) 

For  dye-sensitized  solar  cells:  We  developed  and  elaborated  the  hierarchically  structured 
photoanodes  for  much  enhanced  power  conversion  efficiency.  Although  we  started  with  zinc 
oxide  photoanodes,  we  quickly  moved  to  titania  photoanodes  and  also  studied  the  atomic  layer 
deposition  to  stabilize  the  surface  of  zinc  oxide  photoanodes.  Fig.  C9.  In  addition  to  the  wet 
chemical  processing,  we  explored  other  approaches  including  electrospray  for  the  fabrication  of 
hierarchically  structured  photoanodes. 


Figure  C9:  SEM  images  of  ZnO  nanocrystallites  with  perfect  aggregation  (top  left,  sample 
1)  to  no  aggregation  at  all  (bottom  right,  sample  4)  and  the  1-V  curves  of  ZnO 
nanocrystallites  with  different  degree  of  aggregation  showing  the  markedly  different 
power  conversion  efficiency  (right).  (Zhang,  et  ah,  2008) 


C3:  Design  of  airborne  DSSCs  (Taya) 

In  the  nature,  photosynthetic  membrane  enables  the  efficient  conversion  from  light  to  electron 
and  into  chemical  energy  (Figure  ClOa).  The  photosynthesis  made  the  fossil  fuel  over  200 
million  years  by  sunlight.  The  study  of  molecular  assembly  of  chlorophyll  derivatives  using  the 
light-harvesting  polypeptides  can  be  useful  as  a  model  of  the  natural  conversion  (Figure  Cl  Ob). 
The  mechanism  of  DSSC  is  similar  to  photosynthesis  where  photo  conversion  is  carried  out  by 
the  photo-excitation  electron  transfer  of  natural  dyes. 


Figure  CIO.  Bio-photo  harvester  (a)  photo  of  purple  bacteria,  (b)  show  the  nanostructure  of  hybrid 
structure  of  LHl,  LH2,  folding  material. 


44 


The  DSSC  comprise  four  major  parts:  rranocrystalline  titanium  oxide  (Ti02)  semiconductor  film, 
dye  sensitizer,  redox  electrolyte,  and  counter  electrode  (Figure  Cl  1).  Usually,  the  Ti02  with  Ru 
complex  dye  can  harvest  the  photon  energy  for  the  visible  light  region.  In  the  first  year,  we 
established  the  processing  route  of  DSSC  samples,  and  the  characterization  of  DSSC  data 
scientifically  (Figure  Cl  2),  and  developed  DSSC  with  higher  energy  conversion  efficiency  close 
to  8  %  (Figure  Cl 3). 


Light 


Cov-iitfT'Clatrpdr 


Figure  Cl  1 .  Schematic  model  and  Working 
mechanism  of  DSSC 


Figure  Cl  2.  Photos  of  several  new  equipment. 

However,  the  main  drawback  of  ruthenium -based  sensitizers  is  their  lack  of  absorption  in  the 
red  region  of  the  visible  spectrum.  Phthalocyanine  and  naphthalocyanine  dyes  are  well  known 
for  their  intense  absorption  in  the  red  and  near-lR  regions,  therefore  they  are  an  excellent 
candidate  as  dye  sensitizer  of  DSSC.  In  the  first  year,  we  synthesized  phthalocyanine  {Dye2)  and 
naphthalocyanine  (Dye3)  (Figure  Cl 4).  The  dyes  have  absorption  in  the  red  and  near-IR  regions 
and  adsorb  to  Ti02.  The  function  of  the  carboxylic  acid  group  of  the  dyes  is  to  graft  the  sensitizer 
onto  the  semiconductor  surface  and  to  provide  intimate  electronic  coupling  between  dyes  and 
Ti02. 
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characteristic  of  the  DSSC  sensitized 
with  the  Ru  complex.  Conversion 
efficiency  is  7.7%. 


Ru  complex  Dye2  D>e3 
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Figure  C 14.  Absorption  spectra  of  the  dyes  in  a 
solution 


Ti02  film  which  combined  Ru  and  Dye3  can  absorb  broad  light  waves.  However  the  efficiency 
of  the  DSSC  with  such  a  multi-dye  system  was  still  low  (2.6%),  because  the  energy  transfer 
between  dyes  occurred,  whose  dyes  are  aggregated.  The  aggregation  of  dyes  in  DSSC  could  be 
controlled  by  a  ligand  or  a  coadsorbent  like  a  proteins  or  a  carotenoid  of  photosynthesis. 
Learning  from  this  experience,  we  synthesized  three  unsymmeric  ZnPc  dyes,  PcS2,  PCS 5  and 
PcS6,  shown  in  Fig.  Cl 5. 


Fig.  Cl  5  Chemical  structures  of  3D  Phthalocyanine  dyes  for  use  in  the  sensitizers  of  DSSC(Mori 

et  al,  2010). 

Figure  CI6  shows  the  absorption  spectrum  of  PcS2  and  PcS6  in  solution  and  on  nano-porous 
Ti02  films  without  3a,  7a-Dmydroxy-5p-cholonic  acid  (CDCA)  as  co-adsorbant.  PcS2  on  the 
TiOj  shows  border  absorption  spectrum  around  600  nm,  suggesting  the  formation  of  aggregation 
on  the  surface.  On  the  other  hand,  PcS6  showed  sharp  Q  band  both  in  the  solution  and  on  the 
TiOj,  showing  significant  decrease  of  the  aggregation. 
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Table  1  summarizes  I-V  characteristics  of  the  DSSCs  under  one  sun  conditions  (short-circuit 
currents  open-circuit  voltages  (Tqc),  fill  factors  (FF),  and  conversion  efficiencies  (;;)),  Fig. 
C17  shows  intrinsic  power  conversion  efficiency  (IPCE)  of  those  dyes. 


Wavelength  /  nm 


Wavelength  /  nm 


Fig.  C16  Absorption  spectrum  of  two  dyes,  PcS2  and  PcS6  (Mori  et  al,  2010) 


dye 

Jsc 

(mA/cm2) 

Voc  (V) 

Fill  Factor 

PCE  (%) 

1PCE(%) 

"PcS2® 

5.3 

0.58 

0.74 

2.3 

30 

PcS5^ 

4.8 

0.58 

0.77 

2.1 

52 

PcS6' 

10.4 

0.63 

0.70 

4.6 

78 

N719‘^’‘’ 

14.4 

0.71 

0.67 

6.9 

86 

Table  C2.  1-V  characteristics  of  DSSCs^  with  various  dyes  under  one  sun  conditions  and  the 
values  of  maximum  IPCE.  Ti02  electrode  consisted  of  a  7-8  mm  thick  nanoporous  layer  with 
a  5-6  mm  scattering  layer  with  TiC14  treatment.  The  electrolyte  was  0,1  MLil,  0,6  M  DMPlml, 
0.05  M  12  and  0.5  MtBP  in  acetonitrile.The  measurement  used  an  aperture  mask  with  area  of 
0,17  cm2.  The  area  of  Ti02  was  ~  0.14  cm2. with  10  mM  CDCA  in  the  dye  solution, 
Without  CDCA,  N7I9  is  a  Ru  complex  dye. 


400  500  600  700  800 

Wavelength  /  nm 


Figure  <3 1 7.  IPCE  of  Dye-sensitized  solar 
cells  with  four  different  dyes.  The  cells 
were  the  .same  examined  onTable  C2(Mori 
etal,  2010) 
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In  contrast  to  PcS2,  PcS6  showed  the  high  performance,  77  of  4.6  %  under  one  sun  conditions. 
To  the  best  of  our  knowledge,  this  rj  value  is  the  highest  efficiency  reported  so  far  among  the 
DSSCs  using  ZnPcs  as  a  light  harvesting  dye. 

We  used  commercial  dye  N719  to  design  two  glass-based  DSSCs  and  a  flexible  PET  film  based 
DSSC  where  we  used  standard  solution  iodine  electrolyte,  see  Table  C3.  A  comparison  of  these 
designs  is  shown  in  Table  C4,  The  l-V  curves  of  these  three  DSSCs  are  shown  in  Fig.  C 1 8. 


0.05M 

O.IM 

0.6M 

0.5M 


Iodine 

Lithium  iodide 
DMPIml 


4-tert-ButyIpyridine 


in  acetonitrile  ch^cn 


0.8M  I- 


Table  C3  Electrolyte  solution  used  to  evaluate  DSSC  performance  of  Table  C4. 


Glass  DSSC 

Glass  DSSC 

Flexible  DSSC 

(2.Sninfi  FTO) 

(Umni  ATO) 

rrO-PET) 

Efficiency 

(PCE) 

10.3% 

9.6% 

6.0% 

Thickness 

5  mm 

2,2  mm 

0.3  mm 

Weight 

0,88  g  /  cm^ 

0,33g  / 

0.053  g/cm^ 

Sintering 

Temp. 

-^OO^C 

^130^(: 

Specific  PCE 

12  %(cmVg) 

29  %  (cmVg) 

10%  (cmVg) 

(PCE/wcight) 

1  : 

2.4 

;  9.4 

Table  C4  The  power  conversion  efficiency  (PCE)s  and  thickness  of  substrate,  weght ,  sintering 
condition  and  specific  PCE  of  two  glass  based  DSSCs,  one  flexible  PET  film  based  DSSC 

It  is  noted  from  Table  C4  that  the  flexible  PET  based  DSSC  has  the  highest  specific  PCE  thank 
to  its  lowest  weight  as  compared  with  glass  substrates,  even  if  the  PCE  of  the  flexible  DSSC  is 
the  lowest  among  these. 
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VnrtjiBc  /  V 


(a)  (b)  (c) 

Fig.  Cl  8  1-V  curves  of  three  different  DSSCs  corresponding  to  three  designs  in  Table  C3,  based 
on  three  different  substrates,  (a)  thick  FTO  glass,  (b)  thin  ATO  glass  and  (3)  thin  flexible  PET. 


Even  if  we  demonstated  high  PCE  of  DSSCs  based  on  both  glass  and  flexible  PET  substrate  in 
Table  C4  and  Fig.  Cl  8,  they  are  all  based  on  liquid  electrolyte  of  Table  C3,  thus,  they  may  not 
be  durable  for  longer  flight  time  of  a  MAV  or  UAV  mounted  with  such  DSSCs.  Therefore,  we 
developed  a  new  gel  electrolyte  which  is  in  form  of  gel  (or  solid  state),  thus  least  likely  for 
leakage  from  the  DSSC  package.  Table  C5  shows  the  DSSC  performance  based  on  several 
different  electrolytes,  (i)  solution  electrolyte,  (ii)  ionic  electrolyte  and  (iii)  gel  electrolyte. 


Electrolyte 


Electrolyte  for  DSSC 

operating  temperatufe  limtt 


Dye 


PCE 


Liquid  base 

mbdura  of  acelonitnie  andvalefonttrile  (85:15^ 


Tm  =  -48®C 
Tb=  81-82^ 


Tfn  =  WC 
TlJ=  139-141"C 


{-5CM:  -  80^) 


Q.ey  8MII  D83M  b.  It  1  ^  gusnldfniumltiicicv^ate. 
and  0.5M  tefT^^ylpyn«tn« 


N719 

VW*  11^- 


11.18% 


Md.  K.  riBJHrud(l^et.al.JAC5.1Z7.I6S35  CZOtS) 


Ionic  liquid  base 

1-ethyt-3^rnetnytTnidazDriurTi  tetracyafrodorate  (jEMI-TCBV 


Tm  =  less  llianO®C 

DMIl  /  EMU  / 

Tm  =  81-S3*sCTm  “  79^ 

b  /  N-  E>ut>te}«fizDjmitlau>lel  guanidinmmttiioQ^nate 
(molar  raflo:  16n2#12n.e7/3.33«J.67) 


C103 


(more  than  0«C) 


B3% 


O.  S^i  et  »i. I  Oiem.  d  112,  ITDftT  (2005) 


Gel  base 

1CNf%  PVDF41FP  »tti  34TiemcKyproptonltr1le(MPN> 


T  (ibntdHvu)^  -62®C 
Tm  =  14(M45*<; 


Tm  =  -62.9^ 
Tb=  164-165*C 


(-6<K;  -  1GCFC) 


0.6M  DMPIt.  0.1M  O.iSM  N-meth^nzimldazoie. 

1102  nano  (HJwder  (P25> 


N719 

I  -JT 

.TXTV. 


7.18% 


Z.  Huoet  Bl.,  jD(arfn«w  MsterfiiCs  fi  Sotur  Ceds.  91. 1959  (2007) 


Table  C5  DSSCs  based  on  three  different  electrolytes 


We  optimized  the  gel  electrolyte  so  as  to  increase  the  PCE  of  the  DSSC  where  the  requirement 
of  the  use  of  such  DSSC  in  lower  temperature  (up  to  -  160C)  is  met  with  this  gel  electrolyte,  the 
results  of  which  are  shown  in  Table  C6  where  the  comparison  is  between  the  DSSC  based  on 
liquid  and  gel  electrolytes. 
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Electrolyte  of  Glass  based  DSSCs 


Liquid 

EfUclefic^'  (PCE) 

nUcImess 

5mm 

Weight 

D.88  g  /  cm^ 

4.6 

SinfeHiig 

-500*C 


Gel 

7-8®/© 


i.4initi 

a'lmnt^O  ^  ITO'VirT) 
thinner  glass 

0-15>  g  /  cm^ 

I 


Specific  PCE  (PCE/weigtit) 

12  ^  41  icmV0 

1  i:  3  A 


Gel  electrolyte: 
r  1  1^.  F 


_  -  F  CPt 

rtT>F-HFP 


0.i2M  h  ,  a.2M  EBP 


0-7M  MPTml,  0-lM  TBAl 


~T 


\£ 


Ttef 


Table  C6  Two  different  electrolytes,  (i)  liquid  and  (ii)  optimized  gel  electrolyte. 
Larger  sized  DSSCs 


If  the  proposed  DSSC  is  to  be  used  for  top  surface  of  wings  in  UAV,  use  of  larger  surface  DSSC 
is  desired.  To  this  end,  we  developed  a  processing  route  for  making  larger  sized  glass  substrate 
with  Ni  mesh  pattern,  shown  in  Fig.  C 19.  This  substrate  was  processed  first  making  the  meshed 
groove  lines  by  using  Si-lithography  with  polymer  resist,  followed  by  the  electroplating  of  Mi 
into  the  groove,  then  removing  the  polymer  resist  to  complete  the  final  Mi-meshed  pattern  on 
FTO  glass  substrate.  Similar  scheme  was  applied  to  Flexible  ITO  film. 


Fig.  Cl 9  Large  sized  DSSC  of  12  cm  x  12 
cm  with  Ni  meshed  line  where  unit 
rectangular  cell  size  being  10  mm  x  5  mm. 
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Bioinspired  moth-eye  surface  for  airborne  DSSCs 

In  order  to  harvest  the  solar  energy  over  the  entire  daylight  time  when  sun  is  out,  we  employed 
Bioinspired  design  of  anti-reflection  coating,  so  called,  moth-eye  film  to  be  placed  on  the  top  of 
the  airborne  DSSC  plate. 


(a)  (b)  (c)  (d) 


(c) 


Fig.  C20  Anti-reflection  film,  (a)  synthetic  moth-eyed  film  by  blue-ray  laser  processing  ,(b)  its 
enlarged  view,  (c)  and  (d)  are  the  corneal  surfaces  of  the  Eocene  Fly  (Parker  et  al  1998),  (e)  the 
measured  profile  of  the  moth-eye  of  (a),  (b)  indicating  the  height  and  pitch  being  250  nm  and 
200  nm,  respectively. 


Fig.  C2]  (a)  Field  Experimental  set  up  to  measure  the  effectivenss  of 
using  moth-eye  anti-reflection  (AR)  film  on  glass  substrate  (b)  glass 
substrate  with  and  without  moth  anti-reflection  (AR)  film. 


51 


We  did  study  experimentally  the  transmittance  of  the  film  with  and  without  the  above  moth-eye 
film  under  realistic  sun  tight  in  Seattle.  Fig.  C2I  shows  the  experimental  set  up.  The  results  of 
transmittance  of  the  substrate  glass  with  and  without  moth-eye  film  are  shown  in  Fig.  C22  which 
indicates  the  higher  transmittance  through  the  moth-eye  film. 


Fig.  C22  Transmittance  of  solar  light 
through  glass  with  (red,  top  curve) 
and  without  (blue  bottom  curve) 
moth-eye  AR  film. 


Task  D:  Integration  of  EHSS  components  into  load-bearing  $tructure$(Hahn,Taya, 

Feraboli  and  Chasiotis  ) 

DJ  UCLA  achievements 

Multifunctional  integration  of  solar  and  thermal  energy  harvesting  structures  is  expected  to  play 
a  key  role  in  powering  a  wide  variety  of  devices  and  systems  of  interest  to  the  Air  Force, 
including  sensors  for  structural  health  monitoring  and  unmanned  aerial  vehicles  (UAVs). 
However,  little  work  has  been  reported  in  the  literature  on  the  performance  of  energy  harvesting 
structures  under  mechanical  stress. 


During  the  first  project  year,  we  have  developed  multifunctional  characterization  capability  for 
solar  energy  harvesting  devices  and  performed  a  systematic  investigation  of  the  performance  of 
commercially  available  solar  modules  under  tensile  loading.  Our  effort  resulted  in  one  MS  thesis, 
one  manuscript  submitted  for  presentation  at  a  workshop,  and  one  manuscript  in  preparation  for 
submission  to  an  archival  journal. 

1)  Structural/mechanical  characterization  of  solar 
modules 


We  have  developed  expertise  to  use  SEM,  EDX,  in 
combination  with  optical  microscope  to  identify 
device  structures  at  microscopic  levels.  Since  few 
mechanical  properties  of  solar  modules  are 
available  in  the  literature,  we  have  also  developed 
mechanical  testing  setup/procedures  to  determine 
the  Young’s  modulus  and  fracture  strength  of  solar 
energy  harvesting  modules. 
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Figure  D1  SEM  micrograph  of  a 
commercial  a-Si  solar  cell  module. 


2)  Multifunctional  characterization  of  solar  modules 

We  have  developed  a  specialized  experimental  setup  consisting  of  a  versatile  mechanical  tester 
(Instron  4483),  a  light  illuminator,  a  precision  current-voltage  source-meter,  and  a  data 
acquisition  system  for  strain  measurements.  This  setup  allovv^s  us  to  determine  functional 
(energy  harvesting)  characteristics  of  a  module  under  various  mechanical  loading  conditions. 
Tests  can  be  performed  under  monotonically  increasing  load  or  cyclic  load  conditions. 

From  the  1-V  characteristics  measured  at  different  load  conditions,  we  extract  key  metrics  of 
solar  modules,  including  the  short  circuit  current,  open  circuit  voltage,  fill  factor,  and  maximum 
power  point.  Under  the  monotonic  load  condition  shown  in  Figure  D2  (left),  the  short  circuit 
current  first  decreased  approximately  linearly  with  increasing  strain  for  commercial  a-Si  solar 
cells  we  tested.  Immediately  preceding  fracture,  the  degradation  in  the  short  circuit  current 
accelerated. 

Under  cyclic  loading  conditions,  the  short  circuit  current  density  stayed  relatively  constant  and 
equal  in  both  the  stressed  and  relaxed  state  up  until  approximately  1.4%  strain.  Above  1.4% 
strain,  the  current  density  did  not  fully  recover  when  the  load  was  removed  due  presumably  to 
plastic  deformation. 
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Figure  D2  Stress-strain  relations  and  short-circuit  current  density  of  a-Si  solar 
modules  under  monotonic  (left  column)  and  cyclic  (right  column)  loading 
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3)  Characterization  of  microstructural  changes 


We  used  optical  and  SEM  characterization  to  relate  the 
observed  performance  degradation  to  micro-structural 
changes  in  the  form  of  cracking  under  applied  stress. 
We  custom-built  a  long- working  distance  optical 
microscope  with  x-y-z  translation  capability  that  is 
compatible  with  our  multifunctional  characterization 
facility. 


Figure  D3  Custom-built  long- 
working  distant  video 
microscope. 


Localized  cracks  are  believed  to  initiate  in  the  brittle 
device  layers  of  the  solar  cell  modules  at  a  strain  level 
as  small  as  0.3%.  The  performance  does  not  degrade 
significantly  because  there  is  still  a  clear  current  path  in 
the  device.  When  the  load  is  removed,  the  cracks  can  be 
closed,  allowing  the  device  performance  to  be  reeovered. 

At  strains  larger  than  approximately  1%,  large  cracks 

develop  at  the  interface  between  two  adjacent  segments  through  the  entire  width  of  the  module, 
which  cannot  be  completely  “healed”  after  the  load  is  released  due  to  the  viscoelastic  behavior  of 
the  polymer  encapsulation  layers.  An  image  of  the  fractured  solar  module  show  that  the  polymer 
layers  did  not  fracture  and  that  the  cracks  propagate  from  the  inner  device  layers. 


Figure  D4  (b)  Typical  microscope  image  at  approximately  1.40%  in  the  stressed  state  (c)  and  the 
relaxed  state  showing  the  cracks  getting  smaller  as  stressed  is  released  from  the  cycle  test.  The 
large  crack  extends  throughout  the  entire  width  of  the  solar  module. 
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Figure  D5  (a)  Cross  section  of  an  unstressed  solar  module  where  (b)  is  zoomed  in  on  the  device 
layers,  (c)  Cross  section  after  the  cycle  test  where  (d)  shows  the  cracks  in  all  of  the  device  layers. 

The  arrows  indicate  the  load  direction. 


Accomplishments  and  Discoveries 

•  c-Si  solar  cells  are  very  brittle;  a  failure  occurs  below  0.2%-tensile  strain. 

•  A  functional  performance  of  the  a-Si  solar  cell  was  gradually  degraded  due  to  cracks  on  a 
zinc  oxide  and  amorphous  silicon  layers  above  0.38%-strain.  (Fig.  D6) 

•  After  a  uniaxial  pressure  of  0.83  MPa,  the  capacity  of  solid-state  thin-film  lithium  batteries 
started  to  drop.  For  deflections  higher  than  0.25  mm  for  the  linear  span  length  of  19  mm,  the 
charging  and  discharging  capacity  of  the  solid-state  thin-film  lithium  batteries  decreased. 
These  were  due  to  the  breakage  of  packaging  materials  and  internal  layers  such  as  lithium 
phosphorus  oxynitride  (LiPON)  solid  electrolyte,  and  the  oxidation  of  active  materials. 
Solid-state  thin-film  lithium  batteries  can  be  embedded  into  graphite/epoxy  laminates  if  the 
curing  temperature  does  not  exceed  122“C.  The  functional  performance  of  the  solid-state 
thin-film  lithium  batteries  deteriorates  rapidly  with  mechanical  fatigue  cycles.  The 
maximum  fatigue  strain  that  can  be  applied  without  degrading  the  functional  performance 
up  to  1000  cycles  is  only  0,07%.  (Fig.  D7) 

•  A  conductive  path  was  made  using  inkjet  printing  of  copper  ink.  Resistivity  of  printed 
electrode  using  our  copper  nanoparticles  was  36.7  nfJ  m.  The  printed  electrodes  remain 
fully  functional  up  to  the  maximum  strain  of  1%.  (Fig.  D8) 

•  A  power  laminate  has  been  fabricated  by  first  mounting  a  thin-film  solar  module  and  a  thin- 
film  battery  on  a  BT  core  with  the  desired  printed  conductive  pattern,  and  then  co-curing  the 
resulting  assembly  with  a  graphite/epoxy  composite  laminate.  (Figs.  D9  and  DIO) 
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•  The  energy  harvest! ng/storage  characteristics  of  the  developed  multifunctional  composite 
structure  were  measured  under  the  mechanical  static  loading  condition.  From  the  experiment, 
it  was  found  that  the  power  laminate  was  reliable  until  0.45%  of  the  static  strain.  After 
0,45%  of  the  strain  level,  the  active  layers  of  the  solid-state  th in-film  lithium  battery  was 
failed, 

•  The  Li-ion  polymer  battery  by  Dow  Kokam  was  characterized.  This  battery  was  able  to 
achieve  100%  depth  of  discharge  at  a  low  temperature  of  -ITC  under  1/12  C-rate.  Under 
compression,  the  packaging  of  the  battery  was  able  to  withstand  22.8  MPa,  and  under 
tension,  these  batteries  were  able  to  endure  1 .655%  strain  without  decrease  in  its  capacity. 
The  Li-ion  polymer  battery  should  not  be  embedded  by  co-curing  method  inside  a  CFRP 
composite  laminate  because  Us  capacity  decreased  under  fatigue  loadings.  When  the  Li-ion 
polymer  battery  was  charged,  the  thickness  of  the  Li-ion  polymer  battery  was  increased  by 
48  ±  5  pm  on  average,  compared  to  the  thickness  of  the  completely  discharged  Li-ion 
polymer  battery.  In  contrast,  when  the  Li-ion  polymer  battery  was  discharged,  the  thickness 
of  the  Li-ion  polymer  battery  decreased  by  46  ±  5  pm  on  average,  compared  to  the  thickness 
of  the  completely  charged  Li-ion  polymer  battery. 

•  An  I  PL  system  was  developed  to  sinter  different  materials  instead  of  a  thermal  chamber. 
Both  semiconductor  and  metallic  nanoparticles  can  be  sintered  on  a  composite  laminate 
using  an  intense  pulsed  light.  This  opens  up  an  avenue  for  fabricating  solar  cells  and 
conducting  circuits  directly  on  composite  laminates.  (Fig.  Dll) 

•  The  IPL  can  sinter  the  copper  nanoink.  The  sintered  copper  film  has  a  grainy  structure  with 
neck-like  junctions.  The  film  has  5  pfi-cm  of  electrical  resistivity  which  is  one  third  that  of 
previously  reported  thermally  sintered  copper  nanoink. 

•  A  silver  conductive  pattern  made  by  inkjet  printing  was  sintered  IPL.  Using  a  lumped 
capacitance  method,  it  was  found  that  IPL  irradiation  of  50  J/cm^  increases  the  temperature 
of  the  printed  silver  nanoink  pattern  by  500°C.  The  average  resistivity  of  the  fully  sintered 
silver  pattern  was  49  ±  3  nfi-m.  IPL  sintering  is  more  appropriate  for  protecting  the 
substrate  rather  than  a  conventional  thermal  sintering  process,  (Figs.  D12  and  D13) 

•  CIGS  films  were  prepared  from  CIG  metallic  alloy  and  Se  nanoparticles  using  the  IPL 
technique  without  toxic  selenization  or  vacuum  conditions.  (Fig.D14) 

•  Unidirectional  stretching  was  applied  to  improve  the  optical  and  electrical  properties  of  a 
dilute  CMT  mat.  After  the  dilute  CNT  mat  on  a  polyvinyl  chloride  film  was  stretched  by 
80%  strain,  it  reached  the  average  sheet  resistance  of  213  Q/n. 
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Fig,  D6  Crack  propagation  on  amorphous  silicon  layer  due  to  tensile  loading. 
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Fig,  D7Thin-fiim  solid-state  lithium  battery* 


Fig,  D8  (a)  Optical  microscope  image,  and  (b)  profilometer  results  of  a  printed  electrode. 
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Direction  of  Tensile 


Ekitlom  side 


Fig.  D9  Printed  and  sintered  conductive  circuit  for  power  laminate. 


(ore 

Battery 

■'  j'""" . . 

C7F  composite 


}  Ethylene  Tetrafluoroethylene  (~35  urn) 
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.T700SC/RS-30G  (Total  thickness:  530  um) 


Fig.  D-10  Cross  sectional  view  of  the  composite  laminate  with  integrated  solar  cell  and 
embedded  battery. 


Fig.  D1 1(a)  Actual  IPL  system,  and  (b)  schematic  drawing  of  IPL  system. 
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Fig.  D12  (a)  Original  BT  core,  (b)  BT  core  after  thermal  sintering  at  200_C,  and  (c)  BT 
core  after  IPL  sintering. 


Fig,  DI3  (a)  SFM  image  of  the  surface,  and  (b)  adhesive  tape  test  of  the  silver  nanoink 
pattern  after  three  consecutive  IPL  pulses  at  50  J/cm^. 


Fig.  D14  (a)  Overall  surface  morphology  SEM  image,  (b)  image  at  higher  magnification, 
(c)  cross-section  SEM  image  of  GIGS  film  prepared  on  glass  substrate  after  IPL  pulse  at 
20  J/cml 

D2  University  of  Washington  achievements 
D2‘l  Integration  of  DSSC  to  structures(Taya) 

Taya  group  at  U  W  studied  experimentally  the  residual  solar  cell  performance  of  a  DSSC  which 
was  subjected  to  mechanical  fatigue  loading  where  two  different  bonding  conditions  are  used,  (1) 
bonding  condition  1  :  each  DSSC  is  strongly  bonded  to  the  substrate  surface,  (2)  bonding 
condition  2:  top  polymer  film  covering  all  DSSCs  are  bonded  only  along  its  periphery  while  each 
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DSSC  is  merely  compressed  by  the  top  polymer  film*  In  our  experiment,  we  mounted  only  one 
DSSC  sample  on  the  substrate  which  is  then  subjected  to  bending  strain  test,  see  Fig*  D16. 


adhesive 


adhesive 


\, 


DSSCs  Bonded  birectly  to  Wing 


Bonding  condition  1 


adhesive  adhesive 


DSSCs  heid  in  place  by  Plastic  Moth-Eye  Film 
Bonding  condition  2 


Fig*  D15  Two  different  bonding  of  DSSCs  on  the  substrate,  (1)  each  DSSC  is  bonded  strongly  to 
the  substrate  surface,  (b)  polymer  adhesive  covering  all  DSSCs  is  boned  in  its  periphery  while 
each  DSSC  is  just  mechanically  compressed  to  the  substrate* 


DSSC 


Fig*  D16  Four-point  bending  test 
on  the  substrate  attached  with 
DSSC  at  its  bottom. 


The  results  of  the  cyclic  bending  tests  of  the  DSSC  mounted  on  the  substrate  are  given  in  Table 


Failure  Strain  (  %  £  ) 

Average 

Sample  # 

1 

2 

3 

4 

Bonding 
Condition  1 

Tension 

0.06% 

0.16% 

0.10% 

0.06% 

0.08% 

Compression 

0.10% 

0.06% 

0.16% 

0.24% 

0.11% 

Bonding 
Condition  2 

Tension 

0.24% 

0.24% 

0.24% 

0.24% 

0.24% 

Compression 

0.24% 

0.24% 

0.24% 

0.24% 

0.24% 

Table  Dl  Failure  strain  of  the  DSSC  mounted  on  the  substrate  subjected  to  one  time  static 
bending  test,  indicating  that  the  DSSC  with  bonding  condition  2  has  higher  failure  strain,  better 
design* 

Then  we  made  cycling  bending  test  of  the  glass-based  DSSC  with  bonding  condition  2  to  study 
the  durability  of  the  DSSC  .  The  results  are  given  in  Fig*  D17  to  find  that  the  glass  DSSC  has 
longer  life  with  even  its  PCE  slightly  increasing* 
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PCE  vs.  Bending  Cycles 

Bonding  Condition  2,  Tensile  Loading 


Fig.  D17  Residual  performance  (power  conversion  efficiency  (PCE))  of  the  DSSC  subjected  to 
cyclic  bending  up  to  40,000  cycles  with  bending  strain  of  0.1  %. 

We  also  performed  similar  cyclic  test  on  flexible  PET  DSSC  sample,  the  results  of  which  are 
summarized  in  Fig.  Di8. 


Fig.  D1 8  Normalized  PCE  of  flexible  PET  based  DSSC  subjected  to  cyclic  bending  loading  of 
up  to  5,500  cycles  with  little  degradation  of  the  residual  PCE  of  the  DSSC. 
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D2-2  Integration  of  thin  film  Li-ion  batteries  into  structural  laminated  composites:  limits  of 
process-ability  and  electromechanical  performance  (Feraboli) 

1.  INTRODUCTION 

The  flight  endurance  of  an  electric  propeller-driven  unmanned  aerial  vehicle  (UAV)  is 
significantly  improved  by  reducing  the  weight  more  so  than  increasing  the  battery  capacity 
(Thomas,  2002).  Moreover  the  structure  and  the  battery  each  contribute  typically  20-40%  to  the 
total  UAV  mass  (Thomas,  2002;  Thomas,  2004).  The  integration  of  battery  and  structure  can 
theoretically  reduce  the  total  weight  by  reducing  the  sum  of  the  structure  and  the  battery  weight. 
This  would  be  accomplished  by  exploiting  the  battery  components  as  load  bearing  elements  and 
by  eliminating  battery  fittings  or  supports.  In  addition,  rather  than  bulky,  centralized  batteries, 
the  integration  of  multiple  lightweight  batteries  into  the  structure  enables  distributed  power 
supply  and  storage,  thereby  reducing  the  amount  of  wiring.  However,  to  date,  system 
performance  improvements  achieved  with  such  integration  have  been  documented  only  for  low 
mechanical  stress  demanding  applications  (Thomas,  2005),  and  the  current  battery  technology, 
based  on  micro-structured  lithium  intercalation  compounds,  is  such  that  high  specific  energy,  as 
well  as  specific  power,  and  good  mechanical  properties  cannot  be  contemporary  achieved 
(Thomas,  2005;  Qudway,  2005;  Liu,  2008;  Snyder,  2009(a);  Snyder,  2009(b)).  In  order  to  meet 
the  requirements  for  the  next  generation  of  airborne,  load  bearing  batteries,  thin  film  Li-ion 
batteries  (TFB)  have  been  recently  proposed  (Xu,  2008;  Pereira,  2006;  Pereira,  2008(a);  Pereira, 
2008(b);  Pereira,  2009).  The  thin  film  structure  maximizes  the  specific  contact  surface  between 
electrodes  and  electrolyte,  thereby  increasing  the  electric  energy  and  power  stored  per  unit  mass 
by  increasing  the  fraction  of  reactants  and  the  rate  of  the  electrochemical  reaction  respectively. 
Moreover  scalable  manufacturing  processes  for  nanostructured  materials,  such  as  electrospinning, 
chemical  vapor  deposition  (CVD),  atomic  layer  deposition  (ALD),  are  suitable  for  TFB 
manufacturing.  The  nano-scale  technology  applied  to  the  battery  active  components  can 
potentially  lead  to  several  advantages.  It  allows  further  increase  of  TFB  capacity  and  power  by 
decreasing  the  path  length  of  Li-ion  and  electronic  transportation  (Liu,  2009;  Wang,  2008;  Scott, 
2011).  It  can  also  enable  the  simultaneous  implementation  of  electric  energy  storage  and  load 
bearing  capabilities  through  engineered  composite  electrodes  comprised  of  electrochemical ly 
active  particles  or  fibers  bound  by  a  structural  matrix  with  electrically  conductive  filler.  Lastly  it 
allows  better  accommodation  of  the  cycling  strain  induced  by  Li-ion  insertion  and  removal, 
which  causes  low  durability  and  capacity  fading  during  operation  of  conventional  bulk  material 
electrodes,  in  particular  if  Silicon  is  used  as  anodic  material  for  high  specific  capacity  (Chan, 
2008;  Trevey,  2011;  Son,  2011).  The  low  TFB  thickness,' typically  less  than  300  pm,  facilitates 
the  integration  within  thin  section  structures,  but  it  implies  that  a  significant  weight  fraction  is 
constituted  by  the  packaging  layers,  which  provide  mechanical  support  for  the  battery 
manufacturing  process,  electrical  insulation  and  sealing  of  the  highly  reactive  cell  components. 
Figure  D 19-20. 
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Figure  D19.  Perspective  photo  of  all-solid  state  thin  film  Li-ion  battery  (TFB)  with  dimensions. 
Manufactured  by  FrontEdge  Technology  Inc. 
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Anode  '  U-fttetal  (2  pm) 
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Legend: 
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Figure  D20.  TFB  Cross-sectional  schematic 


This  is  a  deterrent  for  most  large-scale  applications  of  TFBs,  which  are  indeed  employed  as 
power  sources  for  memory  chips,  microeiectromechanical  systems  (MEMS)  and  medical 
implantable  devices.  However  If  the  packaging  is  by  design  an  active  load  bearing  element, 
instead  of  a  passive  battery  mass,  then  this  technology  becomes  an  appealing  multifunctional 
system  for  electric  propelled  air  vehicles.  Low-power  applications  such  as,  for  example, 
integrated  power  source  for  an  embedded  health  monitoring  system  of  smart  composite 
structures,  can  also  be  enabled  without  the  need  of  a  load  bearing  battery  design.  The 
manufacturability  and  functionality  of  the  integration  between  structural  composite  laminates 
and  TFBs  have  to  be  proved  in  terms  of  physiochemical  compatibility  of  TFB  materials  to  the 
composite  curing  environment,  durability  of  the  electrodes,  laminate  structural  integrity  and  TFB 
packaging  integrity  under  applied  loading.  Mechanical  tests  conducted  on  a  commercially 
available  TFB  by  (Pereira,  2008)  show  that  the  battery  is  capable  to  withstand  uniaxial  transverse 


63 


pressure  up  to  830  kPa  without  any  detrimental  effects  on  the  electric  functionality,  thereby 
proving  that  the  TFB  is  compatible  with  the  autoclave  curing  pressure  of  epoxy-based  composite 
materials.  The  aforementioned  result  has  been  corroborated  by  this  laboratory.  The  same  TFB 
has  been  successfully  embedded  in  a  carbon/epoxy  composite  laminate  cured  at  12I®C,  while 
failure  has  been  documented  for  a  177®C  embedding  (Pereira,  2008),  giving  evidence  that  the 
processing  temperature  is  an  influential  parameter  for  battery  survivability.  The  TFB  ability  to 
tolerate  the  exposure  to  curing  temperature,  which  ranges  from  I2i°C  to  I77°C,  in  a  quiescent 
status  (e.g.  with  no  electric  current  flow)  is  neither  fully  covered  by  scientific  literature,  nor  by 
commercial  electronics  standards.  On  the  other  hand  the  ability  to  operate  from  -55°C  to  100°C 
has  been  extensively  investigated  because  of  the  relevant  commercial  applications.  High 
temperature  characterization  of  TFB,  having  the  same  chemistry  of  the  batteries  tested  in  the 
present  study,  has  been  performed  under  charge/discharge  cycling  (e.g.  not  in  a  quiescent  status) 
and  permanent  capacity  reduction  has  been  detected  at  80°C  (Nagasubramanian,  2004).  In  (Van 
Sluytman,  2006)  cycling  at  100“C  and  150®C  causes  permanent  capacity  reduction  associated  to 
a  decrease  in  grain  dimension  of  the  LiCo02  cathode.  When  the  battery  is  processed  at  the  same 
temperatures  in  a  quiescent  status  and  brought  back  to  room  temperature,  no  change  in  capacity 
and  no  cathode  crystallographic  change  are  noted.  However  the  state  of  charge  of  the  quiescent 
battery  is  unknown.  The  increasing  detrimental  effect  of  temperature  exposure  at  increasing  state 
of  charge  (Thomas,  2003)  up  to  thermal  runaway  (Roth,  2003)  has  been  studied  only  for 
conventional  18650  Li-ion  cells,  not  for  TFBs.  The  objective  of  this  study  is  to  experimentally 
assess  the  limits  of  process-ability  of  a  commercially  available  TFB  within  composite  curing 
temperature,  pressure  and  liquid  resin  environment,  in  order  to  enable  the  manufacturing  of 
composite  laminates  with  embedded  or  externally  bonded  TFBs.  The  research  covers  the  121  ®C 
to  199°C  temperature  range  at  two  states  of  battery  charge,  corresponding  to  the  higher  and 
lower  limit  of  the  operating  charge  level.  The  second  part  of  the  research  aims  at  investigating 
the  electromechanical  performance  of  new  structural  configurations,  such  as  external  bonding, 
and  adding  new  loading  conditions,  such  as  combined  curvature  and  strain,  to  the  existing  results 
previously  obtained  by  (Pereira,  2008;  Pereira,  2009).  Differently  from  (Pereira,  2008)  and 
(Pereira,  2009)  the  battery  is  co-bonded  or  secondary  bonded  to  the  composite  without  pre¬ 
encasing  it  in  a  protective  vinyl  ester  case,  but  only  with  a  localized  application  of  silicon 
conformal  coating  at  the  battery  leads  to  ensure  electrical  insulation  and  prevent  shorting  with 
carbon  fibers. 


2.  METHODOLOGY 

2.1  SOLID  STATE,  THIN  FILM  LI-ION  BATTERY 

All  TFBs  utilized  in  this  research  are  manufactured  by  FrontEdge  Technology  under  license 
from  the  Oak  Ridge  National  Laboratory  (ORNL),  Figure  D19.  The  same  batteries  have  been 
utilized  in  (Pereira,  2006;  Pereira,  2008(a);  Pereira,  2008(b);  Pereira,  2009).  The  cathode 
material  is  lithium  cobalt  oxide  (LiCoOa),  the  anode  is  lithium  metal  (Li)  and  the  solid  state 
ceramic  electrolyte  is  LiPON  (Li2.9PO3.3N0, 46)  (Bates,  1992).  This  battery  chemistry  has  been 
electrically  characterized  in  (Bates,  2000).  The  active  components  are  encased  by  two  muscovite 
substrates  bound  by  a  thermoplastic  layer  of  Surlyn  sealant,  leading  to  a  total  TFB  thickness  of 
150  pm,  Figure.D20.  Critical  temperature  thresholds  for  the  TFB  materials  are  summarized  in 
Table  D2. 
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Component 

Material 

Temperature 

Property 

Sealant 

Surlyn 

98°C 

Melting  point 

Anode 

Lithium 

18I°C 

Melting  point 

Electrolite 

Li2,9P03,3No.46 
(Li  PON) 

300°C 

Maximum  operating 

temperature^^ 

Cathode 

LiC02 

700°C 

Annealing  temperature^^ 

Substrate 

Muscovite 

700°C 

Calcination  temperature 

Current 

collector 

Platinum 

1768°C 

Melting  point 

1)  Yu  X.  et  al.  J.Electrochem.  Soc.  147(2),  5 1 7  (2000) 

2)  Bates  J.  B.  et  al.  J.  Eiectrochem.  Soc.  147,  59  (2000) _ _ 

Table  D2.  Relevant  temperature  thresholds  for  TFB  components 

All  the  active  components  are  grown  by  physical  vapor  deposition  (PVD)  performed  in-situ  over 
the  muscovite  substrate.  The  cathode  is  fabricated  through  RF  magnetron  sputtering  and 
annealing  at  high  temperature  (typically  700°C)  to  obtain  a  crystalline  microstructure  with  large 
grains  and  uniform  preferred  crystalline  orientation  to  maximize  ionic  conductivity  (Bates,  2000). 
The  muscovite  substrate  ensures  dimensional  stability  during  the  in-situ  cathode  annealing, 
avoiding  cathode  cracking  or  dis-bonding  due  to  thermal  stresses.  The  patented  process  is 
described  in  (Krasnov,  2011).  The  electrolyte  is  deposited  by  RF  magnetron  sputtering  in  N2 
atmosphere  (Krasnov,  2005)  according  to  the  ORNL  process  (Bates,  1992).  The  Li  anode  is 
deposited  by  thermal  evaporation  (Bates,  2000).  These  materials  are  highly  reactive.  Li-metal 
reacts  with  N2  (Bates,  1992)  and  O2  (Bates,  2000),  while  all  the  three  active  component  materials, 
react  with  H2O  (West,  2002).  Therefore  TFB  has  to  be  hermetically  sealed.  The  battery  cell  is  a 
25.4  mm  by  25.4  mm  square  with  nominal  voltage  of  4.2V,  and  a  capacity  of  I  mAh.  In  order  to 
calculate  the  actual  energy  density,  the  energy  delivered  by  one  battery  during  a  full  discharge 
has  been  measured  according  to  the  method  described  hereinafter,  and  subsequently  the  active 
components  of  the  same  battery  have  been  separated  Irom  the  substrate  and  sealant  and  weighted 
with  a  Mettler  Toledo  XS64  analytical  balance.  The  active  mass,  given  by  the  sum  of  cathode, 
electroMe,  anode  and  current  collectors,  is  0.01 13  g,  which  leads  to  an  energy  density  of  353 
Wh  kg' ,  while  the  energy  density  calculated  with  respect  to  the  total  mass  is  22  Wh  kg'*. 


2.2  MANUFACTURING  OF  THE  MULTIFUNCTINAL  COMPOSITE  LAMINATE 

Battery  survivability  is  monitored  through  discharge  and  charge  cycling  prior  to  and  post  thermal 
processing  test.  The  same  survivability  check  is  additionally  conducted  two  months  following 
the  treatment.  An  automated  circuit  board  featuring  a  charge-discharge  electronic  circuit.  Figure 
D21,  is  connected  to  a  National  Instruments  BNC-2120  connector  and  controlled  via  a  LabView 
program. 
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Analog  switches  Si  and  S2 
controlled  by  the  computer 


Figure  D21.  Schematic  of  the  circuit  used  to  test  a  battery’s  survivability  by  monitoring  charge 
and  discharge  performance  before  and  after  processing 


Figure  D22.  Description  of  the  four  discharge  and  charge  cycles  for  the  survivability  test 

Current  and  voltage  readings  are  collected  every  three  seconds.  Discharging  occurs  under  a 
constant  load  of  3.8  kfl  which  gives  a  discharge  current  of  about  1mA.  The  TFB  is  considered 
fully  discharged  when  the  voltage  reaches  3  V.  The  current  can  be  considered  constant  over  the 
entire  discharge  process,  as  shown  in  Figure  D22,  leading  to  a  discharge  rate  of  IC.  Immediately 
following  a  discharge,  charging  is  performed  at  a  constant  voltage  of  4.2  V,  applied  by  a 
Hewlett-Packard  6632A  System  DC  Power  Supply.  A  shunt  resistance  of  10  D  is  utilized  to 
measure  the  current.  The  battery  is  considered  fully  charged  when  the  current  drops  below  50  pA. 
For  each  battery,  five  charge-discharge  cycles  are  performed;  the  first  is  used  to  condition  the 
battery  and  the  following  four  provide  the  average  discharge  capacity.  Figure  D22.  The  capacity 
is  obtained  by  numerical  integration  of  the  discharge  current,  as  a  function  of  time,  over  a 
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discharge  cycle.  TFBs  are  subjected  to  pressure,  temperature  and  liquid  resin  environETient 
representative  of  composite  curing  cycles.  Batteries  are  exposed  to  a  one  hour  isothermal  hold  at 
121°C,  149°C,  I77°C  and  199°C.  At  each  temperature,  one  battery  is  tested  at  ambient  pressure, 
a  second  is  placed  under  a  26  mm  Hg  vacuum  and  a  third  is  embedded  in  a  50.8x50.8x3.8  mm 
pool  of  Fiberlay  Pro  Glas  1300  series  neat  epoxy  resin  under  a  26  mm  Hg  vacuum.  The  vacuum 
is  established  within  a  flexible  nylon  bag  in  order  to  apply  a  hydrostatic  pressure  to  the  resin. 
Figure  D23. 
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Figure  D23.  Simultaneous  thermal  treatment  of  batteries  under  ambient  pressure,  in  a  vacuum 
bag  and  embedded  in  a  neat  resin  pool,  (a)  Location  of  batteries  before  addition  of  (b)  the 
vacuum  bag 


(a)  Photo  of  thermal  processing  set-up  (b)  Circuitry  schematic 


Figure  D24.  Voltage  drop  across  the  battery  and  thermocouple  voltages  are  reeorded  by  a 
Lab  View  program  via  the  DAQ  board  during  thermal  testing  under  electrical  load,  (a)  Set-up 
before  the  specimens  tray  is  placed  into  the  oven,  (b)  Circuit  schematic  include  TFB  inside  the 
oven  and  resistive  load  (R3)  outside  the  oven 

A  series  of  batteries  are  first  thermally  processed  at  I21°C  in  a  fully  charged  state  and  in  a 
partially  charged  state  to  determine  if  the  state  of  eharge  affects  battery  survivability.  In  order  to 
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avoid  overdischarge,  as  recommended  by  FrontEdge,  the  partial  state  of  charge  is  reached  by 
discharging  to  3.9  Voc  (open  circuit  voltage)  prior  to  thermal  processing.  Based  on  the  results  of 
the  12]°C,  only  partially  charged  TFBs  are  tested  at  higher  temperature.  The  temperature  is 
recorded  by  a  LabView  program,  which  monitors  the  oven  temperature  and  the  resin  temperature 
from  two  independent  thermocouples  via  a  Mational  Instrument  USB-6210  multifunctional  data 
acquisition  (DAQ)  module.  A  Despatch  LAC  bench-top  ventilation  oven  heats  the  specimens  at  a 
rate  of  5®C  min'*  from  room  temperature  to  the  test  temperature.  Moreover  the  test  is  conducted 
under  an  electrical  load  of  108  kH  (R3)  for  better  TFB  failure  detection.  The  current  is  small 
enough  to  consider  the  state  of  charge  constant  throughout  the  test.  Voltages  and  currents  are 
recorded  by  the  DAQ  and  LabView  program  with  a  frequency  of  one  measurement  every  3 
seconds.  Figure  D24.  In  total  18  batteries  have  been  subjected  to  the  thermal  processing  test. 
Based  on  the  aforementioned  tests  results  a  cure  cycle  is  designed  and  validated  through  TFB 
embedding  tests  in  glass  fiber/epoxy  (GFRP)  and  carbon  fiber/epoxy  (CFRP)  composite 
laminates.  The  TFB  is  embedded  at  the  midplane  of  the  laminate.  Two  slits  in  the  top  layers  of 
the  laminates  allow  for  the  Flat  Flexible  Cables  (FFC),  Nicomatic  254PW01E6095  polyester 
coated  single  copper  conductor,  to  exit  the  laminate.  Figure  D25.  FCCs  are  0.25  mm  thick  and 
5.12  mm  wide,  and  they  are  rated  for  operation  up  to  1 50®C  for  short  durations  and  to  1 00®C  for 
continuous  use. 


Figure  D25.  Cross-sectional  schematic  of  a  thin-film  battery  embedded  at  the  mid-plane  of  a  four 
ply  GFRP  laminate.  Two  slits  in  the  top  two  layers  allow  for  the  FFC  to  exit  the  panel 
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Figure  D26.  Photos  of  a  battery  (a)  before  and  (b)  after  embedding  within  a  GFRP  laminate,  (c) 
Pane!  with  the  embedded  battery 

They  are  connected  to  the  TFB  electrodes  using  MG  Chemicals  silver  conductive  epoxy  8331- 
14G.  Toray  AGATE  prepreg  glass  fiber  FGF7781/2510,  8  Harness  Satin  Weave  Fabric  is  used 
for  the  embedding  test.  The  4-ply  laminate  is  cured  at  I32'’C  and  520  kPa  for  two  hours  through 
heated  press  molding.  The  laminate  stacking  sequence  is  [0,9014.  The  battery  is  located  at  the 
panel  mid-span,  Figure  D26,  A  second  embedding  test  is  performed  using  the  Toray  AGATE 
prepreg  carbon  fiber  tape  T700/25I0  with  stacking  sequence  [0/90]2s,  Figure  D27.  The 
aforementioned  fabrication  process  is  adopted  with  the  addition  of  a  localized  application  of  MG 
Chemicals  silicone  conformal  coating  422-55  at  battery  exposed  leads  to  insulate  the  connections. 
Attempts  of  embedding  the  TFB  without  insulated  connections  have  failed  due  to  electrical 
shorting  with  the  CFRP. 


Figure  D27.  Photos  of  a  battery  embedded  in  a  CFRP  laminate  showing  (a)  a  close-up  of  the 
battery  and  (b)  the  whole  panel. 
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2.3  ELECTROMECHANICAL  CHARACTERIZATION 


A  list  of  the  tests  performed  is  shown  in  Table  D3.  All  the  specimens  are  manufactured  by  press 
molding  of  IM7/977-3  prepreg  material.  Based  on  the  results  obtained  in  the  thermal  processing 
tests,  a  2  hours  curing  cycle  at  132°C  and  85  psi  has  been  adopted.  For  the  externally  bonded 
TFB  configuration  the  structural  adhesive  film  3M  AF163-2  has  been  cured  for  1.5  hours  at 
121°C  in  the  vacuum  bag.  The  double  cantilever  beam  test  (DCB)  is  conducted  according  to  the 
ASTM  D5528  standard.  Specimens  are  25.4  mm  wide,  which  is  equal  to  the  TFB  width.  Three 
specimens  have  been  tested  to  determine  the  mode  I  fracture  toughness  of  the  plane  laminate. 
Other  three  specimens  with  a  TFB  embedded  at  the  laminate  mid-plane  are  then  tested.  The 
crack  is  started  at  the  mid-plane  of  the  loaded  end  of  the  specimen,  as  in  a  standard  DCB  test, 
and  propagated  towards  the  TFB,  which  is  located  at  specimen  mid-span. 


Test _ Description _ 

Double  cantilever  beam  Measure  mode  1  interlaminar  fracture  toughness  (GIC)  of 
(DCB)  (ASTM  D5528) _ laminates  with  and  without  embedded  TFB  (3  repetitions) 

Uniaxial  tension  with  Digital  Image  Correlation  (DIC)  monitoring  of  laminate  strain 

full  field  strain  to  detect  strain  concentration  produced  by  embedded  and 
monitoring  (modified  externally  bonded  TFB  with  different  laminate  lay-ups.  (7 

ASTM  D3Q39) _ laminate  stacking  sequences  tested;  2  repetitions). _ 

TFB  full  field  thickness  Test  conducted  on  TFB  alone  to  characterize  battery  thickness 

monitoring  during  vs.  charge  level.  (2  repetitions) 

charge/discharge 

cycling _ 

Uniaxial  tension  with  Monitor  capacity  of  embedded  battery  during  and  following 

TFB  capacity  tension  test  to  identify  strain  at  electric  failure  and  failure 

monitoring _ modes.  (2  repetitions) _ 

Four  point  bend  with  Monitor  capacity  of  embedded  or  externally  bonded  TFB 

TFB  capacity  during  and  following  4PB  flexure  test  to  define  strain  and 

monitoring  curvature  at  ekectric  failure  as  well  as  failure  modes.  (4 

_ configurations;  2  repetitions) _ 

Table  D3.  Summary  of  mechanical  testing  activity 

The  goal  of  the  test  is  to  determine  the  fracture  toughness  of  the  multifunctional  laminate  and  the 
crack  propagation  path  when  crossing  the  TFB.  In  order  to  characterize  the  stress/strain  field  for 
different  laminate  configurations,  uniaxial  tension  testing  with  full  field  strain  monitoring  have 
been  performed.  A  Digital  Image  Correlation  (DIC)  technique  has  been  adopted  to  monitor  the 
surface  strain  on  the  specimen.  The  load  is  applied  up  to  failure  in  order  to  assess  the  mechanical 
strength,  whereas  the  electrical  failure  is  not  considered.  The  procedure  is  compliant  to  ASTM 
D3039,  except  for  a  specimen  width  of  76.2  mm  to  minimize  the  interaction  between  the  edge 
effect  and  the  stress  concentration  caused  by  the  TFB,  Seven  laminate  layups  have  been  tested, 
which  include  externally  bonded  battery,  interlaminate  embedding.  Figure  D28,  and  cut  out 
embedding,  which  consists  in  locating  the  battery  inside  a  cutout  in  the  mid-plane  ply. 
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Ext  e  n  so  meter 


Figure  D29.  Uniaxial  tenstion  test  setup  of  a  composite  specimen  with  embedded  TFB  and 
capacity  monitoring.  The  strain  applied  with  a  test  frame  (only  hydraufic  grips  shown)  is 
measured  with  an  extensometer. 


Two  repetitions  per  configuration  have  been  tested.  The  DIG  has  been  also  employed  for 
monitoring  the  thickness  increase  of  a  TFB  alone  during  a  single,  full  discharge  and  followed  by 
a  full  charge  cycle.  The  Li-ion  insertion  is  expected  to  cause  anode  volume  increase  during 
charge.  Therefore  the  difference  in  thickness  between  the  fully  discharged  and  the  fully  charged 
battery  status  has  to  be  assesses  carefully  in  order  to  avoid  excessive  interlaminar  normal  stress 
inside  the  laminate  and  inside  the  battery  packaging  layers.  Finally  uniaxial  tension  and  four 
point  bending  tests  with  TFB  capacity  monitoring  are  conducted  to  assess  the  electric  integrity. 
Two  tension  test  specimens  are  tested  with  an  interlaminar  embedding  in  order  to  apply  a 
uniform  strain  to  the  battery.  The  specimens  are  304.8  mm  long  and  76.2  mm  wide.  The  battery 
is  connected  to  the  charge/discharge  circuit  by  means  of  FFCs,  Figure  D29.  The  four  point  bend 
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tests  are  conducted  for  applying  curvature  or  combined  curvature  and  in-plane  strain.  This  is 
achieved  by  interlaminar  embedding  and  external  secondary  bonding  respectively:  while  the  first 
leads  to  pure  curvature,  the  second  leads  to  curvature  plus  compressive  or  tensile  strain 
depending  on  which  side  of  the  specimen  the  TFB  is  bonded  on.  For  the  four  point  bending  test  a 
plane  laminate  is  also  tested.  The  mechanical  properties  of  the  TFB  packaging  materials  and  the 
CFRP  are  listed  in  Table  D4, 


Material 

Density 

Modulus  of 

elasticity 

In-plane 

modulus 

elasticity 

^^^of  '  fracture 

toughness 

Strain  to  failure 

Surlyn 

0.96  g/cm^ 

282  Mpa 

106  MPa 

1250  JW 

80000  pstrain 

Muscovite 

2.60  g/cm^ 

176  GPa 

64.3  GPa 

891  mJ/m^ 

unknown 

90°  -  7700 

iM7/977- 

3 

1.63  g/cm^ 

90°  -  8.34  GPa 
0°-  162  GPa 

4.96  GPa 

316  J/m^ 

pstrain 

0°  -  14600 

pstrain 

1)  Compston  P,  et  al.,  J  Mater  Sci  Lett,  21,  383-386  (2002) 

2)  Giese  Jun  R.F.,  Nature,  248,  580-581  (1974) 

Table  D4,  Mechanical  properties  of  TFB  packaging  components  and  prepreg  com  posite  material 


3*  RESULTS 

3,1  LIMITS  OF  PROCESSABILITY 

Typical  temperature  and  voltage  profiles  for  the  thermal  processing  test  are  shown  in  Figure  D30, 
During  ramp-up  the  resin  temperature  shows  a  spike  due  to  the  exothermic  reaction  associated  to 
the  resin  cross  I  inking.  However  the  maximum  difference  between  the  resin  and  the  oven 
temperature  during  ramp-up  never  exceeds  15°C,  Moreover  the  temperature  peak  at  crosslinking 
always  occurs  during  temperature  ramp-up  and  it  never  exceeds  the  isothermal  hold  temperature, 
which  is  therefore  the  maximum  temperature  experienced  by  the  TFB,  The  battery  voltage 
decreases  with  increasing  temperature  and  it  recovers  its  initial  level  when  the  batteries  are 
removed  from  the  oven  at  the  end  of  the  one  hour  long  isothermal  phase.  The  voltage  of  the 
battery  embedded  in  neat  resin,  however,  shows  a  negative  spike  which  corresponds  with  the 
beginning  of  the  exothermic  event.  The  spike  exists  for  all  TFBs  processed  in  the  neat  resin, 
lasting  3  to  9  seconds,  without  having  a  detrimental  effect  on  survivability.  While  the  reason  for 
the  voltage  drop  is  unknown,  electric  conductivity  in  epoxy  resins,  which  could  cause  electronic 
flow  between  the  battery  exposed  leads,  is  known  to  be  caused  by  impurities,  with  Na^  and  Cl" 
ions  normally  being  the  charge  carriers  (Gal lone,  2001), 
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Figure  D30.  (a)Typical  temperature  cycle  for  TFB  thermal  processing  test.  The  temperature 
profile  for  the  oven  and  neat  resin  is  shown  for  a  1  hour  isotherm  at  121®C.  (b)  Voltages  of  the 
three  batteries  shown  in  Figure  D23  undertaking  the  thermal  treatment  recorded  in  (a).  A  voltage 
drop  of  the  battery  embedded  within  the  resin  coincides  with  heat  generation  by  the  resin  due  to 
exothermal  crosslinking  reaction 
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Figure  D3L  A  discharge  profile  from  the  survivability  test  for  a  fully  charged  battery 
subjected  to  a  1  hour  isotherm  at  121°C  under  ambient  pressure.  The  decrease  in  the  discharging 
time  is  proportional  to  the  6%  decrease  in  capacity  retention  after  processing  (see  Table  D3) 

The  voltage  drop  does  not  occur  with  the  embedding  in  the  Toray  2510  resin  system.  Thermal 
testing  at  121'^C  of  fully  charged  batteries  always  leads  to  battery  capacity  reduction,  as 
summarized  in  Table  D5,  leading  to  an  average  capacity  reduction  of  9%,  A  reduction  in 
discharge  time,  proportional  to  the  capacity  reduction,  is  noted,  while  the  cell  voltage  and  power 
are  unchanged,  Figure  D3L 


Fully  Charged _  Partially  Charged 

Capacity  Survival  Capacity  Survival 

Retention  Rate  Retention  Rale 


Ambient 

94% 

1/1 

102% 

1/1 

Vacuum  Bag 

95% 

1/1 

102% 

1/1 

Neat  Resin 

84% 

1/1 

102% 

1/1 

Table  D5.  Capacity  retention  of  batteries  undergoing  a  1  hour  isothermal  hold  at  121^C  in  either 
a  fully  charged  state  or  partially  charged  to  3.9  Vqc  (open  circuit  voltage) 

When  TFBs  are  partially  discharged  to  3.9  V  before  being  tested,  they  withstand  thermal 
processing  up  to  149°C  without  any  detrimental  effect  on  their  electric  performance.  The  detailed 
summary  of  the  results,  shown  in  Table  D6  and  D7,  shows  that  the  short  term  and  long  term 
capacity  after  thermal  processing  oscillates  within  ±3%  with  respect  to  the  baseline  capacity 
measured  before  the  test  The  capacity  variation  is  associated  to  a  slight  offset  of  the  discharge 
current  and  voltage  curves,  Figure  D32.  Either  the  liquid  resin  or  the  vacuum  bag  environment 
does  not  seem  to  have  an  effect  at  this  processing  temperature. 
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Figure  D32.  A  discharge  profile  from  the  survivability  test  for  a  partially  charged  battery 
subjected  to  a  1  hour  isotherm  at  121°C  under  ambient  pressure.  The  battery  had  102%  capacity 
retention  after  processing  (see  Table  D3) 
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Figure  D33.  A  discharge  profile  for  a  partially  charged  battery  subjected  to  a  1  hour  isotherm  at 
177°C  under  vacuum.  The  battery  had  91%  capacity  retention  after  processing,  ultimately  failing 
with  aging  (see  Table  D4) 
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Isothermal  Temperature 

i2rc 

149°C 

\77°C 

199°C 

Ambient 

Vacuum 

i/i 

1/1 

1/2 

0/1 

bag 

1/1 

1/1 

(2)/2 

0/1 

Neat  Resin 

1/1 

1/1 

0/2 

0/1 

Table  D6.  Survival  rate  of  partially  charged  batteries  subjected  to  temperatures,  pressures  and 
liquid  resin  environment  typical  of  composite  manufacturing*  Parentheses  indicate  failure  after 
two  months  of  aging 


Isothermal  Temperature 

]21°C 

149“C 

177°C 

199X 

Ambient 

102% 

(100%) 

100% 

(100%) 

76%  (70%) 

0% 

Vacuum 

bag 

102% 

(100%) 

103% 

(100%) 

87%  (0%) 

0% 

Neat  Resin 

102% 

(100%) 

98%  (98%) 

0% 

0% 

Table  D7,  Average  discharge  capacity  retention  of  survived  batteries  from  Table  D4,  Values  in 
parenthesis  measured  two  months  after  test 
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Figure  D34.  A  discharge  profile  from  the  survivability  test  for  a  partially  charged  battery 
subjected  to  a  1  hour  isotherm  at  177°C  under  ambient  pressure*  The  battery  had  76%  capacity 
retention  after  processing,  reducing  to  70%  with  two  months  of  aging 
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All  the  batteries  tested  at  177°C,  either  at  ambient  pressure  or  in  the  vacuum  bag,  are  affected  by 
partial  or  total  capacity  loss.  Table  D6  and  D7,  associated  to  a  decreased  discharge  time  and  a 
negative  translation  of  the  current  and  voltage  discharge  profile.  Figure  D33  and  D34.  The 
batteries  embedded  in  neat  resin  at  177°C,  as  well  as  all  the  ones  processed  at  199°C  failed 
during  testing.  Three  distinct  types  of  failures  have  been  recognized  and  characterized  based  on 
optical  microscopy  at  anode  side  and  capacity  retention,  Figure  D35(a-d).  Type  I  failure.  Figure 
D35(b,bb),  is  observed  as  a  localized  grey  spot  on  the  Li-anode  with  complete  loss  of  grain 
boundaries.  FrontEdge  attributes  these  observations  as  electronic  failures  consisting  of  local 
breakdown  of  the  electrolyte.  The  grey  spot  is  always  associated  to  bubbling  of  the  overlaying 
Surlyn  sealant  layer  at  failure  location.  This  failure  occurs  in  batteries  processed  at  ambient 
pressure  or  within  the  vacuum  bag.  The  affected  batteries  are  operational  but  with  a  reduced 
capacity.  Type  II  failure  produces  a  neutral  grey  discoloration,  which  contacts  at  least  one  edge 
of  the  active  component,  Figure  D35(c,cc),  The  micrographs  reveal  dark  patches  intermixing 
with  the  pristine  anode.  It  is  believed  that  this  failure  is  caused  by  the  reaction  of  the  Li-anode 
with  contaminants  diffusing  through  the  Surlyn  sealant  and  entering  into  the  battery.  The 
appearance  of  Type  II  failure  is  always  associated  to  total  battery  failure.  This  failure  is  observed 
after  thermal  processing  of  batteries  embedded  in  neat  resin,  or  in  aged  batteries  previously 
processed  at  177°C.  The  Type  III  failure  shown  in  Figure  D35(d,dd)  has  occurred  for  all  the 
batteries  tested  at  1 99°C. 


Pristine  anode  Type  I  taflure 

At-r«c*iv*d,  121°C.  149*0  ITTC  (it  unbtant  & 

in  tha  V9Giium 


Type  II  fallufe  Type  III  failure 

177'^C  (in  n«n  resln|  1 99X 


Figure  D35.  Anode  side  photos  of  the  three  observed  types  of  failures  ((b),  (c)  and  (d))  compared 
to  an  as-received  battery  (a).  Corresponding  dark-field  optical  micrographs  are  shown  in  (aa) 
through  (dd).  Type  I  failure  is  associated  with  a  local  breakdown  of  the  electrolyte  and  bubbling 
of  the  sealant  (Surlyn  M.P.  of98°C).  Reaction  of  the  anode  with  contaminants  results  in  a  Type 
11  failure.  Type  111  failure  occurs  above  the  Li  melting  temperature  (18]°C).  Type  I  leads  to 
partial  capacity  loss;  Type  II  and  III  lead  to  total  capacity  loss. 
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(a)  As-received  (b)  After  testing  (c)  Aged  after  testing 


Figure  D36.  A  battery  thermally  tested  at  177°C  in  the  vacuum  bag  for  an  hour  showing  (b)  the 
formation  of  Type  I  failures  immediately  after  testing,  followed  by  (c)  Type  il  formations  with 
two  months  of  aging.  Immediately  after  testing,  the  cell  exhibited  87%  capacity  retention, 
however,  completely  failed  with  aging 


Figure  D37.  The  appearance  of  bubbles  and  a  disbonded  front  (b)  along  the  edge  the 
active  components  and  (c)  over  the  leads  after  121°C  processing.  In  some  as-received  batteries, 
bubbles  are  apparent  along  the  edge  of  the  active  components,  (d)  shows  the  Surlyn  which  has 
been  squeezed  out  of  between  the  muscovite  layers  (exhibiting  a  slight  offset);  this  is  typical  of 
the  vacuum  bagged  specimens. 

The  failure  occurs  above  177°C  during  temperature  ramp-up,  leading  to  sudden  loss  of  voltage 
and  battery  failure.  The  failure  is  due  to  lithium  melting  (M.P,  18iX),  The  TFB  becomes  black 
at  the  anode  side,  occasionally  showing  gray  spots  which  are  probably  Type  1  failures  formed 
before  the  anode  melting.  The  battery  shown  in  Figure  D36  is  affected  by  Type  I  failure  and 
associated  capacity  reduction  immediately  after  processing  at  177°C  in  the  vacuum  bag,  but 
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completely  fails  by  Type  II  failure  with  aging.  This  indicates  that  the  failure  Type  II  is  time- 
dependent.  Additionally,  a  similar  result  is  observed  when  the  battery  seal  is  artificially 
compromised  in  the  ambient.  Hence  the  low  temperature  thermoplastic  sealant  (M.P.  98°C) 
seems  to  be  the  limiting  factor  for  survivability  from  a  processing  standpoint.  Further 
observation  of  treated  batteries  reveals  sealant  bubbling  and  flowing  even  for  TFBs  tested  at 
12PC,  leading  to  the  formation  of  a  disbonded  front  along  the  edges  of  the  active  components 
(Figure  D37(b))  and  over  the  leads  (Figure  D36(c)).  For  the  TFBs  processed  under  vacuum  bag 
the  Surlyn  is  seen  to  have  been  squeezed  out  of  the  battery  edges  (Figure  D37(d)).  However  the 
functionality  of  these  batteries,  as  tested  for  survivability  after  testing  and  after  two  months  aging, 
is  not  affected.  Inspection  of  as-received  batteries  revealed  that  some  batteries  have  long  strands 
of  bubbles  along  the  edges  of  the  active  components.  Heating  of  the  battery  leads  to  melting  of 
the  Surlyn  allowing  for  mass  and  void  diffusion.  However,  no  correlation  has  been  found 
between  the  appearance  or  location  of  bubbles  and  battery  survivability.  Finally  successful 
embedding  tests  with  full  capacity  retention  are  performed  by  co-curing  the  TFB  within  a 
composite  laminate  through  press  molding  of  prepreg  materials  at  132°C  and  517  kPa,  with  a 
cure  time  of  2  hours,  Figure  D27-29,  Table  D8. 


Capacity 

Survival 

Retention 

Rate 

CFRP 

103% 

1/1 

GFRP 

102% 

1/1 

Table  D8.  Survival  of  partially  discharged  batteries  embedded  in  either  a  GFRP  or  CFRP 
laminate  cured  at  132^^0  and  520  kPa 


3.1  ELECTROMECHANICAL  PERFORMANCE 

When  the  delamination  that  propagates  trough  the  CFRP  laminate  reaches  the  battery  it 
propagates  with  an  unstable  fashion  throughout  the  battery,  Figure  D38.  Further  data  reduction 
of  the  DCB  test  show  that  the  average  fracture  toughness  of  the  laminate,  calculated  over  the  net 
crack  surface,  defined  as  the  CFRP  surface  only,  is  equal  to  the  average  fracture  toughness  of  the 
plane  laminate.  Therefore  the  TFB  behaves  as  a  perfect  defect,  similarly  to  a  Teflon  insert.  A 
expected  from  the  very  low  fracture  toughness  of  the  muscovite.  Table  D4,  the  crack  propagates 
through  the  muscovite  substrate,  Figure  D39.  The  uniaxial  test  under  DIG  monitoring  shows  that 
the  interlaminar  embedding  does  not  lead  to  surface  strain  concentrations.  Figure  D40,  Moreover 
the  strain  distribution  is  reasonably  well  predicted  by  the  Classical  Laminate  Theory  (CLT), 
Figure  D40  and  Table  D9,  The  good  agreement  between  the  experimental  and  calculated  strain 
seems  to  confirm  that  the  TFB  shares  load  with  the  laminate  as  desired. 
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Figure  D38.  Typical  DCB  load  displacement  curves  for  a  laminate  with  embedded  TFB  at  the 
mid-plane 


Figure  D39,  Crack  surfaces  of  DCB  specimen  with  embedded  TFB  showing  the  muscovite 
substrate  left  on  both  sides* 
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Figure  D40.  Uniaxial  tension  test  of  laminate  with  interlaminar  embedding  at  laminate  mid¬ 
plane;  Laminate  stacking  sequence  is  [0/90/TFB/90/0].  (a)  Surface  strain  field  measured  with 
Digital  Image  Correlation  (DIC).  (b)  Comparison  between  measured  strain  at  section  a-a  and 
calculated  strain  with  Classical  Laminate  Theory  (CLT) 
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Figure  D41.  Uniaxial  tension  test  of  laminate  with  externally  bonded  TFB.  Laminate  stacking 
sequence  is  [TFB/(0/45/90/-45)S],  (a)  Surface  strain  field  measured  with  Digital  Image 
Correlation  (DIC).  (b)  Comparison  between  measured  strain  at  section  a-a  and  calculated  strain 
with  Classical  Laminate  Theory  (CLT)  and  shear  lag  theory 


^tpb/^o 

Lay-up 

Experimental 

(DIC) 

CLT 

[0/90ArFB/90/01 

0.75 

0.65 

[02A’FB/02] 

0.89 

0.79 

[0/45/90/-45ArFB/- 

45/90/45/0] 

0.80 

0.74 

^TFB:  average  strain  over  TFB 
average  far  field  strain 

Table  D9  .  Uniaxial  tension  testing  with  interlaminar  embedding.  Comparison  between  measured 
and  calculated  strain  with  classical  laminate  theory  (CLT) 
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Figure  D42*  Micrographs  showing  delamination  failure  of  the  TFB  substrate  under  uniaxial 
tension.  (a)Ful!  specimen  section  with  TFB  lifted  apart  from  the  laminate  throughout  the  entire 
battery  length.  (b)Multiple  delaminations  propagated  through  the  lower  substrate,  (c)  A  thin 
muscovite  layer  left  over  the  adhesive  film  throughout  the  whole  crack  surface. 


(a)  (b) 

Figure  D43.  (a)  DIC  contour  plot  of  z-displacement  at  battery  full  charge  level,  (b)  Average  z- 
displacement  at  the  anode  surface,  measured  at  increasing  charge  levels. 


The  slightly  higher  values  of  measured  strain,  with  respect  to  the  calculated  ones,  might  be  due 
to  the  higher  stiffness  of  the  center  section  of  the  specimen  compared  to  the  side  section.  This 
local  increase  of  stiffness  caused  by  the  TFB,  concentrates  the  stress  towards  the  center  of  the 
specimen  and  decreases  the  average  stress  along  the  edges.  Another  reason  could  be  the 
uncertainty  in  the  muscovite  modulus  of  elasticity  utilized  as  input  for  the  CLT  calculation.  The 
strain  distribution  relevant  to  the  cutout  embedding  is  affected  by  strain  concentrations.  The 
strain  increase  ranges  from  the  16%  of  the  unidirectional  laminate  to  the  38%  of  the  crossply 
laminate.  The  concentrations  are  along  the  transverse  TFB  edges.  The  reason  of  the  strain 
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concentrations  seems  to  be  the  stiffness  mismatch  between  the  TFB  and  the  cutout  ply.  Hence 
the  cutout  embedding  is  not  recommended.  The  externally  bonded  TFB  shows  a  non-uniform 
strain  distribution  because  the  TFB  is  in  a  shear  lug  condition.  In  fact  the  battery  length  is  not 
enough  to  allow  for  the  stress/strain  field  to  stabilize  and  recover  a  plane  stress  condition.  Figure 
D41 .  This  is  also  due  to  the  low  shear  modulus  of  elasticity  of  the  Surlyn  sealant.  The  two  strain 
spikes  located  at  the  battery  edges  are  apparent  strain  caused  by  the  relative  displacement 
between  the  upper  and  lower  TFB  substrate.  From  the  strength  point  of  view  the  TFB  does  not 
hace  detrimental  effects  on  the  laminate  strength,  but  mechanical  failure  occurs  at  the  battery 
substrate  when  the  applied  strain  reaches  4898  pstrain.  Table  DIO,  The  failure  mode  is 
muscovite  delamination,  which  is  due  to  shearing  at  the  battery  edges  and  low  fracture  toughness. 
The  micrograph  of  the  crack  surfaces  demonstrates  that  the  muscovite  adherent  failure  occurs 
across  the  entire  crack  surface,  while  the  adhesion  at  the  adhesive/muscovite  interface  remains 
intact.  Figure  D42.  The  TFB  thickness  increase  during  a  full  charge  has  been  measured  using  the 
Die.  As  shown  in  Figure  D43,  although  the  measurement  is  affected  by  a  fairly  high  noise  level, 
the  trend  has  been  clearly  captured-  The  total  increase  is  about  4  pm,  which  is  equal  to  the  40% 
of  the  total  thickness  of  the  active  components.  Therefore  this  expansion  has  to  be  considered  as 
a  design  load  case  for  the  next  generation  of  structural  batteries.  The  mechanical  tests  conducted 
with  TFB  capacity  monitoring  demonstrate  that  the  applied  strain  and  curvature  does  not  affect 
battery  operation. 


strain  [^strain] 

Figure  D44.  TFB  capacity  measured  under  increasing  applied  strain  cycles  up  to  electric  failure. 
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Figure  D45.  Four  point  bend  test  with  standard  60:1  span/thickness  ratio  at  incipient  laminate 
failure.  Stacking  sequence  [(0/45/90/-45)3/TFB]s 
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Figure  D46.  TFB  capacity  measured  under  cycles  of  increasing  curvature  up  to  electric  failure 


(a)  (b) 

Figure  D47.  (a)  Four  point  bend  test  with  standard  60:1  span/thickness  ratio  at  TFB  failure. 
Stacking  sequence  [TFB/(0/45/90/-45)3S].  The  battery  is  externally  bonded  onto  the  specimen 
surface  that  faces  compression,  (b)  Failure  occurs  by  disbonding  buckling  of  the  upper  TFB 
substrate  and  subsequent  battery  contamination. 
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Strain _ _ 

Rep.  1  5072  pstrain 

Rep.  2  4724  pstrain 

Avg. _ 4898  pstrain 

Table  DIO  .  Uniaxial  tension  test  with  externally  bonded  TFB;  strain  values  at  TFB  mechanical 
failure 


Strain 

Rep.  1 

6000  pstrain 

Rep.  2 

6750  pstrain 

Avg. 

6375  pstrain 

Table  Dll.  Uniaxial  tension  test  with  TFB  capacity  monitoring;  strain  values  at  electric  failure 


Rep,  1 
Rep.  2 
Avg. 

Radius 

curvature 

of  . 

strain 

124  mm 

99  mm 

1 12  mm 

0  [istrain 

0  jistrain 

0  jistrain 

Rep.  1 

504  mm 

-2091  pstrain 

Rep.  2 

557  mm 

-1890  pstrain 

Avg. 

531  mm 

-1990  pstrain 

Rep.  1 

204  mm 

5160  pstrain 

Rep.  2 

224  mm 

4688  pstrain 

Avg. 

214  mm 

4924  pstrain 

Table  D12  .  Four  point  bend  test  with  TFB  capacity  monitoring;  TFB  curvature  and  strain 
values  at  electric  failure 


When  the  uniaxial  strain  is  applied  the  battery  capacity  remains  constant.  Figure  D44,  up  to 
failure,  which  occurs  on  average  at  6375  pstrain.  The  causes  of  the  failure  are  unknown  since  the 
battery,  being  embedded  in  the  specimen,  is  not  accessible.  However  the  strain  at  electric  failure 
is  only  slightly  higher  than  the  strain  at  mechanical  failure  of  the  muscovite  reported  above, 
which  can  be  an  indication  that  in  this  case  a  mechanical  failure  of  the  packaging  has  caused 
TFB  contamination.  Curvature  without  strain  is  applied  through  a  four  point  bending  of  a 
specimen  with  the  TFB  embedded  at  the  laminate  mid-plane.  This  is  the  only  loading  case  where 
the  laminate  mechanical  failure,  which  occurs  at  a  radius  of  curvature  of  159  mm.  Figure  D45,  is 
more  critical  than  the  electrical  failure  which  occurs  at  1 12  mm.  Table  D12.  The  laminate  critical 
radius  of  curvature  is  relevant  to  a  quasi-isotropic  laminate  tested  at  a  span  to  thickness  ration  of 
60:1.  In  order  to  achieve  the  critical  radius  of  curvature  for  the  TFB  the  ratio  has  been  increased 
to  a  non-standard  145:1.  In  the  other  cases  of  combined  curvature  and  in-plane  strain,  electrical 
failure  has  been  obtained  with  a  60:1  ratio  and  the  results  are  shown  in  Table  D12.  As  is  can  be 
see  by  the  very  low  absolute  value  of  the  critical  strain,  the  compressive  case  is  particularly 
critical  because  it  leads  to  early  failure  by  TFB  substrate  disbonding  buckling  and  subsequent 
battery  contamination  failure.  Figure  D47.  This  is  due  to  the  low  modulus  and  large  thickness  of 
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the  sealant,  and  furthermore  it  is  possibly  related  to  the  sealant  bubbling  and  disbonding  that 
occurs  during  curing,  as  documented  in  Figure  D37. 

4.  CONCLUSIONS 

This  study  identifled  the  limits  of  processing-ability  of  solid  state  thin  film  lithium  batteries 
embedded  into  composite  laminates.  Cure  temperature  is  the  most  influential  parameter  for 
battery  survivability  during  composites  manufacturing.  Successful  embedding  tests,  with  full 
capacity  retention,  have  been  performed  with  glass  fiber/epoxy  and  carbon  fiber/epoxy  cured  at 
I32®C  and  517  kPa.  With  proper  procedures  it  is  possible  to  successfully  cure  the  battery  inside 
the  laminate  up  to  149°C.  If  higher  processing  temperatures  are  reached,  either  locally,  due  to 
resin  exothermic  reaction,  or  by  design,  the  battery’s  electrical  performance  will  deteriorate. 
Failures  are  associated  either  to  the  TFB  polymeric  sealant  failure,  or  to  physical-chemical 
degradation  of  the  electrolyte  or  the  Li-anode.  Ongoing  research  is  focused  on  the  mechanical 
and  electrical  characterization  of  laminates  with  embedded  batteries  under  applied  strain  and 
curvature.  The  substrate  (muscovite)  brittleness  and  low  fracture  toughness  is  most  important 
limiting  factor  for  the  battery  integrity  under  mechanical  loading.  The  sealant  (surlyn)  low- 
modulus  is  second-most  important  limiting  factor.  Embedded  TFB  do  not  have  a  detrimental 
effect  on  mechanical  properties,  they  show  good  adhesion  and  load  sharing  with  the  composite 
structure.  The  adhesion  to  the  laminate  is  stronger  than  the  fracture  toughness  of  the  muscovite 
substrate.  The  lectrical  failure  occurs  at  about  50%  of  the  strain  at  failure  of  structural  composite. 
In  addition  the  embedding  method  and  the  relative  stiffness  mismatch  between  battery  (both 
sealant  and  substrate)  and  CFRP  laminate  are  critical  parameters  for  strain  concentration.  A  more 
ductile  substrate  and  a  sealant  with  higher  modulus  would  increase  load  share  and  limits  of 
utilization  of  externally-bonded  or  embedded  batteries.  TFB  shape  can  also  be  optimized  to 
minimize  strain  concentration.  The  successful  embedding  of  all-solid-state  thin-fllm  batteries 
into  composite  laminates  can  find  applications  in  commercial,  military  and  aerospace  industries. 
Based  on  the  findings  of  this  study  a  prototype  of  energy  self-sustained  structure,  which 
incorporates  flexible  solar  arrays  and  TFBs  to  harvest  and  store  energy,  has  been  manufactured. 
Figure  D48. 


Figure  D48.  Working  prototype  of  energy  self-sufficient  CFRP  stiffened  panel  with  integrated 
TFBs  and  flexible  solar  array 

D3  Study  on  the  residual  properties  of  thin  film  EHSS  (Chasiotis) 

Introduction 
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In  this  research  program  we:  (a)  investigated  the  residual  stresses  of  Transparent  Conductive 
Oxide  (TCO)  and  the  p-n  junction  diode  thin  film  layers  in  commercial  thin  film  solar  cells  by 
taking  advantage  of  straight  and  “zig-zag”  blister  formation  and  by  measurements  of  the 
curvature  TCO-Si  and  Si  thin  films  with  1  pm  thickness;  (b)  conducted  a  series  of 
electrochemical  and  mechanical  tests  with  thin-film  composite  Si  (Si  +  acetylene  black  +  PVDF) 
anodes  and  lithium  cobalt  oxide  (LiCoOj  +  acetylene  black  +  PVDF)  cathodes.  In  the  first  task, 
we  observed  that  when  both  the  TCO  and  p-n  junction  were  present,  telephone  cord 
delaminations  were  formed,  while  when  the  TCO  was  removed,  only  straight  delaminations 
formed.  Calculations  of  the  mean  residual  stresses  from  straight  delaminations  of  Si  films 
resulted  in  -463±1 1 5  MPa,  while  telephone  cord  delaminations  of  TCO-Si  layers  provided  a 
residual  stress  of  -742±104  MPa.  Using  strips  of  the  TCO/p-n  junction  and  the  p-n  junction  alone 
the  residual  stress  gradient  of  the  bilayer  was  found  to  be  170±1 5  MPa/pm,  while  for  the  Si 
monolayer  was  172±12  MPa/pm.  In  the  second  task  we  evaluated  the  interlaminar  shear  strength 
between  the  composite  Si  layer  and  the  Cu  substrate  in  an  anode  as  123+48  MPa.  Si  anodes  were 
lithiated  and  delithiated  under  galvanostatic  and  potentiostatic  conditions,  and  the  morphological 
surface  changes  were  optically  recorded.  The  surface  morphology  changed  mainly  during 
lithiation  but  there  were  no  significant  changes  observed  in  the  surface  images  for  the  cathode 
samples  after  lithiation  and  delithiation.  In-situ  observations  showed  a  thick  solid-electrolyte 
interphase  (SEI)  forming  on  the  electrode  surface.  Finally,  the  LiCo02  composite  layer  was 
found  to  be  structurally  weak  under  small  tensile  stresses. 
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RESIDUAL  STRESSES  AND  STRENGTH  OF  THIN  FILM  SI  SOLAR  CELLS 


Thin  film  photovoltaic  cells  require  a  fraction  of 
the  material  used  in  conventional  solar  cells 
while  maintaining  relatively  high  efficiency 
(Fthenakis,  2009).  Thin  film  photovoltaic  cells 
are  comprised  of  a  semiconductor  pm  junction 
diode  which  is  responsible  for  converting  the 
energy  of  photons  to  electricity,  a  transparent 
conductive  oxide  (TCO)  which  acts  as  a 
protective  layer,  and  a  bottom  contact.  The 
circuit  is  completed  with  a  metallic  grid  on  the 
surface  of  the  cell  acting  as  the  top  contact.  The 
most  common  semiconductor  p-n  junction 
material  is  silicon  (Si),  which  is  investigated  In 
this  project,  while  the  TCO  that  provides  a  top 
metallic  transparent  contact,  enables  the 
transport  of  charge  carriers,  and  acts  as  a 
protective  layer  for  the  p-n  junction  was  zinc 
oxide  (ZnO).  A  Kapton  bottom  contact  and  a 
thick  A I  substrate,  with  approximate  thicknesses 
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Figure  D66.  Cross-section  of  solar  cell  thin 
film  showing  the  TCO  at  the  top  (light 
gray),  the  p-n  junction  diode  in  the  middle 
(dark  grey),  and  contact  at  the  bottom  (thin 
light  gray  layer). 


of  100  nm  and  50  pm,  respectively,  completed  the  thin  film  solar  cells  which  were  obtained  from 
commercial  sources.  Figure  D66(a)  shows  a  cross-section  of  the  solar  cell  film.  Quantification  of 
the  mechanical  properties  of  the  thin  film  materials  in  a  solar  cell  is  a  cumbersome  task  due  to 
the  small  film  thicknesses  and  the  formation  of  cracks  when  bending  or  extensional  loads  were 
applied.  Figures  D67(a,b)  show  a  region  of  the  thin  film  photovoltaic  cell  which  was  subjected  to 
bending  and  tension,  resulting  to  fragmentation. 


(a)  (b) 

Figure  D67.  Fragmentation  of  brittle  thin  film  layers  in  a  solar  cell  under  (a)  bending  and  (b)  tension. 
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We  quantified  the  film  residual  mean  and  gradient  stresses,  through  an  analysis  of  the  blister 
profiles  of  delamination  and  by  curvature  measurements  of  freestanding  film  strips. 

Delamination  was  induced  by  immersing  the  photovoltaic  cells  in  acetic  acid  to  etch  the 
polymeric  contact  and  release  the  TCO  and  p-n  junction.  Application  of  bending  to  the 
photovoltaic  films  while  immersed  in  acetic  acid  was  adequate  to  peel  the  top  ZnO  layer  from  Si. 
Thus,  the  interfacial  adhesion  between  the  Si  and  ZnO  layers  was  weak.  When  the  ZnO  layer 
was  removed,  the  Si  film  (still  attached  to  the  substrate)  formed  straight  de lamination  columns. 
Figure  D68,  due  to  the  residual  stresses  in  the  Si  layer.  The  mechanics  of  straight  delamination 
have  been  analyzed  before  by  modeling  the  thin  film  as  a  clamped -clamped  buckled  beam 
(Hutchinson,  1992),  which  leads  to  the  following  expressions  for  the  buckling  and  the  residual 
stresses  (Volinsky,  2003): 


ir^  E  /h^ 

'  12(1- v2)U/ 

(Dl) 

(D2) 

where  E  is  the  arithmetie  mean  of  the  elastic  moduli  of  the  two  thin  layers,  v  is  the  harmonic 
mean  of  Poisson’s  ratios,  h  is  the  film  thickness,  6  is  the  blister  height  and  b  is  the  blister  half¬ 
width.  From  Equations  (Dl)  and  (D2),  Ob  =  160±49  MPa  and  Or  =  466±1 1 0  MPa  (residual  mean 
stress). 

When  the  photovoltaic  films  were  simply 
immersed  in  acetic  acid  solution  without 
the  application  of  bending,  the  ZnO 
remained  attached  to  the  Si  surface. 

However,  the  ZnO  layer  delaminated  at 
edges  where  the  photovoltaic  films  were 
cut,  pointing  out  to  the  weak  ZnO-Si 
interface.  The  same  conditions  that  resulted 
in  straight  delaminations  of  only  the  Si 
layer  in  Figure  D68,  produced  the  “zig¬ 
zag”  delaminated  patterns  shown  in  Figure 
D69(b),  commonly  known  as  telephone 
cord  delaminations.  Figure  D69(a)  is  a 
close  up  of  the  ZnO  surface  showing  its 
grain  structure. 

In  order  to  calculate  the  residual  stresses 
responsible  for  the  formation  of  the  telephone  cord  de  laminations,  the  approximate  values  of  the 
buckling  and  the  residual  stresses  were  calculated  by  using  the  same  method  as  for  the  straight 
delaminations  (Hutchinson,  1992).  These  calculations  averaged  128±30  MPa  and  742±104  MPa 
for  the  buckling  and  the  residual  stress,  respectively. 

The  mean  residual  stress  computed  from  the  aforementioned  delaminations  is  based  on  the  Euler 
buckling  theory  and  it  does  not  account  for  the  effect  of  the  residual  stress  gradient  which  also 
contributes  to  the  formation  of  a  buckling  delamination.  Instead,  the  stress  gradient  can  be 
calculated  from  the  curvature  of  freestanding  Si/ZnO  strips.  While  the  Si/ZnO  bilayer  is  detached 
from  the  substrate,  the  cracks  on  the  ZnO  surface  either  propagate  into  the  Si  layer,  creating  a 
crack  through  the  entire  thickness  of  the  bi  layer,  or  in  the  weak  interface  causing  de  lamination. 


Figure  D68,  Straight  delamination  column 
propagating  to  the  left  of  this  image. 
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Figure  D70(a)  shows  a  Si  mono-layer  strip,  and  Figure  D70(b)  a  Si/ZnO  bi-layer  strip  that  are 
the  result  of  the  aforementioned  processes.  The  residual  stress  gradient  is  then  calculated  using 
curvature  measurements.  For  the  bi-layer  sample,  the  composite  moment  of  inertia  is  computed 
as  (Timoshenko,  1925); 


where  h  is  the  total  thickness,  w  is  the  width  of  the  beam  and  n  is  the  ratio  of  the  elastic  moduli. 
The  stress  gradient  in  the  Si  monolayer  was  found'equal  to  I71±12  MPa/gm,  while  for  the  bi- 
layer  it  was  computed  as  1 70±1 5  MPa/pm.  While  the  two  values  appear  to  be  the  same  it  should 
be  emphasized  that  they  correspond  to  1  pm  and  2  pm  thick  layers  and  as  a  result,  the  latter  were 
subject  to  a  larger  curvature, 
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Figure  D70.  Cross-sections  of  (a)  Si  and  (b)  ZnO/Si  bilayer  films. 


COMPOSITE  ANODES  AND  CATHODES  FOR  LI+  BASED  BATTERIES 


The  objective  of  this  task  was  to  investigate  the  failure  of  composite  anode  and  cathodes  for  thin 
film  batteries.  The  electrode  materials  were  obtained  in  the  Spring  and  in  the  Summer  of  201 1 
from  the  University  of  Colorado  at  Boulder,  Given  the  short  duration  of  this  project  (7  months 
plus  7  months  of  no  cost  extension)  we  focused  on  the  assembly  and  characterization  of  battery 
cells  and  in  situ  observations  of  the  surface  of  the  electrodes  during  lithiation/delithiation. 
Composite  Si  anode  and  LiCo02  cathode  samples  were  prepared  by  the  research  group  of  Prof. 
Se-Hee  Lee  (Associate  Professor,  Mechanical  Engineering)  at  the  University  of  Colorado  at 
Boulder.  They  were  composite  aggregates  spread  over  the  surface  of  a  Cu  or  an  A1  substrate,  as 
shown  in  Figures  D71(a,b).  The  anode  were  comprised  of  a  mixture  of  60%wt.  Si  particles, 
20%wt  acetylene  black  (AB),  and  20%wt  polyvinylidene  diflouride  (PVDF)  binder,  deposited  as 
a  slurry  on  a  (~19  pm)  Cu  foil  (total  thickness  ~50-55  pm).  The  exact  thickness  of  the  composite 
Si  layer  varied  significantly  between  different  samples  and  within  each  sample.  The  cathode 
samples  were  made  in  a  similar  manner  and  were  comprised  of  a  mixture  of  LiCo02  (83%wt), 
AB  (7.5%wt)  and  PVDF  (9.5%wt),  on  an  A1  foil.  Figures  D72,  D73  show  the  surface  details  of 
the  anode  films. 


Experiments  were  carried  out  to  measure  the  interlaminar  shear  strength  between  the  active  layer 
of  the  electrode  (Si+AB+PVDF,  henceforth  referred  to  as  the  composite  Si  layer)  and  the  Cu 
substrate.  The  tests  were  performed  using  a  Bose  ElectroForce®  3200.  In  total,  15  samples  were 
prepared,  shaped  as  long,  narrow  rectangular  pieces.  However,  only  7  samples  provided  reliable 
results  (as  concluded  from  the  load-displacement  curves).  The  sample  width  was  1 .06±0.047  mm 
and  the  length  was  -'lO  mm,  with  an  actual  gage  length  of-S  mm.  The  samples  were  bonded  to 
the  grips  through  the  Si  layer  on  one  side  using  super-glue  (Elmer’s  Fix-all  adhesive),  while  on 
the  other  side  the  entire  film  was  bonded  using  a  12-hour  cure  time  epoxy.  When  subjected  to 
tension,  the  end  bonded  through  the  Si  layer  delaminated  from  the  Cu  layer.  The  force  at  failure 
was  used  to  calculate  the  shear  strength  of  the  Cu  -  composite  Si  interface.  While  the  stress  at  the 
point  of  failure  might  also  include  the  tensile  strength  of  the  Si  layer  alone,  this  was  neglected  in 
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our  calculation  as  the  Si  layer  was  observed  to  have  low  strength.  The  results  are  tabulated  in 
Table  D 13. 


I'  igure  D72.  Fexture  of  Si  anode  materia!  (Si  layer  is  visible). 


Figure  D73-  A.  Si  particles,  B.  acetylene  black  (AB)  and  C,  PVDF  binder. 
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Table  D13.  Shear  strength  of  composite  Si  —  Cu  interface. 


Sample 

Width  (mm) 

Strength  (MPa) 

1 

1.04 

164.58 

2 

1.17 

68.32 

3 

1.09 

75.42 

4 

1.08 

97.92 

5 

1.06 

102.24 

6 

0.95 

183.08 

7 

1.02 

172.86 

Average 

1.06 

123+48 

It  should  be  noted  that  the  interfacial  strength  values  in  Table  D13  are  quite  high  when  compared 
to  common  engineering  material  systems.  This  may  be  the  result  of  the  unaccounted  tensile  force 
in  the  composite  Si  layer,  and  uncertainties  in  the  calculation  of  the  effective  area  subjected  to 
interfacial  shear  stress.  Nevertheless  the  particular  experiment  provides  the  means  to  test  such 
interfaces  and  obtain  quantitative  results  about  the  quality  of  bonding  between  the  electrode 
(PVDF  composite)  and  the  current  collector  (metal). 

In  order  to  measure  strain  in  the  composite  Si  layer  during  lithiation  and  de-lithiation,  the 
electrodes  were  subjected  to  both  conditions  under  different  electrochemical  conditions  and 
charge/discharge  regimes.  Subsequently,  the  surface  topography  was  compared  before  and  after 
each  test  using  a  confocal  laser  microscope.  Comparative  images  from  5  different  samples  are 
presented  in  Figures  D74-D79.  A  variety  of  failure  modes  were  observed,  varying  from  matrix 
surface  cracking  taking  place  during  the  first  lithiation,  to  the  formation  of  a  form  of  SEI  which 
covered  almost  the  entire  electrode. 

The  average  interlaminar  shear  strength  between  the  electrode  layer  and  the  substrate  was  found 
to  be  of  the  order  of  120  MPa  for  the  composite  Si  anode  samples.  This  implies  that  the  metal- 
polymer  composite  layers  were  bonded  well,  thus  loss  in  efficiency  or  failure  of  the  cell  or 
electrode  due  to  debonding  of  these  layers  is  unlikely. 

The  composite  LiCo2  cathode  samples  on  the  other  hand,  could  not  be  tested  for  interlaminar 
shear  strength.  The  specimens  failed  under  tension  before  the  two  layers  could  separate  or 
debond.  The  peak  stresses  at  failure  were  much  lower  than  the  interlaminar  shear  strength  of  the 
composite  Si  samples,  due  to  the  low  strength  of  the  A I  substrate. 
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Figure  D74.  Surface  of  Si  composite  electrode  sample  1  (a)  before  and  (b)  after  lithiation  at  50  x 
magnification  (lithiated  under  potentiostatic  conditions,  with  an  applied  voltage  of  04  V). 
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Figure  D75-  Surface 

after  de-lithiation,  at  lOOx  magnification  (lithiated  for  10  minutes  under  potentiostatic  conditions,  with  an 
applied  voltage  of  OJ  V,  then  fully  de-lithiated  at  1.5  V). 
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Figure  D76.  Surface  of  Si  composite  electrode  sample  3  (a)  before  lithiation,  (b)  after  lithiation  and  (c) 
after  de-lithiation,  at  50x  magnification  (lithiated  for  20  minutes  under  potentiostatic  conditions,  with  an 
applied  voltage  of  0.1  V,  then  fully  de-lithiated  at  1.5  V). 
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Figure  D77.  Surface  of  Si  composite  electTode  sample  3  (a)  before  lithiation,  (b)  after  lithiation  and  (c) 
after  de-lithiation,  at  lOOx  magnification  (lithiated  for  20  minutes  under  potentiostatic  conditions,  with  an 
applied  voltage  of  0. 1  V,  then  fully  de-lithiated  at  1 .5  V), 
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(C)  (d) 

Figure  D78.  Two  regions  on  the  surface  of  Si  composite  electrode  sample  4  {a,b)  before  and  (c,d)  after  lithiation  at  50x  magnification; 
(lithiated  under  galvanostatic  conditions  at  very  low  charging  rate,  --C/300). 


(b) 

Figure  D79.  Surface  of  composite  Si  electrode  sample  5  (a)  before  and  (b)  after  lithiation  at  SO’* 
magnification  (lithiated  under  gaivanostatic  conditions  at  2C/3  charging  rate). 


The  surfaces  of  the  Si  anode  and  the  LiCo02  cathode  supplied  by  the  University  of  Colorado  was 
significantly  rough.  In  the  case  of  Si  based  anodes,  the  peak-to- valley  surface  roughness  was  as 
high  as  70  pm  (ignoring  the  occasional  spikes  often  seen  in  confocal  microscopy  images)  as 
shown  in  Figure  D80.  Whereas,  in  the  LiCo02  cathode  samples,  the  peak-to- valley  surface 
roughness  was  of  the  order  of  45  pm,  as  shown  in  Figure  D8 1 .  Due  to  the  large  surface  features, 
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the  use  of  an  AFM  to  perform  surface  measurements,  as  was  initially  proposed,  was  ruled  out. 
Instead,  the  electrode  surfaces  were  compared  using  confocal  laser  microscopy.  During  the  initial 
electrochemical  experiments,  no  optically  discernible  changes  occurred  on  the  electrodes’ 
surfaces  (both  Si  and  LiCo02).  This  was  attributed  to  incomplete  tithiation/delithiation  during 
the  charge/discharge  cycle.  Experiments  were  then  conducted  using  constant  current  and 
constant  voltage.  While  the  surface  features  changed  drastically  during  lithiation,  very  limited 
changes  took  place  upon  de-lithiation.  Furthermore,  the  Si  particles  rarely  expanded  or 
contracted,  whereas  the  PVDF  binder  changed  dramatically,  as  shown  in  Figures  D74-D79, 
In-situ  imaging  of  the  Si  electrode’s  surface  revealed  a  dynamically  forming  and  expanding 
surface  layer  which  interfered  with  the  optical  analysis  of  the  sample  surface.  After  the  charge- 
discharge  cycle  was  over,  this  layer  was  deposited  onto  the  surface  accounting  for  the  large 
changes  in  the  surface  morphology.  The  charging  rate  determined  the  extent  of  such  changes  on 
the  electrode’s  surface,  and  barely  any  change  was  seen  for  very  low  charging  rates  (~C/100  or 
less).  This  is  expected,  since  the  SEI  would  have  sufficient  time  at  low  charging  rates  to  dissolve 
back  into  the  electrolyte  without  a  residue. 

An  in-depth  analysis  of  the  relationship  and  interdependence  between  the  electrochemical 
performance  and  mechanical  properties  of  the  Si  anode  and  LiCo02  cathode  materials  would 
require  well-prepared  samples  that  would  allow  individual  particles  to  be  observed  by  time 
resolved  microscopy  techniques. 


Figure  D80.  Surface  roughness  of  composite  Si  anode. 
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Task  E:  Design  of  Energy  Storage  System 

This  task  includes  the  subtasks  of  designing  Li-ion  battery  cathodes  and  super  capacitors 
which  have  higher  performance  and  durable  in  airborne  use  conditions. 

E-l:  Design  of  new  electrode  materials  (Cao,  UW) 

Li  ion  battery  research  can  be  divided  into  two  parallel  efforts:  (1)  developing  coherent 
carbon  cryogel  -  vanadium  pentoxide  nanocomposite  electrodes  for  both  enhanced 
specific  power  and  specific  energy  and  (2)  new  chemical  compositions  and  new 
nanostructures.  Three  approaches  have  been  taken  to  incorporate  vanadium  oxide  to 
mesoporous  carbon  cryogels:  (a)  electrochemical  deposition,  (b)  water  electrolysis 
induced  deposition  of  vanadium  oxides  due  to  the  pH  change  at  the  vicinity  of  the  carbon 
surface,  and  (c)  sol  electrophoretic  deposition.  Such  coherent  nanocomposites  would 
offer  excellent  energy  storage  performance;  however,  so  far  no  sufficient  amount  of 
vanadium  oxide  has  been  successfully  incorporated  into  carbon  cryogels  yet.  On  a 
parallel  effort,  MnOx  nanowall  arrays  and  Ti02  nanotube  arrays  have  been  successfully 
fabricated  and  tested  for  Li  ion  intercalation  and  demonstrated  some  promising  results. 
Figure  El  shows  the  SEM  images  of  both  MnOx  nanowall  arrays  fabricated  by  water 
electrolysis  induced  deposition,  and  Ti02  nanotube  arrays  by  means  of  acidic  anodization. 


E-2:  Carbon  Cryogel  supercapacitors  (Cao,  UW) 

The  research  in  the  past  year  has  been  focused  on  the  understanding  of  the  determining 
factors  on  the  capacitance:  microstructure  or  bonding.  Our  prior  research  has 
demonstrated  that  the  capacitance  varies  significantly  with  the  processing  conditions. 
Different  processing  conditions  result  in  different  meso  and  microporous  structures; 
however,  the  capacitance  does  not  find  linear  relationship  with  porous  structure.  Our 
recent  X-ray  Raman  spectroscopy  experiments  have  ruled  out  the  possibility  that  the 
different  processing  conditions  may  result  in  different  chemical  bonds  in  carbon  cryogels 
(see  Figure  E2).  It  is  concluded  that  the  capacitance  in  carbon  cryogel  is  solely  dependent 
on  meso  and  microporous  structure;  however,  it  is  controlled  by  a  combination  of  surface 
area,  pore  size  and  pore  volume. 
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Figure  E2,  X-ray  Rainan  scattering 
spectra  of  carbon  cryogels  fabricated 
under  various  processing  conditions, 
showing  all  the  samples  have 
identical  local  structure  and  chemical 
bonding. 


E-3:  Design  of  high  performance  electrolytes  (Xu,  UW) 

Accomplishments: 

12,5  mS/cm  ion  conductivity  obtained  in  solution  state 

^3  mS/cm  ion  conductivity  obtained  in  gel  state 

2.3  mS/cm  ion  conductivity  obtained  in  solid  state.  See  Fig.  E3 
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Figure  E3  Ionic  conductivity  of  different  compositions  and  states 
(measured  at  room  temperature). 


E-4  Microstructural  Studies  of  Li  Insertioii(Lee,  CU) 

L  Objectives 

Traditional  LI«ion  batteries  employ  carbonaceous  anodes,  but  projected  performance 
targets  for  next-generation  Li-ion  batteries  require  new,  higher  capacity  materials  that  can 
be  electrochemical  I  y  cycled  in  a  stable  manner.  The  most  attractive  candidate  to  replace 
carbon  is  silicon;  it  has  the  highest  known  capacity,  is  relatively  low-cost,  and  is 
physically  abundant.  At  present,  the  drawback  with  silicon  is  that  a  volume  expansion  on 
the  order  of  ca,  300%  occurs  upon  Li  insertion  that  leads  to  rapid  capacity  fade  during 
cycling.  Creative  approaches  to  realize  highly-reliable  electrodes  are  being  pursued  with 
nanometer-scale  materials  and  geometric  architectures;  impressive  results  have  been 
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obtained,  but  a  fundamental  understanding  is  still  lacking.  To  this  end,  this  project  has 
focused  on  understanding  the  complex  phenomena  of  Li  insertion  and  extraction  from  Si, 
and  developing  novel  nanostructured  Si-based  anodes  which  can  address  these  critical 
technical  issues  with  Si  anodes. 

2.  Brief  summary  of  the  CU  team  progress  toward  the  above  objectives 
We  have  reported  the  direct  observation  of  microstructurai  changes  of  Li^Si  electrode 
with  lithium  insertion.  HRTEM  experiments  confirm  that  lithiated  amorphous  silicon 
forms  a  shell  around  a  core  made  up  of  the  unlithiated  silicon  and  that  fully  lithiated 
silicon  contains  a  large  number  of  pores  of  which  concentration  increases  toward  the 
center  of  the  particle.  Chemomechanical  modeling  has  been  employed  in  order  to  explain 
this  mechanical  degradation  resulting  from  stresses  in  the  Li^Si  particles  with  lithium 
insertion.  Beeause  lithiation-induced  volume  expansion  and  pulverization  are  the  key 
mechanical  effects  that  plague  the  performance  and  lifetime  of  high-capacity  Si  anodes  in 
lithium-ion  batteries,  our  observations  and  chemomechanical  simulation  provide 
important  mechanistic  insight  for  the  design  of  advanced  battery  materials. 


Figure  E4.  TEM  images  of  Li2Si.  (a)  BE  image  of  L^Si.  (b)  HR  and  corresponding  EFT 
of  area  A  indicates  that  amorphous,  (c)  HR  and  corresponding  EFT  of  area  B  indicates 
that  Si  crystalline  phase  along  the  [112]  direction. 

Figure  E4a  shows  the  bright  field  (BF)  TEM  image  of  the  one  single  particle  Li2Si 
sample  which  has  been  known  to  exhibit  a  two-phase  region.  Energy  dispersive  X-Ray 
spectroscopy  (EDS)  measurements  shown  in  inset  of  Figure  E4a  reveal  that  samples  are 
mostly  composed  of  Si,  with  a  small  amount  of  Cu  detected.  No  signal  from  oxygen  is 
detected,  which  implies  that  the  TEM  sample  is  never  exposed  to  air.  It  is  observed  that 
the  majority  of  the  surface  (which  contacts  the  Pt)  is  amorphous  phase  while  the  inside 
area  of  the  particle  is  filled  with  crystalline.  Fig.  E4b  shows  HR  TEM  image  and 
corresponding  FFT  of  area  A,  clearly  indicating  the  area  A  is  amorphous  region.  BF 
image,  HR  images  and  FFTs  in  Figure  lb  and  Figure  Ic  support  that  each  area  is  formed 
with  amorphous  and  crystalline  Si  phase  (along  [112]  zone  axis),  respectively.  This 
confirms  that  a  two-phase  region  is  formed  with  amorphous  Li^Si  at  the  surface  and 
crystalline  Si  phase  in  the  core  of  particle  as  previously  proposed.  As  shown  in  Figure  la, 
amorphization  due  to  Li  insertion  in  Si  begins  from  the  surface.  While  it  is  widely 
accepted  that  amorphous  Li^Si  initially  forms  at  the  surface  with  crystalline  Si  being 
maintained  in  the  core  of  particles,  there  has  never  been  precise  observation  of  this 
phenomenon  using  TEM.  Here,  a  LixSi  particle  maintaining  a  crystalline  core  with 
amorphization  at  the  surface  during  Li  insertion  is  confirmed.  This  observation  also  infers 
that  Li  concentration  in  the  Si  particle  is  surface  dominant  which  explains  mechanical 
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failure  of  Si  anodes  upon  Li  insertion.  Thermodynamically  stable  Li-Si  crystalline  phases 
such  as  Lii2Si7,  Li7Si3  and  LioSL  are  not  shown  in  this  region  because  the  energy  barrier 
for  nucleation  which  is  combined  with  surface  and  volume  energy  of  nucleus,  is  higher 
when  it  is  about  to  nucleate  and  thus  makes  meta-stable  amorphous  Li^Si , 

3.  Accomplish  men  ts/New  findings 

1)  Direct  observation  of  microstructural  transformations  of  lithiated  Si  anodes  in  all- 
solid-state  Li  half  cells. 

2)  It  has  been  shown  that  crystalline  Si  undergoes  amorphization  during  lithiation  and 
forms  a  crystalline  phase  after  full  discharge,  known  as  LiisSL.  Two-phase  regions  have 
been  observed  by  HRTEM  showing  amorphous  LixSi  at  the  surface  and  crystalline  Si 
phase  at  the  core  of  particles  for  the  first  time. 

3)  HRTEM  observation  of  fully  discharged  Si  anode  samples  revealed  that  particles 
consisting  of  high  concentrations  of  pores  in  the  Si  particle  are  responsible  for  cracks  and 
pulverization  of  Si  during  Li  insertion. 

4)  The  chemomechanical  modeling  has  been  employed  to  explain  the  generated  stresses 
during  Li  insertion  in  relation  to  observed  porosity  pattern  in  detailed  TEM  study. 

E-5  Electrospinning  fabrication  and  electrochemical  properties  of  LiCoO;  and 
Li(Nii/3Coi/3Mn  1/3)02  for  cathodic  material  of  Li-ion  battery  (Taya) 

we  have  successfully  fabricated  single-wall  carbon  nanotubes  embedded  LiCo02 
nanofibers  by  electrospinning  technique  from  a  precursor  solution  containing  lithium 
acctatc/cobalt  acctatc/PVP  (Polyvinylpyrrolidone)  with  different  concentration  of  single 
wall  carbon  nanotubes  (SWCNT).  XRD  (x-ray  diffaction).  Fig.  E6,  and  SEM  (scanning 
electron  microcope) ,  see  Fig.  ESwere  performed  to  investigate  the  phase  and 


Fig.  E5  SEM 
images  of 
LiCo02- 
SWCNT  with 
(a)(c)  before 
calcine,  (b)(d) 
after  calcine 
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Fig.ES  shows  the  SEM  photographs  of  the  collected  LiCo02-S  WCNT  nanofibers  by 
electrospinning  followed  by  calcination  at  500°C  for  30m in.  All  the  nanofibers  are 
randomly  distributed  on  the  substrate.  A  closer  look  of  pre-calcined  nanofiber  shown  in 
Fig.E5(c)  reveals  that  carbon  nanofibers  could  be  observed  from  a  breaking  point  of 
composite  fibers.  Whereas  Fig.E5(d)  shows  that  nano-sized  LiCo02  particles  was  formed 
an  junction  points  from  which  other  nanofibers  are  branched 


Fig.  E6  XRD  results  for 
LiCoO:  nanofibers 
deposited  on  A1  substrate 
after  calcination  at  500“C 
for  30  min 


Fig.E6  shows  the  XRD  patterns  for  LiCo02  nanofibers  with  different  Li'^/Co^'^/SWCNT 
starting  ratio  deposited  on  Al  substrate  after  calcinations  at  500°C  for  30m in.  As  shown 
in  Fig.E6,  the  Al  substrate  has  very  strong  background  peak  intensity  in  all  samples,  the 
three  peaks  marked  with  star  indicated  the  three  strongest  peaks  of  the  hexagonal  LiCo02. 
Our  previous  report  based  on  thick  layer  of  LiCo02  electrospun  nanofibers  indicates  that 
LiCoOa  fabricated  by  similar  method  exhibited  a  pure  crystalline  LiCo02  after 
calcination  at  500®C  for  30min.  In  our  previous  work,  the  impurity  concentration  in 
SWCNT  (graphite,  iron  oxide)  was  too  low,  there  was  no  peak  in  XRD  corresponding 
these  contaminates. 


Fig.E7  XRD  results  for 
Li(Ni|/3Coi/3Mn  1/3)02  nanofibers 
deposited  on  Al  substrate  after 
calcination  at  different 
temperatures  for  6  hours 
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Fig.E7  shows  the  XRD  patterns  for  Li(Nii/3Coi/3Mni/3)02  nanofibers  after  calcinations  at 
725  to  SOO^C  for  6  hours.  The  patterns  of  three  strongest  peaks  coincide  very  well  with 
LiCo02.  In  fact,  the  peak  pattern  of  Li(Ni|/3Co i/sMn  1/3)02  was  very  similar  to  LiCo02 
and  LiNi02,  except  that  the  peak  position  was  in  between  of  LiCo02  and  LiNi02due  to 
its  lattice  parameter.  Furthermore  CV  (cyclic  voltammetry)  and  charge-discharge 
experiments  were  applied  to  characterize  the  electrochemical  properties  of  composite 
nanofibers.  To  improve  the  capacity  of  cathode  materials,  more  recent  work  of 
Li(Ni!/3Coi/3Mn  1/3)02  nanofiber  based  on  electrospinning  method  has  been  done  and  the 
capacity  has  been  dramatically  improved.  It  is  noteworthy  that  this  nanostructured 
cathode  offers  high  charge-discharge  capacity  even  under  high  charge-discharge  rate 
compared  with  conventional  micro  powder  or  film  cathodes,  therefore  SWCNT-LiCo02 
composite  nanofiber  as  well  as  Li(Nii/3Coi/3Mn  1/3)02  naon fiber  cathode  may  become 
promising  candidates  for  high  energy  density  cathodic  materials  of  lithium-ion  batteries. 

The  morphology  of  Li(Nii /3Coi/3Mn  1/3)02  nanofiber  by  SEM  was  very  similar  to  LiCo02- 
SWCNT  nanofiber.  The  Due  to  resolution  limitation  of  SEM,  identification  of  embedded 
SWCNT  will  be  examined  later  by  TEM.  The  average  diameter  of  these  fibers  was  found 
to  be  200-500nm  after  calcination.  The  diameter  of  the  nanofibers  would  depend  on  the 
PVP  concentration  as  well  as  electrospin  parameters. 


Fig.ES  TEM  image  of  pure  SWCNT  (a)  and  SWCNT  embedded  within  LiCo02  (b) 

Fig.ESshows  the  TEM  image  of  pure  SWCNT  and  a  bundle  of  SWCNT  among  LiCo02 
crystals.  The  SWCNT  is  characterized  by  rough  tube  wall  morphology  with  low  contrast, 
compared  with  higher  contrast,  sharp  crystal  planes  of  LiCo02.  However,  during  the 
observation,  the  concentration  of  SWCNT  was  not  as  high  as  proposed  according  to 
material  starting  ratio,  which  means  that  only  a  small  amount  of  SWCNT  in 
electrospinning  precursor  solution  was  embedded  into  LiCo02  nanofibers.  This  problem 
could  be  attributed  to  highly  aggregated  bundles  of  SWCNT,  which  was  difficult  to 
isolate  even  by  ultrasonication. 


108 


In  order  to  improve  the  capacity  of  cathode  materials,  the  novel  material  with  higher 
capacity  such  as  Li(Ni]/3Coi/3Mni/3)02  was  very  promising  candidate  for  lithium-ion 
battery.  The  discharge  capacities  at  2C  of  Li(Nii/3Coi/3Mni/3)02  nanofiber  with  5wt% 
MWCNT  under  different  calcinations  temperatures  were  tested  and  the  best  result  with 
calcination  temperature  at  SOOX  are  shown  in  Fig.9. 


Fig.E9  Discharge  capacity  of  the 
cathode  based  on 
Li(Ni|/3Coi/3Mni/3)02“M  WCNT 
nanofibers  under  current  rate  2C 


With  the  first  discharge  capacity  of  almost 
l90mAh/g,  it  clearly  indicated  that  Li(Ni[/3Coi/3Mni/3)02-MWCNT  nanofibers  show 
dramatically  capacity  improvement  compared  with  LiCo02-SWCNT  nanofiber.  Even 
with  some  degradation  after  5  cycles,  the  discharge  capacity  is  still  higher  than 
150mAh/g,  indicating  that  Li(Ni|/3Coi/3Mni/3)02-MWCNT  nanofibers  is  a  very 
promising  candidate  for  high  energy  density  lithium-ion  batteries.  In  future,  more  work 
should  be  focused  on  how  to  improve  the  discharge  stability  for  better  cycleability. 


Progress  Report:  Task  F  Circuit  Design  for  EHSS  (Kuga) 

Objectives  and  Proposed  tasks  during  the  1st  year 

During  the  first  year,  we  have  conducted  the  following  projects  under  Task  F:  (1) 
Development  of  HF-VHF  MMD  antennas  for  a  small  UAV,  (2)  development  of  the 
control  circuits  for  ECW  and  other  conductive  polymer  devices,  and  (3)  Integration  of 
energy  harvesting  circuits  and  study  of  energy  efficiency  using  the  testbed. 

(1)  Development  of  FIF-VHF  MMD  antennas  for  a  small  UAV 

The  unlicensed  HF  (27  MHz)  or  VHF  (72  MHz)  bands  can  be  utilized  for  long  distance 
communications  or  for  beyond  light-of-sight  (BLOS)  communications.  BEOS 
communication  is  often  called  upon  when  hills  and  building  are  in  the  landscape.  To 
attach  an  HF-VHF  antenna  to  a  small  UAV  (the  wing  span  is  less  than  2  m),  we  have 
studied  the  modified  meanderline  dipole  (MMD)  antennas.  The  length  of  the  existing 
dipole  antennas  was  not  suitable  for  small  UAVs.  For  size  efficiency,  we  designed  the 
MMW  antenna  such  that  this  antenna  would  fold  into  the  top  and  bottom  surfaces  of  the 
main  wing.  This  folding  would  allow  it  to  be  embedded  into  the  main  wings.  To  optimize 
the  antenna’s  performance,  the  design  must  include  realistic  air-frame  materials  such  as 
composite  wings. 
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Fig.F'l  An  example  of  a  small  UAV  (bottom)  Fig,  F-2  The  120“  MMD  antenna  on 

and  layout  of  the  MMD  antenna  (top).  dielectric  material  { £•,  =  4 ) 


Fig.  F-1  shows  an  example  of  a  small  UAV  as  well  as  the  layout  of  an  MMD  antenna. 
The  width  of  the  trace  is  5  mm  and  the  thickness  is  0,5  mm  of  the  MMD  antenna.  The 
total  length  of  the  MMD  antenna  is  1200  mm  while  the  total  width  is  140  mm.  The  folded 
length  in  the  MMD  antenna  is  20  mm.  This  antenna  is  designed  at  50  MHz  and  the 
simulated  antenna  responses  are  shown  in  Fig.  F-3. 


KrequfFvcytMtfi) 


Fig.  F-3  (left)  Simulated  Sn  of  the  MMD  antenna  as  shown  in  Fig.  F-1.  (right)  Simulated 
E-plane  radiation  patterns  of  the  MMD  antenna  at  50  MHz  as  shown  in  Fig.  F-1 


(2)  Development  of  the  control  circuits  for  ECW  and  other  conductive  polymer  devices 
We  have  developed  a  driving  circuit  for  ECW  which  can  accommodate  solar  cells. 


DPDT  switch 
reverses  polarity 


(3)  Integration  of  energy  harvesting  circuits  and  study  of  energy  efficiency 

To  test  the  performance  of  the  energy  harvesting  system,  we  have  developed  the  testbed 

using  existing  devices. 


Fig.  F-5  Testbed  for  energy 
harvesting  system 


Task  G.  Multimode  Energy  Harvesting  for  ISR/MAV  Missions  (Inman) 

Goals  Our  subcontracted  task  is  to  examine  the  use  of  piezoelectric  materials  for 
harvesting  ambient  energy.  Our  goals  this  period  included  examining  the  state  of  the  art 
in  piezoceramic  based  harvesting  and  discovering  the  best  predictive  models  for 
harvesters.  The  follow  on  goal  is  to  form  these  predictive  models  into  useful  design  rules 
and  to  transition  these  to  AF  applications.  An  additional  goal  is  to  characterize  ambient 
energy  sources  for  mechanical  harvesting. 

Main  findings 

The  main  discoveries  and  inventions  are  summarized  as  follows: 

Piezoelectric  cantilever  harvesters  as  currently  existing  at  the  beginning  of  this  MURl 
where  modeled  in  correctly.  So  in  the  first  year  of  the  MURI  new,  higher  fidelity  models 
where  derived  and  experimentally  clarified  which  corrected  the  existing  models. 


Figure  G-1  Plots  showing  the  error  of  over  prediction  of  early  energy  harvesting  models 
and  the  corrective  mode!  developed  under  the  MURI  which  agrees  with  experimental 
data.  The  red  dashed  line  is  the  previous  theory  based  on  uncoupled  modeling  and  the 
green  line  is  the  improved  model,  which  agrees  with  experiments* 

During  the  second  year  we  addressed  structural  integration  of  harvesting  materials.  The 
idea  developed  was  that  each  wing  has  a  span  The  spar  has  a  designed  volume,  mass  and 
stiffness.  We  kept  these  three  parameters  fixed  while  we  adjusted  the  amount  of 
piezoelectric  material  (PZT)  in  the  spar  to  optimize  the  energy  we  can  extract  as  the  spar 
vibrations.  Our  initial  flight  test  resulted  in  50  micro  watts  but  this  was  not  optimized. 
Our  colleagues  at  Colorado  the  used  this  model  to  optimize  the  distribution  of  PZT 
throughout  the  spar.  We  provided  them  with  a  mechanics  model  that  they  used  their  code 
on.  We  also  performed  some  flight  tests  in  order  to  determine  how  much  vibration 
energy  was  available  in  a  UAV.  We  also  discovered  the  concept  of  strain  nodes  which 
leads  to  cancellation  of  power  out  and  restricts  the  design  space. 


Figure  G-2  Piezoelectric  harvesting  devices  integrated  in  the  body  of  a  small  UAV  as  a 
cantilever  (top  photo)  and  layered  into  the  spar  (bottom  picture). 
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During  the  third  year  we  investigated  the  concept  of  using  '"L-shaped”  piezoceramic 
transducers  to  try  and  improve  the  amount  of  energy  obtained  from  a  given  input 
vibrations.  We  also  investigated  using  both  piezoceramic  transduction  and 
electromagnetic  transduction  in  the  same  device.  The  PZT  produces  high  voltage  but  low 
current  where  as  the  electromagnetic  transduction  produces  high  current  but  low  voltage 
at  small  levels  of  vibration.  Hence  the  combination  provides  a  high  enough  voltage  to 
charge  a  battery  and  a  high  enough  current  to  increase  the  charge  time. 


Figure  G-3  On  the  right  is  the  L-shaped  configuration  used  to  introduce  nonlinear 
behavior  and  on  the  left  is  the  combined  piezoelectric  and  electromagnetic  harvester. 

During  the  forth  year  we  integrated  thin  film  batteries  into  our  wing  spar  along  with  PZT 
for  harvesting  and  developed  the  concept  of  a  "'seif-charging'’  structure.  We  looked  and 
strength  tests  and  characterized  the  multifunctionality  of  the  self  charging  spar  in  terms  of 
weight  strength  and  performance. 

During  the  last  year  several  new  discoveries  where  made  and  several  new  devices 
designed  and  fabricated.  First  it  was  discover  through  a  series  of  experiments  and 
analytical  derivations  that  the  amount  of  energy  harvested  from  mechanical  motion  could 
be  increased  an  order  of  magnitude  by  introducing  nonlinear  behavior  into  a  cantilever¬ 
harvesting  configuration.  This  was  further  expanded  to  structural  systems  by  using  by¬ 
stable  composite  materials.  Last  we  derived  the  mathematical  models  for  and 
experimentally  validated  the  harvesting  of  energy  through  aeroelastic  phenomena. 
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Figure  G-4  Experiment  to  validate  the  theory  that  energy  harvested  from  light  through 
wing  vibration  produces  usable  electrical  energy  and  increases  the  flutter  speed. 

Our  second  group  of  achievements  consisted  of  designing,  modeling  and  fabricating  self 
charging  structures  which  gather  energy  from  both  solar  arrays  and  from  mechanical 
vibrations.  This  system  level  approach  to  multifunctional  structures  forms  a  multi-mode 
harvesting  and  storage  device  consisting  of 
] .  A  piezoceramic  energy  harvesting  device 

2.  A  thin  film  solar  panel 

3.  A  thin  film  battery 

4.  Electronics  to  combine  both  forms  of  energy  into  a  single  battery 

The  result  was  fashioned  into  a  wing  spar  of  a  model  airplane  and  flown.  The  spar  itself 
was  tested  in  the  laboratory  illustrating  the  proof  of  concept.  Strength  analysis  was 
performed  along  with  a  measure  of  multifunctionality  based  on  the  flight  formulation  put 
forth  by  Thomas,  et.  al. 


Piezocerainic 
^  in  Kapton  film 

Tbin-film 

ninum  ^ 

tructure 


Aluminum 
•>  substructure 

Figure  G-5  Integrated  self  charging  with  spar  including  harvesting  element  as  well  as 

storage  element. 
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Chapter  5.  Interactions  with  DoD  lab  and  industry 

a)  DoD  labs 

Taya  of  UW  visited  J.  Thomas  of  Naval  Research  Laboratory  to  discuss  the 
thermoelectric  generator  mounted  on  a  small-sized  UAV. 

CU  team  interate  with  G.  Reich,  Adaptive  Structures,  Vehicles  Directorate,  AFRL, 
Dayton,  OH. 

b)  DoD  industry 

Taya  of  UW  discussed  with  Bradley  J  Mitchell, Boeing  Company,  regarding  possible 
technology  transfer  of  the  thermoelectric  generator  to  Boeing  commercial  airplanes. 

Taya  gave  a  talk  at  GM  headquarters  on  Energy-harvesting  system  including  airborne 
applications 

Chapter  6.  Honors/awards 

Keynote  Address:  K,  Maute  and  M.L.  Dunn,  2012  Society  of  Engineering  Science 
Annual  Meeting,  Atlanta,  GA, 

T.  Hahn,  World  Fellow,  International  Committee  on  Composites,  July  2009 
T.  Hahn,  Scala  Lecturer,  International  Committee  on  Composites,  August  201 1 
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